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Preface 


Thii  document  repreienti  a  great  deal.  To  aome  extent,  it  repreienti  the  perional 
accomplishment  of  individual  research.  However,  a  more  fundamental  fact  is  apparent. 
This  thesis  represents  the  educational  opportunities  and  superb  research  facilites  which 
have  been  available  to  me  at  AFIT.  Closer  scrutiny  reveals  a  still  more  basic  fact.  This 
thesis  represents  the  outgrowth  of  my  affiliation  with  the  highly  professional  AFIT  faculty 
and  with  my  classmates  here.  Ultimately,  however,  this  document  represents  the  years  of 
dedicated  research  which  preceded  it. 

In  reviewing  the  draft  of  this  document,  1  noted  that  the  references  (which  have 
proven  to  be  vitally  important  to  me)  collectively  span  twenty  years.  This  is  a  particularly 
important  fact  in  my  estimation.  First  it  demonstrates  that  some  “timeless”  works  have 
been  produced  which  underpin  all  research  of  this  type.  Additionally,  it  points  out  how 
strongly  current  developments  can  affect  the  direction  and  emphasis  of  those  engaged  In 
independent  research  such  as  this. 

I  do  not  claim  (to  borrow  the  immortal  Words  of  Newton)  that  I  have  “seen  farther 
than  other  men”,  but  in  my  attempt  to  see  beyond  my  own  horizons,  I  have  certainly  “stood 
on  the  shoulders  of  gisuits.”  Those  whose  influence  has  contributed  to  the  work  represented 
here  may  never  know  it,  but  I  wish  to  thank  them  profoundly  for  their  effort.  My  respect 
for  their  vision,  my  appreciation  of  their  dedication,  and  my  outright  amazement  at  their 
accomplishments  have  all  grown  without  bound  during  my  tour  at  AFIT. 

I  wish  to  extend  special  thanks  to  my  thesis  committee  members,  Dr  Peter  Maybeck 
and  Col  Stan  Lewantowice,  and  to  my  thesis  advisor,  Capt  Randy  Paschall,  for  their 
contributions  to  this  research.  Their  collective  experience  and  advice  proved  to  be  an 
extraordinarily  rich  resource  which  never  diminished  in  quality  or  quantity. 

I  also  want  to  thank  Captains  Britt  Snodgrass  and  Jack  Taylor.  Their  selfless  assis¬ 
tance  throughout  this  effort  proved  invaluable  in  helping  me  to  “And  the  handles.” 

In  a  somewhat  sentimental  tradition  of  my  own,  l  want  to  include  the  names  of  my 
classmates  here.  Their  contributions  have  come  in  many  forms  -  not  always  academic.  Let 
friends  be  friends  forever. 


ii 


Tom  Cox,  Ken  Crosby,  Roger  Evans,  Tom  Flynn,  Maurice  Marlin,  Russ  Miller, 
Rob  Moyle,  Barbara  Niblett,  Paul  Rohi,  Dan  Sims,  Jack  Taylor,  Mark  VonBokern,  Dave 
Wheaton 

I  will  remember  all  of  you,  always,  and  I  wish  you  the  beat. 


Though  the  contributions  of  those  mentioned  above  were  many,  and  the  magnitude 
of  my  appreciation  la  great,  there  la  one  who  has  contributed  far  more  than  they,  and  to 
whom  I  offer  my  deepest  gratitude.  I  will  not  compare  her  gentle,  understanding  nature 
to  others,  for  none  can  compare.  She  will  always  be  the  light  and  happiness  in  my  life. 
1  dedicate  this  work  to  my  devoted  wife,  Valerie,  and  to  our  unborn  child  she  bears  so 
proudly.  I  love  them  both  beyond  measure. 


Richard  D.  Stacey 
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Abstract 


To  quantify  the  performance  abilities  of  existing  or  proposed  navigation  systems,  the 
U.S.  Air  Force  hat  for  the  last  several  yean  compared  the  performance  of  the  lyttem  under 
test  to  the  performance  of  a  baseline  navigation  system  known  at  the  Completely  Integrated 
Reference  Instrumentation  System  (CIRIS).  C1RIS  obtains  a  highly  accurate  navigation 
solution  by  combining  Ae  output  from  three  major  subsystems:  Inertial  navigation  sys¬ 
tem  (INS)  information  V/irometrlc  altitude  information,  and  range  and  range-rate  data 
from  ground  transponders  which  have  been  precisely  surveyed.  Although  the  navigation 
solution  produced  by  CIRIS  is  highly  accurate,  it  will  soon  be  inadequate  as  the  standard 
against  which  future  navigation  systems  car  be  tested.  This  research  proposes  an  alterna¬ 
tive  to  CIRIS  -  a  hybrid  Navigation  Reference  System  (NRS)  which  is  designed  to  take 
advantage  of  a  newer  INS  (the  LN-93),  certain  features  of  the  current  CIRIS,  and  certain 
features  of  the  Global  Positioning  System  (GPbj.  Analysis  is  conducted  using  a  Kalman 
filter  development  package  known  as  the  Multhnodet  Simulation  for  Optimal  Filter  Eval¬ 
uation  (MSOF'E).  Both  a  large  order  truth  model  for  the  NRS  (In  which  a  full  24  satellite 
constellation  is  modeled)  and  a  full-order  Kalman  filter  are  developed.  Results  suggest 
that  the  proposed  NRS  (with  GPS  aiding)  provides  a  significantly  Improved  navigation 
solution  as  compared  to  CIRIS. 
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A  NAVIGATION  REFERENCE  SYSTEM  (NRS)  USING 
GLOBAL  POSITIONING  SYSTEM  (GPS) 

AND  TRANSPONDER  AIDING 


/.  Introduction 

Thii  thesis  conititutei  an  extension  of  previous  work  in  which  th«  goal  has  been  to 
deilgn,  Implement,  and  analyze  an  improved  Kalman  filter  for  the  Completely  Integrated 
Reference  Initrumentation  Syitem  (CIRIS).  Developed  at  the  Central  Inertial  Guidance 
Test  Facility  (CIGTF),  Holloman  AFB,  NM,  CIRIS  function*  at  a  high  fidelity  navigation 
system  standard  or  reference  against  which  various  INS  mechanizations  are  tested. 

LI  Background 

In  order  to  quantify  the  accuracy  of  existing  or  proposed  navigation  systems,  It  is 
necessary  to  compare  the  navigation  solution  of  the  system  under  test  to  the  solution  pro¬ 
duced  by  a  reference  system  such  as  CIRIS  [1?,  18]  .  Additionally,  the  reference  system 
must  be  at  least  an  order  of  magnitude  more  accurate  when  compared  to  the  navigation 
systems  which  are  currently  tested  at  CIGTF,  as  well  as  those  which  will  be  tested  there 
in  the  near  future  [17,  28].  CIRIS  has  served  as  the  high  quality  navigation  system  lest 
standard  since  becoming  operational  in  1975.  However,  due  to  expected  improvements 
in  future  navigation  systems,  CIRIS  may  no  longer  be  adequate  as  a  navigation  system 
standard.  Consequently,  Incorporation  of  Global  Positioning  Syitem  (GPS)  aiding  to  Im¬ 
prove  the  CIRIS  navigation  solution  is  being  pursued.  It  Is  believed  that  GPS  aiding  will 
improve  the  navigation  solution  accuracy  well  beyond  the  accuracy  which  is  obtained  by 
the  current  implementation  of  CIRIS  [18,  32,  33]. 

Plans  ej;ist  at  CIGTF  to  Implement  an  upgraded  system  (similar  to  CIRIS)  which  is 
called  the  Advanced  Reference  System  (AIIS).  The  AHS  Is  designed  to  function  In  much 
the  same  way  that  CIRIS  currently  does,  but  is  implemented  with  a  state-of-the-art  inertial 
navigation  system  (INS)  and  Is  aided  by  GPS.  One  candidate  INS  being  considered  for  use 
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Ill  tln<  ATIS  |i  the  “strap-down”  Litton  LN-93  [17|.  The  accuracy  and  reliability  of  (lie  t  ing 
laier  gyro  (RLG)  Inertial  systems,  luch  as  the  LN-93,  make  them  better  candidates  for 
use  in  ARS  than  the  older  LN-39  glmbaled  INS  which  is  used  in  CIRIS.  The  newer  INS, 
along  with  the  addition  of  GPS  aiding,  will  ensure  that  ARS  is  a  more  reliable  test  stan¬ 
dard  against  which  future  navigation  systems  may  be  tested.  However,  the  performance 
specifications  of  ARS  remain  unchanged  with  the  exception  of  some  improvement  In  the 
vertical  channel. 

J.S  Problem  Statement 

In  order  to  accommodate  future  INS  testing  requirements,  a  reference  system  of  much 
higher  accuracy  is  required.  Such  accuracy  may  be  obtained  with  a  properly  Integrated 
system  including  INS,  OPS,  and  existing  transponders.  A  fundamental  first  step  In  this 
process  Is  the  development  of  a  truth  model  and  a  full-order  extended  Kalman  filter  which 
Incorporates  the  LN-93  INS  system  model,  transponder  aiding  (as  implemented  In  CIRIS), 
and  GPS  aiding,  must  be  designed,  implemented,  and  analyzed  [17,  I8j.  The  truth  model 
proposed  in  this  thesis  is  called  the  Navigation  Reference  System  (NRS). 

1.3  Summary  of  Applicable  Previous  Research 

Although  a  substantial  body  of  knowledge  exists  with  respect  to  the  GPS-INS  inte¬ 
gration  problem,  published  research  which  relates  specifically  to  GPS-CIRIS  or  GPS-AllS 
integration  is  limited.  Only  Solomon’s  research  specifically  relates  to  GPS-CIRIS  integra¬ 
tion  [32,  33].  Other  research  which  is  relevant  to  this  thesis  Includes  the  CIRIS  research 
performed  by  Snodgrass  [31]  and  a  group  project  completed  by  students  In  the  AFIT  nav¬ 
igation  sequence  In  which  the  LN-93  INS  error  model  Is  used  in  an  INS-OPS  integration 
problem  [11], 

1.3.1  GPS-INS  Integration.  Cunningham  addressed  the  problem  of  the  "filler- 
driving-fUter"  instability  which  can  occur  when  two  systems  (INS  and  GPS,  in  this  case) 
are  used  to  aid  one  another  while  each  of  them  employs  it  s  own  optimal  Kalman  filter  [7,  8]. 
Cunningham  assumed  a  local-level  (north,  east,  down)  INS  platform  using  three  gyros  and 
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three  accelerometers,  The  QPS  receiver  was  &  four-chnnnel  let  capable  r>f  simultaneously 
obtaining  range  meaiurementi  from  four  GPS  satellites  (alio  called  ipace  vehicles  or  SVs). 
GPS  range-rate  measurements  were  not  modeled  due  to  the  assumption  of  the  high-dynamic 
conditions  (which  would  preclude  an  actual  GPS  receiver  from  accurately  calculating  range- 
rate),  Cunningham  developed  truth  models  for  both  the  INS  and  GPS.  Kalman  filters 
based  on  the  truth  models  were  combined  to  produce  a  joint-solution  Kalman  filter  which 
wai  believed  to  be  the  belt  poulble  itructure  for  such  an  Integration  [7,  8]. 

Cunningham  concluded  that  a  combined  INS-OPS  mechanisation  provides  the  user 
with  the  high  accuracy  of  the  GPS  during  benign  conditions,  and  provides  the  beneficial 
characteristics  of  the  INS  during  high  dynamic  maneuvers  [7].  Kis  conclusion  is  consistent 
with  other  authors  who  discuss  the  INS-GPS  integration  problem  [6,  22]. 

1,3,2  GPS-CIR1S  Integration,  In  addition  to  being  directly  related  to  the  GPS- 
CIRIS  Integration  problem,  Solomon's  work  represents  a  significant  extension  of  Cunning¬ 
ham's  work  [7,  8]  because  the  QPS  error  model  is  improved  considerably  and  the  CIRIS 
extended  Kalman  filter  Is  redesigned  to  enhance  CIRIS  performance  [32,  33].  Solomon’s 
preliminary  GPS  error  model  [32]  (and  extensions  developed  in  research  conducted  as  part 
of  the  navigation  systems  class  project  referenced  previously  [11])  forms  the  basis  for  the 
GPS  error  model  used  In  this  thesis, 

The  thrust  of  Solomon’s  work  was  to  produce  an  Improved  truth  model  for  CIRIS, 
to  produce  a  reduced-  order  Kalman  filter  for  CIRIS,  and  to  Integrate  GPS  aiding  for  the 
CIRIS  navigation  solution.  In  his  thesis  [33],  ho  completed  a  127-state  CIRIS  truth  model, 
He  also  produced  a  reduced-order  Kalman  filter  which  had  70  states.  In  a  special  study 
[32],  Solomon  also  assembled  a  (stationary  SV)  GPS  error  model  and  implemented  It  In 
software. 

Solomon  concluded  that  the  70-state  Kalman  filter  emulated  the  127-state  truth 
model  very  well.  (Performance  indicators  for  the  70 -state  filter  appear  to  be  slightly 
better  than  a  127-state  Kalman  filter  which  was  bused  on  the  truth  model.  This  may  be 
attributable  to  lack  of  observability  of  certain  states  In  the  truth  model  [17,  25].) 
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J,H,3  CIRIS  Kalman  Filter  Improvements,  Solomon  mode  ll  clear  in  hi  a  con¬ 
clusions  that  he  believed  validation  of  the  models  he  had  developed  was  more  important 
than  simply  tuning  the  filters.  As  a  result,  SnodgraBs  had  as  his  thesis  objectives  to  im¬ 
prove  the  structure  and  efficiency  of  the  filter  software,  and  to  validate  a  refined  version 
of  the  70-state  Kalman  filter  originally  developed  by  Solomon.  Snodgrass  achieved  the 
first  oojective  after  considerable  software  development.  The  second  objective  was  satisfied 
by  comparing  the  refined  70-state  Kalman  filter  performance  to  truth  data  in  the  form 
of  position  and  velocity  data  collected  on  the  CIGTF  test  track.  Additional  comparisons 
were  made  between  the  improved  Kalman  filter  and  position  data  obtained  from  the  laser 
ranging  system  at  the  Yuma  test  range.  Due  to  the  amount  of  time  needed  to  accomplish 
the  primary  objectives  discussed  above,  Snodgrass  did  not  work  specifically  on  the  GPS 
integration  problem  [30] . 

Snodgrass  concluded  that  the  70-state  Kalman  filter  performed  better  than  the  older 
version  of  the  CIRIS  filter  in  most  cases.  The  exception  is  in  the  case  of  range-rate 
measurement  processing.  The  new  filter  range-rate  residuals  were  unacceptably  large  at 
certain  times  in  the  flight,  resulting  in  a  large  percentage  of  the  measurements  being 
rejected.  SnodgrasB  concludes  with  a  recommendation  that  GPS  integration  (within  a  full 
order  truth  model)  be  the  next  step  at  AFIT. 

1.4  Research  Objectives 

References  cited  in  the  previous  section  [32,  33,  31,  11]  form  the  foundation  upon 
which  much  of  this  research  is  based.  The  major  objectives  of  this  thesis  are  essentially  to 
consolidate  and  extend  previous  work  as  follows: 

1.  Revise  the  barometric  altimeter  error  model  in  improve  the  fidelity  of  the  vertical 
channel  aiding  model  contained  in  the  I.N-93  documentation  and  compare  the  per¬ 
formance  of  the  new  model  to  the  old  model. 

2.  Write  software  to  calculate  satellite  positions  of  a  24-GPS-SV  consteliution  for  a  t  wo 
hour  period. 

3.  Assemble  a  high  fidelity  truth  model  (NRS)  which  Integrates  GPS  aiding. 
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4.  Produce  a  fvll  -order  Kalman  Alter  bated  ou  the  truth  model. 

5.  Conduct  Monte  Carlo  analyses  to  determine  performance  characteristics  of  the  con- 
figurations  shown  below. 

(a)  Evaluate  performance  of  the  original  93~stB.te  INS  error  model  for  an  alignment 
and  for  a  flight  profile,  Initial  conditions  are  identical  to  those  used  In  the  Litton 
validation  work. 

(b)  Evaluate  performance  of  the  revised  (new  baro  model)  96-state  INS  error  model 
for  an  alignment  and  for  a  flight  profile.  Initial  conditions  are  Identical  to  those 
used  in  the  Litton  validation  work. 

(c)  Evaluate  performance  of  the  revised  06— ma-te  INS  error  model  for  an  alignment 
and  for  a  flight  profile.  Initial  conditions  are  those  for  Holloman  AFB,  NM. 

■  i 

(d)  Evaluate  performance  of  a  reduced  72-state  INS  error  model  for  an  alignment 
and  for  a  flight  profile.  Initial  conditions  are  those  for  Holloman  AFB,  NM. 

(e)  Evaluate  performance  of  a  98-state  reduced  INS  plus  ground  traneponder  error 
model  for  an  alignment  and  for  a  flight,  profile.  Initial  conditions  are  those  for 
Holloman  AFB,  NM. 

(f)  Evaluate  performance  of  the  128 -state  NRS  error  model  (which  Iniegrates  INS, 
RRS,  and  GPS  error  mode;*)  for  an  alignment  and  for  a  flight  profile.  Initial 
conditions  are  those  for  Holloman  AFB,  NM. 

1.5  Assumptions 

The  problem  described  in  the  problem  statement  requires  some  assumptions.  First, 
the  LN--93  error  model  is  accepted  as  the  truth  model  fur  the  INS  subsystem  or  (lie  NRS, 
with  tiie  notable  exception  of  the  baro-altit,ude  portion  of  that  model.  The  baro-altitude 
model  is  replaced  with  a  revised  model  as  discussed  in  Chapter  III.  Although  there  is 
some  verification  that  results  comparable  to  those  obtained  by  Litton  are  achieved  in  tills 
research,  exhaustive  INS  error  model  validation  is  not  undertaken. 
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The  error  model  for  the  ground  transponder  rauge/range-rate  system  (RRS)  from 
the  Snodgrass  thesis  [31]  is  used  as  the  error  model  for  the  RRS  subsystem  of  the  NRS. 
RRS  initial  conditions  and  tuning  parameters  are  also  extracted  from  [31]. 

A  good  error  model  for  the  GPS  error  sources  Is  extracted  from  separate  articles  by 
Cox  [6],  Martin  [21],  and  Milliken  and  Zoller  [27]  and  presented  in  [11].  The  error  Btates  in 
the  model  presented  in  [11]  are  used  as  the  basis  for  the  GPS  subsystem  in  the  NRS  model, 
although  different  initial  conditions  and  tuning  parameters  based  on  personal  experience 
are  used  in  some  cases  to  improve  performance  of  the  Kalman  filter. 

1.6  Research  Approach 

This  section  provides  a  general  outline  of  the  approach  used  to  accomplish  each  of 
the  tasks  listed  in  the  research  objectives  section.  The  correspondence  between  items  in 
this  section  and  those  In  Section  1.4  are  NOT  one-to-one;  some  tasks  overlap  with  one 
another. 

1.  Revision  of  the  barometric  altimeter  error  model  is  necessary  because  improving 
model  fidelity  will  enhance  the  system  (truth)  model  [17,  31].  The  revised  model  is 
based  on  the  Litton  LN-93  error  model  documentation  and  discussions  with  Lewan- 
towicz  [17]. 

2.  A  truth  model  and  full-order  benchmark  Kalman  filter  are  constructed.  The  truth 
model  is  assembled  as  follows: 

(a)  The  Litton  LN-93  documentation  is  the  basis  for  the  INS  error  model  [20].  The 
single  barometric  altimeter  state  in  the  LN-93  model  is  replaced  by  the  revised 
model  discussed  above.  Bnro-nltimeler  error  dynamic  equations  are  developed 
in  Chapter  III. 

(b)  The  RRS  transponder  portion  of  the  NRS  error  model  is  based  on  the  work 
performed  by  Snodgrass  [31]. 

(c)  A  revised  GPS  error  model  is  constructed  for  the  NRS.  GPS  5F  motion  is 
modeled  in  this  effort  whereas  stationary  5Rs  u>ere  assumed  previously  [11,  33]. 
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Software  which  explicitly  calculates  time- varying  SV  position  for  a  full  C! PS 
constellation  is  written. 

3.  Analyses  are  conducted  using  the  Multimode  Simulation  for  Optimal  Filter  Eval¬ 
uation  (MSOFE)  software  package  [5].  MSOFE  is  run  in  a  variety  of  modes  and 
configurations  in  order  to  establish  a  performance  baseline  for  the  (INS-GPS)  truth 
model.  This  objective  is  important  because  there  has  been  no  such  analysis  performed 
on  the  GPS  truth  model  to  date.  This  analysis  also  provides  a  baseline  against  which 
full-order  and  reduced-order  filters  may  be  comparod. 

1.7  Computer  and  Software  Usage 

Simulations  performed  as  part  of  this  research  require  a  substantial  amount  of  com¬ 
puting  power.  This  need  stems  from  the  calculation-intensive  nature  of  time  domain  sim¬ 
ulations  of  dynamic  systems  in  general  and  the  size  of  the  NRS  integration  problem  in 
particular.  Both  the  truth  model  and  the  full-order  Kalman  filter  each  have  128  states. 
The  filter  covariance  matrix  is  treated  as  an  upper  diagonal  array  hi  software,  yet  still 
contains  8,256  elements  (each  of  which  must  be  integrated  continuously).  A  “dedicated” 
VAX-station  III  (operating  at  approximately  3  MIPS)  is  inadequate  for  the  task.  Turn 
around  times  for  the  flight  simulations  are  on  the  order  of  several  days.  A  tremendous 
improvement  is  achieved  by  hosting  and  running  the  NRS  simulation  on  a  Micro Vax  III 
within  which  a  15  MIPS  coprocessing  board  1b  installed.  However,  the  problem  is  still 
nontrivial.  More  than  18  hours  are  needed  to  complete  a  10-run  Monte  Carlo  alignment 
simulation  with  128  states  in  the  truth  model  and  128  states  in  the  filter.  A  10-run  Monte 
Carlo  flight  simulation  (presented  in  Chapter  III)  is  expected  to  take  8  to  9  days. 

In  addition  to  the  software  which  is  written  explicitly  for  this  thesis,  three  software 
packages  of  a  more  general  nature  are  used  to  support  the  research  conducted.  MSOFE 
(the  Multi-mode  Simulation  for  Optimal  Filter  Evaluation),  written  by  Dr  Neal  Carlson 
and  Mr  Stan  Muslck,  functions  as  a  “shell”  within  which  proposed  or  existing  Kalman  filter 
designs  may  be  implemented  and  tested  (5).  PROFGEN  (short  for  PROFile  GENerator), 
also  written  by  Musick,  is  used  to  create  time  histories  of  variables  such  as  aircraft  position 
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ill  jf£  *  ipace,  attitude,  and  acceleration  [1].  The  time  histories  created  in  IIiIn  nmmi«r 
serve  as  the  “truth”  (or  nominal)  aircraft  trajectory  in  MSOFE  simulations  of  the  type 
conducted  in  this  thesis,  PROFGEN  is  used  to  create  the  aircraft  trajectory  in  the  two- 
hour  flight  simulation  which  is  discussed  in  more  detail  in  Chapter  Ill.  Data  which  are 
produced  during  the  MSOFE  “runs"  are  post-processed  and  plotted  using  MATRIX* 
[14).  Analyses  of  several  barometric  altitude  models  and  some  observability  analyses  are 
also  performed  using  MATRIX*, 

1.8  Summary  of  Planned  Research 

While  this  thesis  is  intended  to  consolidate  and  extend  the  work  begun  by  Cunning¬ 
ham,  Solomon,  and  Snodgrass  [7,  32,  33,  31],  it  departs  from  previous  efforts  somewhat, 
CIRIS  accuracy  as  a  reference  is  expected  to  be  inadequate  in  the  very  near  future,  If  it 
is  not  already.  CIGTF  is  in  the  process  of  designing  a  new  system  such  as  AllS  which 
provides  performance  comparable  to  CIRIS,  but  which  records  INS,  GPS,  and  transponder 
data  which  can  be  post-processed  with  high-accuracy  algorithms.  Therefore,  the  primary 
goal  of  this  thesis  is  to  design,  implement,  and  test  a  system  truth  model  and  a  benchmark 
Kalman  filter  which  is  called  the  the  Navigation  Reference  System  (NRS).  To  accomplish 
the  thesis  objectives,  a  system  truth  model  is  assembled,  implemented,  and  analyzed  using 
the  MSOFE  software  package. 

1.9  Thesis  Overview 

Chapter  II  provides  brief  descriptions  of  the  subsystems  which  compose  the  CIRIS, 
ARS,  and  NRS.  Reference  frames  and  general  Kalman  filter  theory  are  also  presented. 
Chapter  III  discusses  the  relevant  INS  theory.  The  LN-93  error  model  Is  presented  and  Its 
role  in  the  NRS  is  discussed,  The  revised  baro-altimeter  model  Is  developed  anti  discussed. 
A  full  development  of  the  RRS  error  and  measurement,  models  is  presented  in  Chapter  IV. 
A  similar  development  for  GPS  is  included  in  Chapter  V,  along  with  the  development  of 
GPS  SV  position  calculations. 

Chapter  VI  presents  the  results  of  simulations  and  analyses  for  both  the  truth  model 
and  the  benchmark  Kalman  filter.  Five  configurations  are  evaluated  and  compared.  Chap- 
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tor  Vlf  summarizes  the  theaii  effort.  Conclusions  are  drawn  from  the  results  presented  lu 
Chapter  VI  and  recommendations  are  made  for  future  research  efforts. 

Appendix  A  includes  the  complete  152  error-state  vector  in  a  set  of  tables  which  are 
arranged  to  emphasize  the  natural  divisions  in  the  model.  Appendix  B  is  obtained  from 
the  Litton  reference  [20]  and  shows  the  non-zero  elements  in  the  various  sub-blocks  in 
the  LN-93  dynamics  matrix  F,  and  the  non-zero  elements  of  the  noise  process  matrix  Q, 
which  is  also  extracted  from  the  Litton  reference  [20]. 

Appendix  C  includes  the  data  plots  from  the  “baseline”  case  of  an  INS  operating 
with  baro-aidlng  only.  Two  cases  are  shown:  first,  the  results  from  a  simulated  8-minute 
alignment  are  presented,  and  second,  the  results  from  a  two-hour  fighter  flight  simulation 
(which  is  explained  in  greater  detail  in  Chapter  III).  These  simulations  are  initialized  with 
parameters  which  were  used  by  Litton  in  establishing  the  baseline  performance  of  the 
LN-93  INS  [20],  The  initial  latitude  is  45  degrees  north,  and  longitude  is  zero  degrees. 
Initial  altitude  is  zero  feet.  Appendix  D  includes  plots  for  the  96-state  INS  (which  Includes 
a  new  baro-altimeter  model).  The  same  two  cases  are  shown  (i.e.  alignment  and  flight 
runs).  The  same  initial  latitude,  longitude,  and  altitude  as  discussed  above  are  used  in  this 
case.  Appendix  E  presents  plots  for  the  96-state  INS  (with  the  new  baro-altimeter  model). 
Both  the  alignment  and  flight  runs  are  repeated,  but  with  initial  latitude,  longitude,  and 
altitude  parameters  set  to  the  appropriate  values  for  Holloman  AFB,  NM. 

In  Appendix  F,  results  from  a  reduced  INS  error  model  are  presented.  The  INS  error 
model  is  reduced  to  72  states.  The  72-state  truth  model  error  behavior  is  compared  to 
that  of  the  98-state  truth  model. 

Appendix  G  Includes  plots  for  the  98-state  (reduced)  INS  and  RRS  configuration. 
The  same  two  cases  are  shown  (i,e.  alignment  and  flight  runs),  again  with  initial  conditions 
set  for  Holloman  AFB.  Appendices  H  and  I  each  contain  plots  for  the  128-state  INS/llRS 
and  GPS  configuration.  In  both  cases,  the  alignment  and  flight  run  results  are  shown,  and 
initial  conditions  are  those  for  Holloman  AFB.  Appendix  H  presents  resuits  from  operating 
NRS  with  GPS  aiding  only,  whereas  Appendix  i  demonstrates  NRS  performance  with 
measurements  from  both  the  RRS  and  GPS  subsystems. 
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II.  Subsystems,  Reference  Frames,  and  Kalman  Filtering 


2.1  Important  Subsystems  in  CIRIS,  ARS,  and  NRS 

Several  concepts  related  to  CIRIS  and  ARS  operation  are  central  to  understanding 
the  purpose  and  need  for  this  thesis  research.  Therefore,  It  is  useful  to  begin  with  a  brief 
description  of  the  subsystems  which  make  up  CIRIS  and  ARS.  NRS  subsystems  will  also 
be  described.  The  subsystems  which  are  defined  include  both  gimbaled  and  strapdown 
INS  mechanisations,  the  barometric  altimeter  and  central  air  data  computer  (CADC),  the 
radio  frequency  range/range-rate  transponder  system,  and  the  Global  Positioning  System 
(GPS).  The  details  provided  in  this  thesis  for  each  of  these  subsystems  are  only  at  that 
level  required  for  clear  understanding  of  the  subsystem  function  in  the  CIRIS,  ARS,  or 
NRS  “integrated”  systems.  It  is  not  within  the  scope  of  this  thesis  to  provide  extensive 
operational  or  functional  details  for  these  subsystems. 

2.1.1  Gimbaled  INS.  Britting  describes  inertial  navigation  as  the  process  of  de¬ 
termining  the  position  and  velocity  of  a  vehicle  (bucIi  as  an  aircraft,  ship,  or  space  vehicle) 
with  respect  to  a  specified  frame  of  reference  by  the  UBe  of  calibrated  electro-mechanical 
devices  installed  within  the  vehicle  (3).  A  typical  inertial  navigation  system  (INS)  makes 
use  of  the  output  of  gyroscopes  and  accelerometers  to  determine  position  and  velocity 
precisely,  Specifically,  Inertial  systems  perform  the  following  functions: 

1.  Instrument  a  known  reference  frame, 

2.  measure  specific  force  and  extract  accelerations  in  the  reference  frame, 

3.  perform  integrations  of  accelerations  to  obtain  velocity  and  position. 

A  gyroscope  (usually  referred  to  as  a  gyro)  is  a  device  which  traditionally  employs  a 
rapidly  spinning  mass  to  create  a  strong  angular  momentum  vector.  Since  changes  in  the 
angular  momentum  of  the  gyro  are  proportional  to  applied  torque,  the  gyro  can  maintain 
a  known  spatial  orientation  if  appropriate  torque  control  Is  applied.  Gyros  are  typically 
used  as  sensing  elements  in  closed-loop  servo  systems  which  operate  to  maintain  the  gyros’ 
spatial  orientation.  Thus,  spinning  mass  gyros  have  traditionally  been  the  devices  of  choice 
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to  Instrument  tin*  throe  dimensional  Cartesian  coordinate  frame  which  Is  referred  to  ns  the 
inertial  frame  in  aircraft  navigation  [3].  However,  some  navigation  systems  use  a  more 
recently  developed  inertial  instrument  known  as  the  ring  laser  gyroscope  (RLG)  [20]. 

RLG  construction  typically  consists  primarily  of  an  optical  cavity,  a  laser  device, 
three  mirrors,  a  prism,  and  a  pair  of  photodetectors  [29].  According  to  Savage,  the  RLG 
operates  as  follows.  Two  beams  of  (laser)  light  are  propagated  (opposing  one  another) 
around  the  path  enclosed  by  the  optical  cavity  [29].  If  the  cavity  is  rotating  in  an  inertial 
sense,  the  propagation  times  of  the  two  light  beams  is  different.  The  delay  manifests  itself 
in  the  form  of  a  phase  shift  between  the  two  beams,  and  the  phase  shift  is  detected  by  a  pair 
of  photodetectors  [29].  The  magnitude  of  the  phase  shift  provides  a  direct  Indication  of  the 
speed  of  angular  displacement  of  the  instrument  [29].  Devices  of  this  type  are  extremely 
reliable  due  to  the  absence  of  moving  parts  [29],  and  with  recent  developments  In  production 
techniques,  are  capable  of  accuracy  approaching  that  of  their  glmbaled  counterparts  [20]. 

Specific  force  is  measured  by  accelerometers.  The  most  common  accelerometers  to 
date  have  been  devices  which  are  sophisticated  variations  of  the  simple  pendulum  [3].  The 
motion  of  the  Internal  mass  is  related  to  the  inertially  referenced  motion  of  the  Instrument 
(case)  by  Newton’s  second  law  of  motion.  However,  to  obtain  the  correct  measure  of  inertial 
acceleration,  the  effects  of  local  gravity  must  be  removed  from  the  measured  specific  force 
[3,  18].  Consequently,  the  navigation  system  must  be  capable  of  distinguishing  between 
the  local  gravity  field  and  accelerations  which  are  being  applied  to  the  accelerometer  as 
a  result  of  flight  dynamics.  Otherwise,  significant  errors  will  appear  in  the  velocity  and 
position  calculations  and  these  errors  will  be  compounded  in  subsequent  system  updates 
[3,  29]. 

The  gyros  and  accelerometers  described  above  are  integrated  into  a  single  system 
called  the  INS.  An  example  of  a  widely  used  gimbalnl  INS  (as  described  above)  is  the 
Litton  LN-39.  The  LN-39  role  in  CIRIS  will  be  discussed  in  a  later  section. 

2.1.2  Strapdown  INS .  The  strapdown  INS  (also  referred  to  as  a  strapped-down 
INS  in  the  literature  [29])  is  conceptually  quite  similar  to  the  gimbaled  INS  referenced 
in  the  previous  section.  The  “strapdown”  system  obtains  its  name  from  the  fact  that 
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there  arc  no  platform  gimbals  In  the  system,  For  example,  the  LN  93  Is  mechanized 
by  Installing  three  single-degree-of-freedom  (SDOF)  ring  laser  gyros  (RLGs)  and  three 
SDOF  accelerometers  to  a  rigid  structure  (in  the  case  of  the  LN-93,  a  rigid  platform 
within  the  INS  “black  box").  An  internal  computer  maintains  a  current  computation  of 
aircraft  attitude  based  on  measurements  obtained  from  the  six  inertial  sensors.  Although 
strapdown  systems  were  previously  less  accurate  than  gimbaled  INS  mechanizations  [3], 
they  are  becoming  far  more  accurate  and  have  already  surpassed  the  precision  achieved 
by  some  of  their  gimbaled  counterparts  [18],  Because  it  uses  high  accuracy  RLGs,  the 
LN-93  is  a  prime  example  of  the  state-of-the-art  in  strapdown  INS  technology.  It  Is  for 
this  reason  that  the  LN-93  is  the  INS  selected  for  use  in  the  NRS. 

2,1.9  Baro- Altimeter  and  Central  Air  Data  Computer  (CADC).  A  well  docu¬ 
mented  shortcoming  of  any  INS  is  the  instability  which  (in  the  absence  of  aiding  infor¬ 
mation)  results  in  unbounded  error  growth  in  the  vertical  position  and  vertical  velocity 
channels  [3,  10,  16].  This  Inherent  instability  is  controlled  by  vertical  channel  aiding.  Such 
aiding  is  frequently  accomplished  with  vertical  position  information  provided  from  either  a 
barometric  altimeter  or  a  CADC.  Either  of  these  two  approaches  provides  vertical  channel 
stability  by  providing  additional  measurements  in  the  form  of  pressure  altitude  informa¬ 
tion.  This  external  altitude  Information  has  the  effect  of  stabilizing  the  vertical  channels 
which  would  otherwise  be  divergent  if  position  were  calculated  solely  on  the  basis  of  infor¬ 
mation  (measurements)  obtained  from  inertial  instruments  [3].  The  baro-altimeter  role  In 
the  NRS  will  be  discussed  in  more  detail  in  Chapter  Ill. 

2.1.4  Range/ Range- Rate  Transponder  System.  The  range/range- rate  transponder 
system  (RRS)  is  designed  specifically  to  augment  the  CIRIS  INS  in  order  to  Improve  the 
CIRIS  navigation  solution  accuracy.  The  RRS  system  consists  of  an  interrogator  unit 
which  is  carried  on  board  the  aircraft  which  carries  CIRIS,  and  approximately  forty  (40) 
fixed-site  ground  transponders  whose  positions  have  been  precisely  surveyed  [28].  The 
Interrogator  transmits  moderate  rate  digital  interrogation  codes  which  are  received  by  any 
of  the  ground  transponders  that  are  within  broadcast  range.  If  the  transponder  recognizes 
its  own  identifier,  it  will  respond  with  a  return  signal.  Range  measurements  are  obtained  by 
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comparing  the  phase  of  the  outgoing  (interrogator)  signal  with  the  Incoming  (transponder) 
signal.  Because  any  electromagnetic  signal  experiences  delays  when  propagated  through 
a  medium,  the  delay  which  is  measured  by  the  phase  shift  comparison  provides  a  direct 
Indication  of  slant  range  from  the  CIRIS  aircraft  to  the  t  ransponder  being  Interrogated  [28]. 
Note  that  the  RRS  transponder  system  is  alio  a  source  of  vertical  Information  which  aids 
in  improving  the  overall  navigation  solution  of  CIRIS,  ARS,  or  NRS.  The  RRS  transponder 
subsystem  is  discussed  in  detail  in  Chapter  IV. 

2.1,5  Global  Positioning  System  (GPS).  GPS  is  a  navigation  system  which  ii 
based  on  world-wide  coverage  of  a  constellation  of  24  space  vehicles  (called  SVs),  or  satel¬ 
lites  [12].  According  to  Milliken  and  Zollcr  [27],  the  GPS  provides  position  and  velocity 
information  to  users  located  anywhere  in  the  world,  with  anticipated  accuracy  on  the  order 
of  10  meters.  Normally,  the  user  needs  to  acquire  and  maintain  lock  on  four  SVs  in  order 
to  obtain  a  complete  navigation  solution.  However,  Milliken  and  Zoller  also  state  that 
fewer  SVs  are  needed  if  the  user  has  other  Information  (such  as  baro-altimeter  measure¬ 
ments)  available  for  supplementing  the  GPS  information.  The  GPS  navigation  solution  is 
obtained  in  a  manner  similar  to  the  RRS  transponder  system  described  previously. 

The  GPS  pseudo-range  between  the  user  and  each  SV  is  computed  based  on  knowl¬ 
edge  of  time  (the  master  GPS  clock)  and  the  unique  signal  format  which  is  broadcast  by 
each  SV.  Once  the  four  ranges  are  known,  a  recursive  algorithm  is  solved  to  compute  the 
uset’s  position  [27],  The  GPS  subsystem  model  is  presented  in  Chapter  V. 

2.2  Systems  Descriptions 

Subsystems  which  play  crucial  roles  in  the  CIRIS,  ARS,  and  NRS  systems  are  de¬ 
scribed  in  previous  sections.  The  complete  systems  are  described  below. 

2.2.1  CIRIS.  Since  1986,  the  0585th  Test  Group  at  Holloman  Air  Force  Base 
has  conducted  component  and  system  testing  to  analyze  and  verify  performance  of  state- 
of-the-art  inertial  navigation  systems.  Various  systems  have  been  used  as  the  standard 
against  which  systems  under  test  were  compared.  The  most  recent  of  the  systems  to  be 
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titllcd  the  toil  standard,  CIRIS  has  been  uied  with  great  eucceea  ilnce  107b.  Uidng  u 
Kalman  Alter,  CIRIS  optimally  combines  data  from  several  subsystems  to  obtain  a  highly 
accurate  navigation  solution  which  has  (until  recently)  met  or  exceeded  all  requirements 
for  navigation  reference  accuracy.  CIRIS  is  configured  to  be  carried  aboard  either  cargo 
or  fighter  aircraft.  Additionally,  CIRIS  may  be  operated  In  a  mobile  test  van.  This  mode 
is  used  lor  low  dynamic  tests  in  order  to  save  costs  associated  with  flight  testing. 

In  Its  current  configuration,  CIRIS  combines  INS,  CADC,  radio  frequency  RRS 
transponder  system  information,  and  a  time  referencing  system  which  Is  connected  to  the 
Holloman  Inter-Range  Instrumentation  Group  (IRIG)  time-keeping  system.  Currently,  two 
versions  of  CIRIS  are  available.  The  first  (CIRIS  I)  uses  the  Litton  LN-15  INS,  and  the 
second  (CIRIS  II)  uses  the  Litton  LN-39  INS,  Additionally,  the  CIRIS  II  system  uses  an 
updated  CADC  which  is  the  same  type  as  those  used  on  the  F-16  aircraft.  (CIRIS  II  is 
not  fully  operational  at  this  time  due  to  difficulty  with  data  collection  software  [28].) 

The  CIRIS  navigation  solution  has  been  shown  to  be  accurate  to  within  14  feot  in 
the  horizontal  direction,  and  40  feet  in  the  vertical  direction.  However,  CIRIS  accuracy 
will  soon  be  eclipsed  by  new  navigation  systems  under  development.  The  new  navigation 
systems  will  take  advantage  of  Improved  gyro  and  accelerometer  technologies  as  well  as 
GPS  aiding  [17,  18].  When  such  systems  are  produced,  they  will  rival  or  surpass  CIRIS 
accuracy.  Although  such  accuracy  in  navigation  systems  is  in  fact  desirable ,  there  is  no 
reliable  means  of  testing  the  new  systems  unless  a  superior  test  standard  can  be  devised. 
Thus,  the  need  for  the  ARS  (and/or  NRS)  is  established. 

2.2.2  ARS.  The  Advanced  Reference  System  (ARS)  is  under  development  at 
Holloman  AFB  and  is  intended  to  provide  a  more  reliable  reference  system  than  CIRIS, 
even  though  comparable  performance  capabilities  are  expected  [17].  When  ARS  becomes 
operational,  performance  of  future  navigation  systems  maybe  validated  by  comparing  them 
to  the  (post-processed)  solution  from  ARS,  in  much  (he  same  way  that  current  systems 
are  compared  to  CIRIS  I  or  II. 

ARS  will  also  employ  a  Kalman  filter  which  optimally  combines  dBta  from  several 
subsystems  including  an  INS  (type  undetermined),  a  CADC  (type  undetermined),  and 
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the  WIS  which  it  uiud  by  CIRIS  I  and  II.  ARS  will  also  make  use  or  OPR  pseudo  range 
information  which  ii  now  available.  Became  of  the  higher  fidelity  component!  and  the 
integration  of  GPS  aiding,  ARS  accuracy  ahould  improve  significantly  over  that  of  CIRIS 
[28].  ARS  it  being  packaged  in  an  AIM-9  missile  pod  in  order  to  facilitate  external  carriage 
of  the  reference  system.  This  form  of  packaging  allows  quick  installation  and  removal  of  the 
ARS  reference  package  on  aircraft  which  are  used  for  flight  testing  of  navigation  systems. 

2,2.8  NRS,  The  Navigation  Reference  System  (NRS)  includes  a  Kalman  filter 
which  combines  subsystems  Including  an  INS  (with  a  baro-altimeter  model),  the  RR5 
used  by  CIRIS,  and  GPS.  The  NRS  uses  the  LN--93  INS  and  a  revised  baro-altimeter 
model  as  well.  As  a  result,  NRS  should  be  a  viable  test  standard  for  navigation  system 
testing,  and  in  particular,  may  possibly  be  used  as  a  “post-processing"  Kalman  filter  and 
smoother  [17]  for  the  ARS.  (Such  use  Is  contingent  upon  which  subsystems  are  chosen  for 
ARS.  If  incompatible  subsystems  are  used,  the  NRS  model  is  likely  to  be  unsuitable  for 
ARS  applications.) 

2.3  Reference  fYamc a 

A  navigation  “solution’’  has  significance  only  If  the  coordinate  frame  in  which  the 
solution  is  expressed  is  clearly  understood.  While  the  preceding  statement  may  seem 
obvious,  it  cannot  be  overemphasised,  Consider  that  the  Litton  LN--93  documentation 
defines  the  earth  frame,  the  true  frame,  the  computer  frame,  the  platform  frame,  the 
sensor  frame,  the  gyro  frame,  the  accelerometer  frame,  and  the  body  frame  (20). 

Hence,  there  is  an  urgent  need  for  precise  notation  which  aids  In  the  communication 
process.  There  is  also  a  need  for  the  ability  to  transform  quantities  which  are  expressed  in 
one  frame  Into  appropriate  expressions  in  another  frame.  This  need  gives  rise  to  transfor¬ 
mation  matrices. 

Solutions  to  both  of  the  needs  discussed  above  nre  provided  by  Brlttlng  [3].  The 
notation  is  Introduced  as  needed.  Prior  to  defining  the  transformation  matrices,  several 
reference  frame  definitions  are  needed.  Figure  (2,3)  depicts  the  central  reference  frames 
which  Brlttlng  has  chosen  to  define.  They  are  discussed  hi  the  next  section. 


2-6 


Figure  2.1.  Coordinate  Frame  Geometry  [3] 

2,4  Reference  Frame  Definitions 

The  coordinate  frame  definitions  below  are  those  which  are  significant  to  this  thesis 
[3].  Figure  (2.3)  depicts  the  first  three  frames  which  are  defined.  Figure  (2.3)  uses  the 
following  notation; 

A  =  Longitude  with  respect  to  inertial  frame 

A/  =  Longitude  with  respect  to  ECEF  frame 

la  =  Initial  longitudinal  displacement  of  ECEF  from  Greenwich 

meridian  (usually  defined  as  ZERO) 

L  =  Geographic  frame  latitude 

h,-  =  Geocentric  frame  latitude 

w,,,  =  Earth  angular  rate  (rarl /sec) 


1.  Inertial  frame:  an  orthogonal  R:l  coordinate  system;  its  origin  is  coincident  with  the 
earth's  center  of  mass  and  the  frame  is  oriented  as  follows.  The  y,  plane  lies  in 
the  earth’s  equatorial  plane  and  does  not  rotate  with  respect  to  the  fined  stars.  The 
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$i  axis  alwayi  polnti  toward  Ariel.  The  s,  axis  projects  from  the  earth's  center  nf 
man  directly  through  the  North  pole.  (This  frame  li  depicted  by  the  [0,  ]/,  x]  frame 
in  Figure  (2.3.) 

2.  Earth-Centered-Earth-Flxed  (ECEF)  frame:  an  orthogonal  R:>  coordinate  system; 
its  origin  is  coincident  with  the  earth's  center  of  mass.  The  x,-  ,  y,  plane  lieu  in 
the  earth's  equatorial  plane  and  the  xr_  axis  is  usually  chosen  to  be  aligned  with  the 
Greenwich  meridian.  The  s„  axis  projects  from  the  earth’s  center  of  mass  directly 
through  the  North  pole.  (This  frame  rotates  at  exactly  the  earth  rate,  u>ln  ,  about 
the  i„  axis.) 

3.  Geographic  (or  navigation)  frame:  an  orthogonal  coordinate  system;  its  origin  is 
at  the  location  of  the  INS  (or  the  user),  and  its  axes  are  aligned  with  the  North, 
East,  and  Down  directions  [N,E,D].  This  description  implies  that  the  INS  platform  is 
torqued  to  maintain  the  [N,E,D]  orientation.  (It  must  be  noted  that  not  all  platforms 
are  torqued  in  this  manner.  Some  are  allowed  to  “wander”  about  the  s-axls,  and 
are  torqued  in  the  level  axes  only.  Platforms  that  are  not  torqued  about  the  vertical 
axis  are  called  wander  azimuth  systems  [3].) 

4.  Body  frame:  an  orthogonal  frame;  its  origin  is  at  vehicle  (i.e.,  aircraft)  center  of 
mass.  Its  axes  are  the  vehicle’s  roll,  pitch,  and  ,yaw  axes  [R,P,Y].  Britting  points 
out  that  the  origin  of  the  body  frame  rarely  (If  ever)  coincides  with  the  origin  of  the 
navigation  frame  [3]. 

5.  Platform  frame:  an  orthogonal  frame  that  “can  be  thought  of  as  three  fiducial  lines 
which  are  physically  inscribed  on  the  platform.”  [3] 

6.  Gyro  frame:  a  (possibly)  non-orthogonal  frame;  its  axes  are  defined  by  the  input  axes 
of  the  three  gyroscopes  Installed  In  the  platform.  The  degree  nf  non-nrthngoiinllty 
depends  upon  manufacturing  and  physical  Installation  considerations  [3,  16]. 

7.  Accelerometer  frame:  a  (possibly)  non-orthogonal  frame;  its  axes  are  defined  by  tho 
three  input  axes  of  the  accelerometers  installed  in  the  platform.  Again,  the  degree  of 
non-orthogonality  depends  on  manufacturing  and  Installation  considerations  [3,  16J. 
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Willi  Drilling'*  definition*  In  hand,  the  Litton  reference  frame  definition*  are  stum*- 
what  eaiier  to  decipher.  The  Litton  reference  frame*  a*  specified  in  [20]  are  interpreted 
with  the  help  of  the  Britting  definition*  as  follows: 


1.  Earth-Centered-Earth-Fixed  (ECEF)  frame:  an  orthogonal  R:I  coordinate  system; 
its  origin  is  coincident  with  the  earth’s  center  of  mas*.  The  z,.  ,  plane  lies  in  the 
earth'*  equatorial  plane  and  the  xf  axis  li  chosen  to  be  aligned  with  the  Greenwich 
meridian.  The  ye  axis  project*  from  the  earth's  center  of  mass  directly  through  the 
North  pole.  (This  frame  rotates  at  exactly  the  earth  rate,  ,  about  the  p,  axis.) 
NOTE:  This  frame  is  not  the  ECEF  frame  described  by  Britting.  The  transformation 
between  the  two  frames  is: 
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2.  True  frame:  a  level  R :l  coordinate  system  located  at  the  user's  actual  latitude  and 
longitude.  If  a(  =  0°  ,  then  the  x,  axis  points  to  true  East,  the  y,  axis  points  to 
true  North,  and  the  z <  axis  points  exactly  Up,  resulting  in  the  Litton  true  frame 
[U,  N,  U],.  (This  is  the  error-free  case  of  Brittlng’s  navigation  reference  frame.) 

3.  Computer  frame:  a  level  R 1  coordinate  system  located  at  the  user’s  indicated  latitude 
and  longitude.  If  ct,..  -  0°  ,  then  the  axis  points  East,  the  y,:  axis  points  North,  and 
the  zK  axis  points  Up,  resulting  in  the  Litton  computer  frame  [22,  JV,  [/]„.  This  frame 
is  defined  by  coordinates  which  are  calculated  by  the  computer  that  is  integral  to 
the  INS.  Because  of  the  presence  of  uncertainties  in  the  computer  frame  calculations, 
discrepancies  exist  between  quantities  (i.e.  latitude  and  longitude)  calculated  in  tills 
frame  and  the  same  quantities  in  the  true  frame. 

4.  Platform  frame:  this  frame  is  important  in  the  system  error  dynamic  equations 
where  the  platform  [mlsjalignment  with  respect  to  the  true  frame  must  be  taken  into 
account.  The  result  of  phytically  installing  an  INS  (situating  the  "black  box”  in 
an  aircraft,  for  example)  is  often  a  slight  misalignment  of  the  platform  frame  with 
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respect  to  the  true  frame.  The  reiultlng  error  angles  are  defined  by  n  skew-symmetric 
transformation  <f>i 
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5.  Body  frame:  Litton  defines  the  body  frame  in  identically  the  same  way  as  does 
Brlttlng.  For  the  LN-93,  the  sensor  frame  is  defined  to  be  identically  the  same  as 
the  body  frame. 

0.  Gyro  frame:  Litton  also  uses  the  Brlttlng  definition  for  the  gyro  frame.  This  frame 
is  important  in  the  system  error  dynamic  equations  where  gyro  misalignments  with 
respect  to  the  navigation  frame  (true  and/or  computer  frame)  must  be  taken  into 
account.  The  physical  construction  of  the  INS  “black  box,”  including  strapdown 
mechanizations  using  RLGs,  is  never  perfect;  the  result  is  that  the  gyro  frame  is  non- 
orthogonal  to  a  certain  extent  [3, 18]  although  this  is  negligible  for  some  applications 
[18]. 


Note  that  Litton  assumes  the  body,  sensor,  and  platform  frames  are  coincident  [20]. 
The  same  assumption  is  maintained  for  this  study  because  INS  platform  location  with 
respect  to  the  body  frame  is  specific  to  each  aircraft  type.  When  the  platform  and  body 
frame  origins  are  NOT  coincident  (as  in  most  applications)  the  “lever  arm"  effect,  discussed 
in  [31]  must  be  taken  into  account.  However,  in  the  interest  of  generality  and  efficiency, 
the  lever  arm  effect  is  not  considered  in  this  study.  The  gyro  and  accelerometer  frame 
errors  ace  considered  and  are  appropriately  reflected  in  the  LN-93  error  model. 


2.JJ  Reference  Frame  'rransformation  Matrices.  The  LN-93  INS  is  a  local-level, 
wander  azimuth  platform.  Consequently,  user  position  1b  specified  in  terms  of  latitude, 
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longitude,  altitude,  and  wander  angle,  Generally,  wander  angle  is  the  angle  about,  the 
n-f :ame  t  axil  which  reiulti  because  the  INS  platform  is  not  torqued  to  maintain  precise 
[ENU]  orientation.  In  this  mechanization,  the  platform  is  "torqued”  only  in  the  “level”  axes 
in  order  to  maintain  the  local  level  orientation.  (The  LN-93  h  a  Btrapdown  mechanization 
which  employs  RLGs,  Consequently,  "torquing”  takes  place  only  in  software.) 

It  is  often  necessary  to  express  vectors  such  as  position,  attitude,  velocity,  or  accel¬ 
eration  in  terms  of  several  different  reference  frames.  As  an  example,  the  Litton  LN-93 
INS  error-state  model  describes  position  error  in  terms  of  an  error-angle  space  vector, 
[W.r,  S0V,  60, ,  fh]1 ,  where  6$x  is  the  error  angle  about  the  local  level  E  axis,  S9U  is  the 
error  angle  about  the  local  level  N  axis,  60-  is  the  error  about  the  local  level  U  axis,  and 
fh  is  the  altitude  error.  Even  though  Litton’s  definition  is  clear,  If  the  error-angle  vector  is 
to  have  physical  meaning,  it  must  be  transformed  into  a  vector  in  navigation  error  space, 
[6<j>,  f\,  fa,  fh  l7 .  where  6</>  is  the  error  In  latitude,  6 A  is  the  longitude  error,  fa  is  the 
wander  azimuth  error,  and  fh  is  again  the  altitude  error.  A  transformation  matrix,  C", 
permits  compact  transformation  of  the  error-angle  vector  into  an  equivalent  expression  in 
navigation  error  space. 

Brltting  developed  general  transformation  matrices  for  many  reference  frame  pairs. 
Those  transformation  matrices  which  are  of  primary  interest  in  this  research  are  shown 
below.  The  first  four  matrices  are  attributed  to  Brltting  [3]  and  the  last  matrix  is  developed 
by  Litton  [20]: 


c; 
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true- to- earth 

C'n 
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navigation-to-earth 
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true-to-navigatlon 
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error  anglc-to-navigation 

Usage  of  these  transformation  matrices  is  illustrated  in  the  following  example  in  which  a 
vector  which  is  written  in  the  j -frame  is  transformed  to  an  equivalent  expression  in  the 
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cos  A  cos  a,  —  sin /..sin  Ailna,  —  (n?»  Asht  a,  4-  sin /.sin  A  cos  a*)  sin  A  cos/. 

coil  lino,  cos  L  cos  a,  sin  L 

(*in  A  cos  o.  A-  cos  A  »in  L  nin  a, )  sin  A  sin  o,  -  cos  A  sin  L  cos  o,  cos  A  co*  L 

(2.5) 

cos  A  --  sin  A  sin  L  sin  A  cos  L 
0  cos  L  siuL  (2*0) 

-  sin  A  -  cos  A  sin  L  cos  A  cos  L 

cos  a,  -  sin  a,  0 

CJ1  =  sin  a,  cos  a,  0  (2.7) 

0  0  1 

cos  0  sin  1 p  cos  0  cos  i p  sin  9 

C/  =  sin  <t>  sin  9  sin  ip  +  cos  <p  cos  ip  sin  <p  sin  0  cos  ip  -  cos  <j>  sin  ip  -sin^cos  9  (2.8) 

cob  <p  sin  0  sin  0  -  sin  <p  cos  ip  cos  <p  sin  9  cos  ip  -I-  sin  ipcoatp  -  cos  f  cos  9 

In  the  equation  for  C\'  above,  <p  is  the  aircraft  roll  angle,  0  Is  the  aircraft  pitch  angle,  and 
ip  is  the  aircraft  heading  angle. 

To  conclude  the  discussion  on  reference  frame  transformations,  the  matrix  C'J  is  pre¬ 
sented.  C'J  converts  a  vector  written  in  the  Litton  error  -angle  space,  [W,;,  69u,  69.,  bh  ]f 
into  one  in  navigation  error  space  [<50,  6 A,  6a,  6k}1,  The  transformation  matrix  is  shown 
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below: 


C 


tl 

II 


-  cob  n 
sin  a  sec  4> 
-  sin  a  tan  <j> 
0 


sin  a  0  0 

cos  a  sec  <f>  0  0 

~  cos  a  tan  ^10 
0  0  1 


This  matrix  is  used  in  the  transformation 


’  6<t>  ' 

'  M,  ' 

a 

r 

69y 

•  St, 

os 

5a 

. 

56, 

6h  _ 

6k 

(2.9) 


(2.10) 


2,5  WGS-81  Geodetic  System 

In  light  of  the  numerous  reference  frames  which  are  discussed  in  the  previous  section, 
it  is  apparent  that  the  designer  may  choose  any  one  of  a  number  of  coordinate  systems 
in  which  to  define  position.  The  current  most  popular  geodetic  reference  frame  (ECEF) 
is  defined  by  the  World  Geodetic  Survey  of  1084  (WGS-84).  According  to  Kumar  [15], 
the  WGS-84  reference  Is  the  state-of-the-art  in  geodesy,  reflecting  the  latest  and  best 
modeling  techniques  for  characterizing  the  earth’s  surface. 

The  WGS-84  model  is  used  exclusively  as  the  ECEF  reference  for  this  thesis  effort. 
Several  of  the  key  parameters  from  the  standard  are  Included  in  Table  2.1  [9,  31].  The  JVn 
V,.,  and  Z,  coordinates  are  simply  those  which  define  the  orthogonal  ECEF  axes  of  the 
WGS-84  model.  The  earth’s  angular  rate  is  represented  by  w,e,  and  the  equatorial  and 
polar  radius  of  the  spheroid  are  given  by  A  and  B,  respectively.  The  flattening  parameter 
is  defined  as  [3]: 

f  =  ~r  (2.n) 

where  A  and  B  are  defined  above,  and  eccentricity  is  [3]: 


e  = 


(2,12) 
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Table  2.1.  WGS  84  Parameter* 


Parameter 

Definition 

Value 

ECEF  Coordinate  Frame  Axes 

not  applicable 

w.v 

Angular  Rate  of  the  Earth 

7.292115.  IF*  s~* 

A 

Semimajor  Axis 
(Equatorial  Radius) 

0378137  m 

B 

Semiminor  Axis 
(Polar  Radius) 

6356752.3142  m 

e 

First  Eccentricity 

0.0818191908426 

Flattening  (Elllptlcity) 

0.00335281066474 

Equatorial  Acceleration 
of  Gravity 

2.6  Kalman  Filter  Theory 

Most  estimation  problems  which  are  not  purely  academic  are  fraught  with  uncer¬ 
tainty.  The  problem  of  navigation  is  a  superb  example.  A  navigation  “solution"  is  ob¬ 
tained  by  solving  some  rather  complex  equation*  which  depend  on  several  quantities  that 
are  known  to  varying  degrees  of  accuracy.  As  a  result,  the  accuracy  of  the  solution  Is 
limited.  It  is  the  rote  of  the  Kalman  filter  to  account  for  the  uncertainties  associated  with 
such  a  problem,  and  to  provide  an  optimal  solution  (under  certain  conditions).  Although  it 
is  far  beyond  the  scope  of  this  thesis  to  provide  extensive  details  on  the  topic  of  Kalman  fil¬ 
tering,  it  is  essential  to  present  a  small  amount  of  the  basic  theory  along  with  the  equations 
which  describe  the  Kalman  filter  algorithm. 

Maybeck  describes  a  Kalman  filter  as  an  “optimal  recursive  data  processing  algo¬ 
rithm"  which  depends  upon  three  basic  assumptions  [22|.  If  a  linear  (or  linearized)  system 
is  driven  by  white,  Gaussian  noise,  the  Kalman  filter  solution  is  “optimal  with  respect  to 
virtually  any  criterion  [22].” 

By  optimally  combining  data  (measurements),  dynamic  characteristics  and  statistical 
properties  of  the  system  and  measuring  devices,  and  initial  conditions,  the  Kalman  filter 
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produces  estimates  of  the  modeled  quantities  which  are  more  accurate  than  estimate* 
based  on  the  individual  measurements.  The  development  which  follows  Is  a  summary  of 
important  Kalman  Altering  concepts  as  presented  by  Maybeck  [22,  23]. 

S,  6. 1  Linear  Kalman  Filtering,  Before  Kalman  Altering  is  undertaken,  the  designer 
must  develop  a  mathematical  description  for  the  system  of  interest.  When  it  is  possible  to 
do  so,  a  system  it  modeled  at  a  set  of  linear  differential  equations  of  the  form  [22]t 


x(t)  -  P (t)x(t)  +  B(!)u(t)  +  G(t)w(<)  (2.13) 


where 


x 

F 

B 

u 

G 


“state”  vector 

homogeneous  state  dynamics  matrix 
control  input  matrix 
deterministic  control  input  vector 
driving  noise  input  matrix 
white  Gaussian  driving  noise  vector 


Because  the  deterministic  control  term  Bu  is  zero  in  this  research,  it  will  be  ignored 
hereafter.  Also  note  that  for  this  thesis,  G  =  I.  The  mean  of  the  white  Gaussian  driving 
noise  vector  is: 


mw  -  £{w(t)}  =  0 


(2.14) 


and  the  noise  strength  is  Q(t): 


E  {w(t)w7'(t  +  r)}  =  Q(t)*(r)  (2.1B) 


While  Equation  (2.13)  is  written  in  terms  of  ‘whole’  value  state  variables,  the  models 
used  in  this  thesis  are  those  of  error  states.  This  choke  of  state  variables  results  in  simpler 
dynamic  equations  [3],  and  Equation  (2.13)  may  be  rewritten  as  [22]: 


tf*(t)  =  F  (t)Sx{t)  +  B(t)u(<)  +  G(t)w(t)  (2.16) 
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where  x(t)  has  boon  replaced  by  the  error  state  vector  dx(t),  and  all  other  quantities  retain 
their  previous  definitions.  The  topic  of  error  states  ii  explored  more  fully  in  the  section  on 
extended  Kalman  filters. 

As  previously  stated,  the  Kalman  filter  incorporates  measurement  Information  from 
external  measuring  devices.  Irrespective  of  the  type  of  measuring  device,  the  equation 
which,  is  used  to  describe  linear  measurements  is  of  the  form; 

«(!,)  =  H(t,)x(ti)  +  v(t,)  (2.17) 

or,  in  the  case  of  error-state  models: 

«*(*,)  =  H(t,)tfx(t,)  +  v(t,)  (2.18) 

where,  in  both  case#  above,  H  Is  the  observation  matrix  and  v  is  a  discrete-time  zero-mean 
white  Gaussian  measurement  noise  vector  whose  covariance  is  [22]: 

/ 

£{v(t, )vr(<,)}  = 

The  Kalman  filter  “propagates”  the  error  state  and  its  covariance  from  the  instant 
In  time  Immediately  following  the  most  recent  measurement  update,  t*,  to  the  Instant  in 
time  immediately  preceding  the  next  measurement  update,  t“+ ,  by  numerical  integration 
of  the  following  equations  [22]: 


R  far  t,  ~  tj 
0  for  ti  £  tj 


(2.19) 


&(t/t,)  =  P(f)5c(t/t.)  (2.20) 

P {t/ti)  =  F(t)P(t/<i)  +  P(t/t, )*''  ('/',)  -I-  G(t)Q(t)Gr(t)  (2.21) 

with  initial  conditions: 

S(t/ti)  =  x(tD  (2.22) 

P(*/*i)  -  P(C)  (2.23) 
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ai  provided  by  the  meaeurement  update  cycle  at  time  t, , 

After  propagation,  £(<)  and  P(t)  arc  "updated”  (meaning  that  itate  eitimatei  are 
revUed  baied  on  new  Information).  The  pi  vote',  element  in  the  update  equation*  ihown 
below  i*  the  time  varying  Kalman  filter  gain  K(t).  The  K(t)  matrix  auigni  "weight*” 
to  the  “new  information”  (new  Information  consists  of  the  difference  between  the  actual 
meaeurement  and  the  filter’*  estimate  of  the  measurement)  based  on  known  measurement 
noise  statistics  and  filter-computed  covariances  from  the  previous  time  step.  This  process 
Is  designed  to  improve  the  estimate  of  each  element  of  the  state  vector.  The  update 
equations  are  [22]; 


K(t,)  =  P(t,r)Hr(t,)[H(t,)P(t-)H7(t,)  +  R(t)]- 

(2.24) 

*i*t)  ■---  *(*,")  +  K(t,)[ss,  -  H(t)8(<-)] 

(2.2B) 

P(tt)  =  P(tr)-K(t,)H(t,)P(tr) 

(2.28) 

Although  the  algorithm  shown  above  is  generally  applicable  to  any  problem  which 
lends  itself  to  a  Kalman  filtering  solution,  It  is  not  necessarily  the  algorithm  which  Is  used 
in  practice.  It  is  often  advantageous  to  use  a  form  of  the  algorithm  known  as  the  U-D 
form.  In  the  U-D  algorithr  ,  the  filter  covariance  matrix  is  not  propagated  as  a  square 
array.  The  matrices  below  representing  the  pre-  ami  post-measurement  filter  covariances, 
respectively,  are  explicitly  computed  instead  [22]: 

P(<r)  =  U(tr)D(t,-)U(t-)  (2.27) 

P(^)  =  U(t+)D(t/)U(t+)  (2.28) 

where  the  U  matrices  are  upper  triangular  and  unitary  (and  thus  contain  ones  along 
the  main  diagonal),  and  the  D  matrices  are  simply  diagonal  [22],  This  form  offers  several 
advantages  including  numerical  stability,  improved  precision,  and  guaranteed  nonnegativity 
of  the  computed  covariances  [22j.  It  is  the  U-D  form  of  the  Kalman  filter  algorithm  which 
is  Implemented  in  the  MSOFE  software  [5]. 
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2.6,2  Linearized  and  Extended  Kalman  Filtering.  Unfortunately,  not  nil  problem* 
are  adequately  described  as  linear  systems  driven  by  white  Qausslan  noise.  In  many  cases, 
the  most  appropriate  model  is  nonlinear.  The  navigation  problem  at  hand  falls  squarely 
into  the  nonlinear  category.  Fortunately,  a  method  exists  whereby  a  nonlinear  system  may 
be  treated  in  much  the  same  manner  as  a  linear  one  for  a  particular  class  of  problems. 
Suppose  that  the  nonlinear  system  may  be  described  by  [23]: 


*(t)  «  f  (x(t),  u(t),  t]  +  G(t)w(t)  (2.29) 

In  this  case,  the  state  dynamics  matrix,  f[>,  v],  is  a  nonlinear  function  of  the  state 
vector  x(*),  time  t,  and  the  control  input  (assumed  to  be  eero  in  this  research).  The 
white  Gaussian  noise  is  defined  exactly  as  in  Equations  (2.14)  and  (2.15).  In  addition,  the 
measurement  equation  may  also  be  a  nonlinear  function  of  time  [23]: 

»(*.)  =  h[x(l,),t,]  +  v(t.)  (2.30) 

The  noise  vector  v  is  again  zero-mean  and  its  covariance  is  described  by  Equation  (2.19). 

Recalling  that  a  system  must  be  linear  in  order  to  satisfy  the  assumptions  that 
yield  the  optimality  of  a  Kalman  filter,  the  nonlinear  Equations  (2.29)  and  (2.30)  must  be 
linearised.  The  following  approach  is  summarized  from  Maybeck  [23]: 

1.  Assume  that  a  nominal  state  trajectory,  x„(t),  may  be  generated  which  satisfies 
Xr.(<u)  =  and 

x„(t)  =  f[x„(t),u(t),t]  (2.31) 

where  ([*,•,■]  is  specified  In  Equation  (2.20).  nml  u(t)  =  0. 

2.  The  “nominal"  measurements  which  accompany  the  nominal  trajectory  are: 

M«0  *  h[x, ,(/,), f(]  (2.32) 
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3.  The  "perturbation"  of  the  itatc  is  obtained  by  subtracting  the  nominal  trajectory 
from  the  original  nonlinear  equation: 

WO  -  *»(*))  =  ~  flx„(t),u(t),t]  +  0(t)w(<)  (2.33) 


4.  The  equation  above  may  be  approximated  to  first  order  by  a  Taylor  series  expansion; 


a  F[(|Xn(0]  +  G(i)w(t) 


(2.34) 


where  6x(t)  represents  a  first-order  approximation  of  the  process  [x(>)  -  x, ,(•)),  and 
F[t;x„(t)]  is  a  matrix  of  partial  derivatives  of  /  with  respect  to  lt.s  first  argument, 
evaluated  along  the  nominal  trajectory  [23]: 


F[<ix„(t)]  = 


iXaXit(f) 


(2.35) 


6.  The  perturbation  measurement  equation  is  derived  In  like  fashion  and  Is  expressed 
as  [23]: 

(5*„(<)  =  H[t,;x, ,(<,)]  +  v(f,)  (2.30) 


where 


H[f {*„(<)]  = 


8h\x{t),t)\ 


ix=xN(/) 


(2.37) 


With  the  “error-state"  model  in  hand,  It  is  possible  to  return  to  the  linear  filtering  theory. 
An  estimate  of  the  whole-valued  quantities  of  interest  is  obtained  from  [23]: 


S(f)  --  x„(t)  +  6x(t) 


(2.38) 


The  expression  above  for  the  linearized  Kalman  filter  is  useful  provided  that  the  lineariza¬ 
tion  assumption  is  not  violated.  However,  if  the  nominal  and  "true"  trajectories  differ 
by  too  large  an  amount,  unacceptable  errors  mav  result  [23].  It  is  for  this  reason  that 
extended  Kalman  filtering  Is  useful  in  many  cases  where  perturbation  techniques  alone 
do  not  suffice.  Extended  Kalman  filtering  allows  for  relincariztng  about  newly  declared 
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nominal*  at  each  sample  time,  to  enhance  the  adequacy  of  the  lhienrtiotion  prom*,  mnl 
thui  of  the  reiultlng  Alter  performance  ai  well  [25]. 

The  extended  Kalman  Alter  equation*  are  summarised  below.  The  reader  la  referred 
to  Maybeck  [23]  for  the  detail*  of  derivation.  The  measurement  equation  for  an  extended 
Kalman  Alter  is: 

«(<,)  =>  h[x(t,),i,]  +  v«.)  (2.39) 

where  v(>)  1*  once  again  tew-mean  with  covariance  given  by  Equation  (2.19).  Measure¬ 
ments  arc  incorporated  Into  the  extended  Kalman  Alter  via  the  following  set  of  equations 
(23]: 


K(tf)  -  P(<r)H7'[(t1);S(<r )]  {H[(t,');X(tr)]P(ir)Hr[(t,)|«(tr)J  +  &(<)}“'  (2.40) 

*(^)a*(<r)  +  K(t,){*.  -h[K(t(-)i  (*,)(}  (2.41) 

p(<<+)  -  m)  -  K(t,)H[(t,)i«(t-)]p(t-) 


whore 


XkX(»(" 


(2.42) 


(2,13) 


The  state  estimate  and  covariance  are  propagated  from  t,  to  t,+  i  by  Integrating  the  fol¬ 
lowing  equations  [23]: 

k(t/t,)  =  f[jt(t/t,),u(t),t]  (2.44) 

£(</<,)  =  F[t;S(t/t,)]P(t/tl)  +  P(f/f,)F7  (t;X(f/t,)]  4-  Q(t)q{t)Gr (t)  (2.43) 


where 


0f[x(t),  u(t).  f] 

b* 


(2,40) 


Ffo*(*/f,)]  = 

and  the  Initial  condition*  are: 

*(*/«.)  =  *(C)  (2.47) 

P(t/t,)  »  P(f,'  )  (2.48) 

The  equations  shown  above  for  the  extended  Kalman  Alter  are  programmed  Into  the 
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MS  OPT,  shell  for  the  problem  defined  by  tide  thesis.  It  is  the  feet  that  the  extended  Kalman 
lllter  li  relinearized  about  each  successive  estimate  of  the  state  [£(<)]  wldch  “enhances  the 
'validity  of  the  aeeumptlon  that  deviations  from  the  reference  (nominal)  trajectory  are  small 
enough  to  allow  linear  perturbation  techniques  to  be  employed"  [23]. 

8,7  Summary 

This  chapter  introduces  several  topics  which  are  related  to  the  NRS  integration  prob¬ 
lem.  Major  systems  such  as  the  INS,  the  RAS  'Transponder  system,  and  the  QPS  are 
described  briefly  and  their  roles  as  subsystems  in  the  NAS  are  outlined,  Reference  frames 
are  discussed  In  moderate  detail  and  several  important  frames  are  defined.  Coordinate 
transformations  are  also  supplied.  The  WGS-84  Geodetic  System  Is  presented  and  key 
parameters  are  tabulated.  Finally,  a  very  brief  outline  of  some  of  the  Kalman  filter¬ 
ing  theory  concepts  crucial  to  this  thesis  are  presented.  Equations  for  the  linearized  and 
extended  Kalman  Alters  are  shown  and  tiie  validity  of  using  the  extended  Kalman  Alter 
algorithm  for  the  nonlinear  problem  in  this  thesis  is  explained. 
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Ill,  INS  and  Baro~Altimeter  Moduli 


3.1  Introduction 

The  INS  chosen  for  use  in  the  proposed  Navigation  Reference  System  (NIIS)  Is  the 
Litton  LN-93.  The  LN-93  Is  a  state-of-the-art  strapriown  INS  which  uses  RLGh  in  its 
mechanisation  rather  than  the  older  gimbaled  gyros.  The  LN-93  characteristics  promise 
a  considerable  advantage  over  mechanisations  used  in  previous  navigation  standards  such 
as  CIRIS  (he.,  the  LN-1G  used  in  CIRIS  I  and  the  LN-39  used  in  CIRJS  II).  The  LN-93 
error  model  is  presented  below,  Additionally,  the  need  for  a  revision  to  the  baro-altimeter 
model  is  presented  and  the  revised  model  is  developed. 

3.2  Litton  LN-93  INS 

Litton  specifications  stipulate  that  the  LN-93  horizontal  velocity  errors  must  be 
smaller  than  2.5  feet  per  second  (RMS)  after  two  hours  of  navigation  following  a  gyrocom¬ 
pass  alignment,  and  that  horlsontal  position  error  must  be  less  than  0.8  nautical  miles  per 
hour  during  the  same  period  [20],  Actual  statistics  provided  in  the  Litton  documentation 
show  that  for  a  period  of  ten  hours  of  (static)  navigation,  the  LN-93  is  accurate  to  within 
2  nautical  miles.  For  a  2  hour  “flight”  simulating  a  fixed  flight  profile  (which  is  discussed 
in  Section  3.4),  the  LN-93  maintains  horlsontal  position  accuracy  to  better  than  1  nautical 
mite. 

Figure  3.1  depicts  the  basic  baro-alded  INS  subsystem.  The  INS  subsystem,  aided 
by  a  baro-altimeter,  forms  the  foundation  upon  which  the  full  NRS  model  is  built.  NRS 
will  eventually  incorporate  RRS  and  GPS  models,  as  presented  in  Chapters  IV  and  V, 
respectively. 

3,2,1  Error  2 types  Litton’*  comprehensive  dynamic  error  model  for  the  LN-93  is 
composed  of  93  states.  The  error  sources  modeled  fall  into  one  of  several  error  source 
types.  First,  the  random  constant  represents  a  variable  which  may  take  on  any  value 
within  prescribed  limits  (defined  by  a  random  value  chosen  from  a  Gaussian  distribution 
of  appropriate  scale)  and  thereafter  retains  its  initial  value  [22].  Variables  of  this  type  are 
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Figure  3.1.  Basic  Baro* Aided  LN-B3  INS 
often  referred  to  au  biases  and  are  modeled  mathematically  by  [22]: 

i(0  =  0  ,  and  P{t)  =  P»  (3.1) 

Approximately  64  bias  error  state  variables  are  included  In  the  Litton  model.  These  are 
lumped  Into  subcategories,  Including  gyro  bias  states,  accelerometer  bias  states,  trend 
states,  and  gyro  compliance  states. 

An  error  variable  type  not  represented  in  the  LN-93  model  is  the  random  walk  or 
Brownian  motion  variable.  Although  the  Litton  modct  does  not  include  such  variables, 
Maybeck  recommends  that  variables  of  this  type  be  used  In  the  Kalman  filter  Implementa¬ 
tion  of  random  bias  states  when  there  exists  a  tendency  for  the  covariances  of  those  states 
to  drop  to  zero  [22].  This  practice  (adding  small  magnitude  “pseudo-noises”  to  states 
whose  truth  model  includes  no  noise  source)  prevents  covariances  of  those  states  from 
dropping  to  zero,  thereby  Improving  the  estimation  potential  for  those  states.  Variables 
of  tills  type  are  characterized  by  a  tendency  to  wander  about  In  an  unpredictable  fashion, 
and  are  dezcribcd  a«  “the  output  of  an  integrator  driven  by  white  Gaussian  noise”  [22]. 
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They  are  represented  mathematically  by  [22]: 


z(t)  a=  w(t)  ,  and  P(t )  =  P„  (3.2) 

where  w  represents  zero-mean  white  Gaussian  noise  of  strength  q.  Variables  of  this  type 
are  not  used  in  this  error  model.  However,  they  are  generally  useful  when  Kalman  filters 
containing  bias  states  are  implemented  [22]. 

Another  error  source  type  included  in  the  LN-93  model  is  the  first-order  Markov 
process.  This  variable  is  also  characterized  by  a  tendency  to  wander  about,  but  with  a 
certain  degree  of  correlation  from  one  time  to  another.  This  process  is  described  as  the 
output  of  a  first  order  lag  which  is  driven  by  zero-mean  white  Gaussian  noise  [22],  The 
first-order  Markov  variable  is  modeled  by  an  integrator  with  a  negative  feedback  loop.  The 
feedback  gain  is  the  reciprocal  of  the  time  constant  for  the  variable  of  interest  [22].  Its 
mathematical  representation  is  [22]: 

i(t)  =  -  i  *(<)  +  u>(<)  ,  and  P(t)  ~  Po  (3.3) 

where  w  represents  zero-mean  white  Gaussian  noise  of  strength  q ,  and  r  is  the  time 
constant  associated  with  the  state  of  interest.  Error  variables  which  are  included  in  this 
category  are  gyro  and  accelerometer  drifts  and  thermal  transients  associated  with  these 
inertial  instruments.  (Additionally,  the  baro-altlmeter  state  Included  by  Litton  falls  into 
this  category.  This  is  discussed  in  more  detail  later.) 

The  last  category  of  errors  in  the  Litton  model  are  those  which  are  termed  “general” 
[20].  These  error  states  are  complex  combinations  of  several  states.  Some  states  which 
fall  Into  this  category  depend  on  as  many  as  40  to  50  other  states.  Navigation  position, 
velocity,  and  platform  tilt  error  states  are  examples  of  such  complex  interdependency,  The 
“lower”  13  states  (termed  the  general  error  states)  are  coupled  to  one  another  and  to  oilier 
states.  The  Litton  linearized  error  model  is  presented  In  the  form  of  a  general  equation  of 
the  type  presented  in  Chapter  II.  The  generalized  INS  error  state  equation  is  of  the  form 
[3,  22]: 

Sx(t)  =  F(<)fx(i)  +  Gw(<)  (3.4) 
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where 


(3.5) 


r(t)  -  gx  lx«x, 

and  Sx(t)  is  the  93  x  1  time  varying  ste*  error  vector,  F(f)  is  a  93  X  93  time  varying  system 
dynamics  matrix,  u(t)is  a  zero  vector,  w(i)  is  a  93  x  1  white  noise  vector,  and  G  -  I. 

3.2,2  Error  State  Vector,  Litton  partitions  the  error  state  vector  into  six  sub¬ 
vectors  as  shown  below  [20].  Note  that  each  of  the  subvectors  is  expiessed  as  one  of 
the  error  types  discussed  above. 

Sx  —  [  SxJ  6xJ  Sx SX'f1  6xrJ  Sxar  (3.6) 

where  dx  is  a  93  X  1  column  vector  and: 

Sx,  represents  the  “general”  error  vector  containing  position,  velocity,  attitude,  and  ver¬ 
tical  channel  errors. 

Sx-j,  consists  of  gyro,  accelerometer,  and  baro-altlmeter  correlated  errors,  and  “trend” 
states.  These  states  are  modeled  as  first  order  Markov  processes  In  both  the  truth 
(system)  model  and  in  the  Kalman  filter. 

Sx:t  represents  gyro  bias  errors.  These  states  are  modeled  as  random  constants  in  the  truth 
model  and  are  modeled  as  random  walks  (with  small  magnitude  pseudo-noises)  in 
the  Kalman  filter. 

ix,(  is  the  accelerometer  bias  error  states.  These  states  are  modeled  in  exactly  the  same 
manner  as  the  gyro  bias  states. 

Sxr>  depicts  accelerometer  and  gyro  initial  thermal  transients.  The  thermal  transient  states 
are  first  order  Markov  processes  in  the  system  mid  Kalman  filler, 

Sxu  models  the  gyro  compliance  errors.  These  error  states  are  modeled  as  biases  in  the 
system  model  and  as  random  walks  in  the  Kalman  filter. 
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The  system  state  space  deferential  equation  is  given  as: 


’ «.  ' 

'  F„ 

Fis 

F.a 

F|.i 

Fir, 

F  io 

d*. 

’  ' 

w, 

0 

Fas 

0 

0 

0 

0 

wa 

0 

0 

0 

0 

0 

0 

**:» 

0 

»  n 

•  4-  - 

0 

0 

0 

0 

0 

0 

S*4 

0 

**  8 

0 

0 

0 

G 

F5r, 

0 

$x5 

0 

0 

0 

0 

0 

0 

0 

Sx  o 

0 

Ail  states  In  the  93  state  INS  error  vector  are  defined  in  Appendix  A,  Tables  A.l 
through  A. 4.  All  non-zero  elements  of  the  dynamics  matrix,  F,  and  the  process  noise 
matrix,  Q,  are  obtained  from  [20]  and  are  presented  in  Appendix  B. 

3.3  Baiv- Altimeter  Model  Revision 

Although  the  Litton  error  model  is  extensive  and  reasonably  thorough  in  its  depletion 
of  Inertial  instrument  error  sources,  it  is  somewhat  inadequate  in  its  depletion  of  the 
error  sources  related  to  the  baro-altimeter.  Other  authors  discuss  several  error  sources 
normally  linked  with  baro-altitude  [16,  28,  34]  which  appear  to  have  been  omitted  in 
the  Litton  LN-93  model.  The  LN-93  error  model  contains  only  a  single  state  for  the 
barometric  altimeter.  For  this  reason,  a  revised  baro-altitudr  error  model  is  developed 
and  embedded  in  the  INS  error  model.  The  new  baro-altitude  model  includes  states  for 
correlated  noise  effects,  bias  error,  and  scale-factor  error.  The  revised  baro  model  coefficient 
for  the  correlated  noise  state  is  extracted  from  the  Litton  documentation  [20].  In  the  case 
of  the  other  (new)  baro  states,  the  coefficients  are  extracted  from  [34].  These  error  sources 
are  combined  in  the  manner  discussed  below  to  formulate  a  more  complete  error  model 
upon  which  to  base  subsequent  work. 

3.3.1  Litton  Vertical  Channel  Aiding.  As  noted  above,  the  Litton  LN-93  error 
model  contains  only  one  state  with  which  to  model  error  sources  Intrinsic  to  the  baro- 
altimeter.  The  state  is  a  first  order  Markov  process,  and  is  modeled  mathematically  as 
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[20): 


i(i)  =  a(t)  +  w(t) 


(3.8) 


where  w  represents  zero-mean  white  Gaussian  noiBe  of  strength  x  /?*/,„  is 

the  inverse  correlation  time  for  the  baro-altimeter,  and  fflu  is  the  variance  of  the  baro- 
altimeter  correlated  noise  as  specified  by  Litton  [20].  Litton  specifies  a  correlation  time  of 
600  seconds  and  <Tfhf  =  100. Qft  for  the  two  hour  fighter  run.  The  initial  covariance  of  the 
barn-altimeter  state  is  assumed  to  be  P0  =  lO  '/t3.  The  state  is  represented  in  Figure  3.2. 


Figure  3.2.  LN-93  Baro- Altimeter  Error  Model 

The  single  state  baro  model  aids  the  vertical  channel  states  (i.e.  INS  altitude  and 
vertical  velocity)  in  order  to  prevent  the  vertical  channel  instabilities  discussed  in  Chapter 
II.  The  baro  error  state  is  coupled  directly  to  four  states  in  the  dynamics  equations,  and 
is  indirectly  coupled  to  more  than  a  dozen  others.  Figure  3.3  depicts  the  vertical  channel 
error  model  included  in  the  Litton  reference  [20).  The  baro  error  Shu  enters  the  INS 
vertical  channel  error  model  on  the  right  side  near  the  middle  of  the  diagram.  It  is  notable 
that  Shu  la  approximated  by  6hr,  the  correlated  barn  alt  itude  error,  in  the  Litton  dynamics 
model.  Thus,  even  though  the  Litton  document  discusses  several  error  sources,  (anti  in 
fact  defines  total  baro-altimeter  error  to  be  a  function  of  bias  error,  scale  factor  error,  and 
correlated  noise  error)  these  error  sources  are  not.  modeled  (without  explanation)  and  a 
single  state  6hr,  is  used  to  represent  the  baro-altimeter  model. 

Because  the  altimeter  error  magnitude  is  significant  to  many  of  the  INS  states,  tin 
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Figure  3,3.  LN-93  Vertical  Channel  Error  Model  [20] 

fidelity  of  the  baro  error  model  Is  of  considerable  concern.  It  is  precisely  to  Improve  the 
its  fidelity  that  the  baro  error  model  is  revised  as  described  in  the  next  section. 

3. 3, ,2  Revised  Baro-Altitude  Model,  The  single  baro  state  included  in  the  Litton 
error  model  Is  adequate  to  represent  one  of  the  major  sources  of  baro  altitude  error.  Its 
relevant  equations  and  parameters  are  exactly  as  shown  in  Equation  (3.8).  However,  other 
baro-altimeter  error  sources  are  present  and  should  he  modeled  in  a  high-fidelity  truth 
model. 

Another  significant  source  of  error  in  the  baro-altimeter  which  must  be  taken  into 
account  is  the  error  due  to  instrument  bias,  This  error  source  is  well  modeled  as  a  random 
constant  whose  Initial  one  <r  value  is  chosen  to  be  15  feet  in  this  study,  The  choice  is  based 
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on  a  baro-altimeter  mOtbucr  trained  in  [34]. 

^  Another  slgnifican.  .  ./^  <>ttrce  included  in  the  revised  model  is  the  scale  factor  error 

inherent  in  baro  altimeters.  This  error  is  a  multiplicative  combination  of  aircraft  altitude 
|  and  the  random  constant  scale  factor  state  as  shown  below.  The  initial  one  o  value  for  the 

*  scale  factor  state  is  0.01  [34].  Although  the  error  contribution  for  this  state  is  negligible 

I 

at  low  altitudes,  it  can  be  significant  source  of  error  at  higher  altitudes, 

Tlw  last  error  state  In  the  revised  baro-altimeter  model  is  the  first  order  lag  which 
sums  and  filters  the  previously  described  errors.  The  output  of  this  state  represents  the 
1  total,  lagged  baro-altimeter  error  6hf).  It  is  shown  along  with  the  other  revised  baro  model 

t  states  in  Figure  3.4. 

I 


Figure  3.4.  Revised  Baro- Altimeter  Error  Model 

The  revised  baro-altimeter  model  shown  in  Figure  3.4  is  incorporated  into  the  NRS 
model.  It  occupies  states  23,  24,  25,  and  26  in  the  NRS  state  vector  as  shown  in  Appendix 
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A.  State  3!)  In  the  total  baro-altlmeter  error,  Shu.  State  24  Is  the  correlated  bnrn  altitude 
error,  state  25  it  the  biai  error  itate,  and  state  26  li  the  scale  factor  error  itato.  The 
performance  and  validation  of  the  revised  baro  model  are  discussed  in  Chapter  VI. 

3.4  Simulated  Fighter  Flight  Profile 

Litton  has  analysed  the  error  model  for  the  LN-93  in  several  modes  [20].  One 
important  example  is  performance  analysis  of  the  LN-93  when  it  is  “stressed1*  in  a  manner 
typical  of  a  fighter  mission  flight  profile.  In  other  words,  when  the  LN-93  is  subjected  to 
rapid  acceleration,  velocity,  and  attitude  changes  typical  of  a  fighter  mission  flight  profile, 
what  are  the  performance  characteristics  of  the  INS?  In  an  attempt  to  answer  this  question, 
the  flight  profile  shown  in  Figure  3.4  is  simulated  In  software, 

Trajectory  data  for  the  flight  profile  shown  In  Figure  3.4  are  calculated  using  PROF- 
GEN  [1].  The  trajectory  data  Include  variables  such  as  latitude,  longitude,  altitude,  ac¬ 
celerations,  velocities,  and  attitude.  These  data  are  used  in  the  MSOFE  [6]  simulation 
as  “nominal1’  quantities  about  which  the  truth  (system)  model  is  relinearised  after  each 
integration  step  [22].  The  extended  Kalman  filter  is  relluearlsed  about  its  best  estimate  of 
the  trajectory  quantities. 

3.5  Summary 

This  chapter  introduces  the  main  error  model  types  which  are  generally  important  to 
Kalman  filtering  and  specifically  important  to  the  Litton  LN-93  error  model.  The  LN-93 
error  model,  which  consists  of  93  states  as  provided  by  Litton,  is  introduced.  Because  it 
has  only  one  baro-altimeter  error  state,  the  LN-93  error  model  is  revised  to  improve  its 
fidelity;  as  a  result,  a  four-state  baro-altlmeter  model  is  used  in  the  overall  INS  model. 
The  flight  profile  which  is  used  for  subsequent  Kalman  filter  evaluation  Is  presented  and 
briefly  discussed. 
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Figure  3.5.  Fighter  Mission  Flight  Profile  [20} 


IV,  Range/Range-Raie  Transponder  System 


1,1  /nfrodvciion 

The  Range/Range-Rate  System  (RRS)  it  a  navigation  aiding  system  which  comprises 
a  significant  part  of  CIRIS  and  ARS.  Consequently,  the  RRS  is  used  in  the  Navigation 
Reference  System  (NRS)  as  well.  Navigation  Information  is  obtained  by  "interrogating" 
ground  transponders  and  subsequently  processing  the  electromagnetic  (EM)  signals  which 
the  transponders  emit.  The  information  obtained  allows  high  qusdity  range  and  range-rate 
measurements  to  be  calculated  by  the  RRS  interrogating  hardware  [31].  Using  these  range 
and  range-rate  measurements,  refinements  to  the  NRS  position  and  velocity  estimates  are 
then  possible. 

The  dynamic  error  model  for  the  RRS  transponder  system  is  contained  in  [31,  33]. 
Portions  of  that  work  are  summarised  in  this  chapter.  Note  that  the  lever-arm  effect  dis¬ 
cussed  in  [31]  is  not  Included  in  this  work.  Although  the  lever-arm  effect  is  very  important 
to  consider  in  actual  hardware  applications,  the  assumption  of  INS  and  RRS  interrogntor 
antenna  collocation  does  not  reduce  the  validity  of  the  conclusions  which  may  be  drawn 
from  this  work,  The  assumption  of  collocation  is  made  in  order  to  maintain  generality 
but  has  the  advantage  of  reducing  computational  loading.  For  applications  in  which  the 
collocation  assumption  is  not  valid  [such  as  the  case  in  which  actual  data  are  used],  a 
transformation  must  be  applied  to  translate  and  rotate  the  RRS  measurement  into  the 
navigation  frame.  One  such  transformation  matrix  is  developed  by  Snodgrass  [31].  How¬ 
ever,  the  transformation  matrix  developed  by  Snodgrass  [31]  assumes  a  [N,W,U]  navigation 
coordinate  system  (as  used  in  the  LN-30  INS)  rather  than  the  [E,N,U]  system  which  Is 
assumed  In  the  Litton  error  model  for  the  LN  03. 

Figure  (4.1)  depicts  the  addition  of  the  RRS  transponder  subsystem  to  the  bnro- 
aided  INS  diagram  presented  In  Chapter  III.  This  is  another  step  in  the  path  toward  the 
full  NRS  model  (which  eventually  incorporates  GPS  as  well). 
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Figure  4.1.  Baro-Alded  INS  with  RRS  Transponder  Aiding 
4.2  RRS  Range  Mcaaurementi 

In  CIRIS  (or  the  proposed  NRS),  RRS  range  measurements  aid  in  estimating  posi¬ 
tion  errors  of  the  reference  INS.  The  RRS  range  measurement  is  derived  from  the  time 
delay  detected  between  the  time  at  which  the  reference  hardware  (CIRIS  or  NRS)  sends 
an  interrogation  signal  and  the  time  at  which  a  reply  from  the  transponder  Is  received. 
This  temporal  difference  is  multiplied  by  the  speed  of  light  [and  divided  by  two,  to  account 
for  the  "round  trip”  of  the  electro-magnetic  (EM)  signal]  to  obtain  an  uncorrected  range 
measurement.  Correction  factor*  are  then  applied  in  order  to  compensate  for  delays  Intro¬ 
duced  by  the  propagation  of  EM  signals  through  the  atmosphere  and  to  correct  for  errors 
Introduced  by  equipment  calibration  biases  [31]. 

Figure  (4.2)  shows  an  earth-user- transponder  system  In  two  dimensions.  A  single 
transponder  is  shown  for  clarity  in  this  discussion.  However,  six  to  ten  transponders  are 
typically  used  during  reference  navigation  system  Kalman  filter  updating. 

The  true  positions  of  the  transponder  and  the  user  are  shown  along  with  the  true 
range  Ri  and  the  (uncorrected)  range  measurement  RMWS.  (which  Is  the  range  from  the 
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SXT 


Figure  4.2.  Two  Dimemlonal  Earth-User-Transponder  Syitem 

ground  transponder  to  the  uier,  or  more  specifically,  to  the  user’*  INS  platform).  It 
should  be  noted  that  both  the  true  range  and  the  uncorrected  range  are  collnear  with 
the  true  uier  and  transponder  positions.  The  vectors  SXT  and  5XU  are  the  transponder 
and  user  position  error  vectors,  respectively.  The  moment-arm  effect  between  the  INS  and 
HRS  receiving  antenna  is  not  represented  in  the  Agure;  It  is  assumed  for  this  discussion 
that  the  RRS  interrogator  antenna  Is  collocated  with  the  INS.  The  (uncorrected)  range 
measurement  as  obtained  from  RRS  is; 

HjtHy  —  R<  -f  dHalm  4"  dRfcr  +  V  (d.l) 

where 


TL„,)S  w  RRS  range  measurement,  from  transponder  to  user 
B,i  =  True  range,  from  transponder  to  user 

SR, ,ir,  =  Range  error  due  to  atmospheric  delay 

SRi„  =  Range  error  due  to  equipment  calibration 

v  =  sero-mean  white  Gaussian  measurement  noise 
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Equation  (4.1)  U  a  model  for  the  range  as  determined  by  the  RRS  truth  model.  It 
includei  the  true  range  [which  in  never  precisely  known]  along  with  termi  which  reflect 
lourcei  of  error  and  uncertainty.  RRS  error  sources  and  models  are  discussed  later. 


4.3  Range  Calculation  From  INS  Data 

In  order  to  formulate  a  difference  measurement  as  discussed  in  Chapter  II,  two  sources 
of  range  Information  must  be  obtained.  The  first  Is  the  RRS  range  measurement  which  is 
modeled  by  Equation  (4,1).  Another  range  indication  Is  computed  from  the  INS  indicated 
position  and  RRS  (surveyed)  positions.  (Both  the  indicated  INS  position  and  transponder 
surveyed  position  contain  uncertainties  which  must  be  considered.)  In  this  approach,  the 
user  (INS)  Indicated  position  Is  represented  by  an  RA  vector  expressed  In  the  Litton  ECEF 
as:  (> 

{*"  1 

y„  (-1.2) 

*U  J 

while  the  true  RRS  transponder  poiltioh  is  represented  in  the  Litton  ECEF  by: 


*, 


(4.3) 


Then  the  calculated  range  from  the  user  (INS)  to  the  transponder  is  given  by: 


^ins  ~  |X„  —  X, 


Equation  (4.4)  may  be  rewritten  as: 
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Tf  tb®  equation  above  li  “perturbed"  to  reflect.  uncertainties  in  uiier  and  transponder  po¬ 
sitions,  then  a  first-order  Taylor  lerlei  may  be  written  to  approximate  the  range  [3,  22], 
The  truncated  (flrit-order)  lerles  it  of  the  form  [33]: 
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After  lubitltuting  Equation  (4.5)  into  Equation  (4.6)  and  evaluating  the  partial  derivative*, 
the  INS-derived  range  approximation  become*  [31,  33]: 
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At  thli  point,  the  difference  meaiurement  may  be  formed  as: 

Sr.  -  Ef/vs  -  Rut,* 


•  4'*,,  - 

'Vr  -■  Vn 

'6ytl  - 

'*7'  -  ■*(-] 

\rZ,,\ 

\Ri«s\ 

si 
1  ! 

sd 

•  SxT  4- 

Vr  ~Vv' 

’  Sy  ,  4- 

Z  y  -  Z„ 

L  U*,n,I 

l*nvJ 

1^/fV.vl  • 

'S*u 

6zr 


-  [1]  SRulm  -  [l]tffl!()f  +  e 


(4.8) 


Note  that  the  true  whole- valued  range  ( Ri )  formerly  preient  in  both  individual  range 
meaiurement*  ii  cancelled  in  the  differencing  operation  [3,  22].  Alio  note  that  the  bracketed 
coefficient!  in  the  equation  above  will  be  prominent  in  the  H  matrix  development  diicuiied 
in  Section  4.7. 

In  order  to  form  the  difference  measurement,  it  in  assumed  above  that  the  INS  and 
transponder  coordinates  are  exprened  in  the  same  frame.  However,  this  azaumption  and 
reality  are  at  oddi.  The  RES  transponder  whole- valued  quantities  are  generally  expressed 
in  the  liGEF  frame  (as  defined  by  Brittlng).  A  simple  transformation  (presented  in  Chapter 
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T1)  ii  used  to  express  the  transponder  positions  in  the  Litton  ECF.F  frame.  Litton  [HO] 
expresses  the  whole- valued  TNS  position  vector  as  a  triplet  [latitude,  longitude,  altitude] 
which  can  be  converted  to  the  (Litton)  ECEF  frame  by  [11]: 
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terrestrial  longitude 
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a« /■/i’'~Kr«nrr5 

earth’s  equatorial  radius 
eccentricity  of  earth’s  ellipsoid  1 


Although  both  position  vectors  are  now  (apparently)  expressed  In  the  same  frame,  Equa¬ 
tion  (4.9)  depends  explicitly  on  the  fact  that  the  user  position  is  known  in  terms  of  latitude, 
longitude,  and  altitude.  The  discussion  in  Chapter  III  presents  the  LN-93  error-angle  vec¬ 
tor  [  59,,,  Sffu,  50,  ]  which  may  be  transformed  into  latitude,  longitude,  and  altitude.  The 
approach  Is  presented  in  detail  in  Section  4.6. 

4-4  RRS  Error-State  Model  Equations 

The  RRS  error  state  vector  is  composed  of  26  elements  (shown  in  Table  A. 6,  Ap¬ 
pendix  A).  The  RRS  states  occupy  numbers  S,xs  +  1  through  5,lVi<.  +  26  in  the  NRS 
error  state  model,  where  SINS  represents  the  total  number  of  states  used  to  model  the  INS 
subsystem  in  NRS.  (The  variation  of  Slrih  is  dim  uased  in  Chapter  6.)  The  llrst  two  RRS 
states  are  simple  random  constant  (bias)  states  which  model  the  effects  of  user  hardware 
[RRS  interrogator]  range  and  range-rate  calibration  errors,  respectively,  Coupling  of  these 
states  into  the  RRS  measurement  equation  is  presented  in  Sections  4.5  and  4.7.  The  error 
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itl.(iie  model  equation  for  these  states  is  [31]: 


where 


i&br 

*ljl> 


0  0 
0  0 


Hr 


®fc  II 


(4.10) 


wi,r  ~  range  equivalent  of  interrogator  bias 
~  velocity  equivalent  of  interrogator  bias 


The  initial  state  estimates  and  covariances  for  these  states  are  [31]: 


and 


MM  1 

0 

'  MM  J 

0 

.  i 

(4.11) 


P|ir,(m(M 


l/t*  0 
0  10  -'fit/sec* 


(4.12) 


Willie  the  two  states  discussed  above  apply  to  all  RRS  measurements,  there  exist 
two  sources  of  errors  which  are  unique  to  each  individual  transponder.  First  is  the  error 
due  to  R:i  trausponder  surveyed  position  uncertainty  (x,  y,  z  components  in  ECEF  frame), 
and  second  is  the  error  due  to  atmospheric  propagation  delays  between  the  user  and  each 
individual  transponder.  The  three  position  error  sources  are  well  modeled  by  random  bias 
states  (as  discussed  in  Chapter  Ill),  and  the  atmospheric  error  states  are  represented  by 
first  order  Markov  processes  (also  discussed  In  Chapter  III).  Then  for  each  transponder, 
four  states  are  used  to  define  error  sources  [31]: 


i, 

’  0 

0 

0 

0 

X, 

in 

.  _ 

0 

0 

0 

0 

y> 

►  4’  ^ 

i, 

0 

0 

0 

0 

wt! 

■ 

,  GR'ilmi 

0 

■ 

0 

0 

....  J. 

*  r  . 

Walmt 

(4.13) 
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The  subscript.  i  In  the  equation  above  denotes  the  i,h  RRS  transponder  and  r  =  700  see 
is  the  RRS  transponder  atmospheric  error  state  time  constant.  The  initial  conditions  for 
these  states  are  [31] : 

Xr,|M,fl/»n(^n)  —  0  (4.14) 

25  ft'2  0  0  0 

0  25/^  0  0 

0  0  25/tJ  0 

0  0  0  lOO(PPM)2 

and 

E  atm}-  0  (4.16) 

0  0  0  0 

E{**.u.*,ahn{t)'"*,u,i,nh,i(t  +  r)}  =  0  0  0  0  5(T)  (4.17) 

0  0  0  0 

.o  o  o  Sr). 

with  cr'fllm  ~  10~ ,n.  Once  again,  the  set  of  equations  above  apply  to  a  single  transponder. 
There  are  six  such  sets  of  equations  for  RRS  transponders  which  are  used  In  this  thesis. 
The  error  vector  is  specified  in  Appendix  A, 

4-5  Predicted  RRS  Measurement  Equation 

The  Kalman  filter  combines  range  measurements  [z]  that  are  generated  in  the  truth 
model  with  its  own  estimate  of  the  measurements  [2]  to  calculate  an  optimal  estimate  of 
the  state  vector.  The  Kalman  filter  update  equation  which  makes  use  of  the  actual  and 
estimated  measurements  is: 

x+  ~  x'  +  K(z  z  )  (4.18) 

where 

x+  =  filter  estimate  of  state  vector  just  after  measurement 

x~  =  filter  estimate  of  state  just  before  measurement 

K  =  Kalman  filter  gain  matrix 
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incoming  measurement  vector 

Kalman  filter  prediction  of  incoming  measurement  vector 


z 

r 


When  measurements  are  processed  sequentially,  the  measurement  vector  2  becomes 
a  series  of  scalar  range  measurements  z.  The  Kalman  filter  prediction  of  each  incoming 
measurement  is  also  in  the  form  of  sequential  scalar  ranges  measurements  z.  The  indicated 
positions  of  the  transponder  and  user  are  modeled  by: 


x,.  ss  xr  +  6x.r  (4.19) 

X,.  =  x„  -i  Sxv  (4,20) 


where 


X., 

ss 

surveyed  transponder  position  vector 

X,. 

ss 

true  transponder  position  vector 

6x,. 

s= 

true  error  in  surveyed  transponder  position 

x„ 

= 

indicated  user  position  vector  (LN-93) 

= 

true  user  position  vector 

Sx,t 

true  error  in  indicated  user  position 

Subtracting  the  Kalman  filter  position  error  estim.  .•  I>zr  and  from  the  applicable 
equations  yields  the  niter’s  best  estimates  of  the  transponder  and  user  positions: 


xr  -  SxT 

(4.21) 

x.f  +  Sxr  -  6x, 

(4,22) 

x7  +  enthvnliimvrror 

(4.23) 

xr  -  f>X, 

(4.24) 

k(;  6xr  -  fix. 

(4.25) 

x,,  + estimationerror,, 

(4.26) 
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Tlit*  Kalman  filter  estimate  of  the  upcoming  measurement  is -given  by i 

Ti  ®  Rn;s  -  Rr  (4.27) 

or 

Sz  ~  sj(x,,  -  ar)2  4-  (y,,  -  yr)'1  f  (?,,  -  £r)z  -  -  Of,r  (4.28) 

where 

Sr,  ffr,  tr  —  filter  estimate  of  transponder  position 

®u>  Vi/ »  =  filter  estimate  of  uier  position 


4.6  Tranaformaiion  0}  Uaer  Position  Error 

In  Section  4.3  it  is  assumed  that  user  position  and  position  errors  are  known  in  the 
ECEF  frame.  Tliia  section  describes  the  transformations  needed  to  process  measurements 
in  the  ECEF  frame  and  update  the  INS  position  error  states  in  error-angle  space  (the 
LN-93  error-model  space).  The  LN-93  Truth  Model  and  Error  Budget  document  [20] 
states  that  globally  valid  equations  which  transform  60, r,  60]n  and  60,  to  6<ps  S,\,  and  6a 
are  generally  quite  complicated  and  seldom  needed.  For  non-polar  regions,  the  following 
approximate  transformation  is  used: 


f  6<f> 

60u  sin  a  -  69 cos  a 

6X 

1  *  ’ 

{SOy  cos  a  +  60x  sin  a)  sec  <p 

6ot 

V 

60,  -  6X  sin  <p 

(4.29) 


Tli*  equation*  above  may  be  expreiaed  a*  a  (ingle  matrix  equation: 


'  6<j>  ' 

-  cos  a 

sin  a 

0 

0  ’ 

60,,.  ' 

’  69,  ' 

5A 

sin  a  sec  0 

cos  a  sec  0 

0 

0 

69v 

1 

69, 

» 

i  ::z 

c;; 

< 

8a 

-  sin  a  tan  0 

cos  a  tan  0 

1 

0 

so. 

69, 

.  6h 

0 

0 

0 

1.  _ 

6h 

6h 

where  Cft  i*  the  error-angle  space  to  navigation  error  space  transformation  matrix  defined 
in  Chapter  II.  If  user  position  is  expressed  in  the  ECEF  frame  as: 


-  f(0,A,a,/i)  = 


\ 

0 

(Re  +  h)  cos  0  sin  A  1 

0 

(  Re  +  h)  sin  0  j 

He  +  h 

(Rtf  + h)c os  0  cos  A  1 

(4.32) 


(where  C"  is  as  yet  undefined,  but  is  intended  to  transform  wander  azimuth  vectors  to 
navigation  frame  vectors),  then  user  position  may  also  be  expressed  as: 


Ixr 
y' 

zr 


(4.33) 


Equations  (4.32)  and  (4.33)  may  be  expanded  In  a  Taylor  series.  Because  this  is  a  linearized 
error  model,  it  is  assumed  that  the  higher  order  terms  (li.o.t.’s)  are  negligible  [3,  10,  22]. 
Then  the  series  becomes  x((  =  x,tt  +  Ax,,,  and  the  error  vector  term  may  be  written  as: 


’ 

’  50  ' 

6xK 

6y(’ 

fif (0,  A,  q,  h) 

ex,. 

< 

X  >1  mu 

5A 

6a 

6zr 

6h 

(4.34) 


where  the  user  position  vector  Is  expressed  in  terms  of  latitude,  longitude,  wander  angle, 
and  altitude  as: 


X,. 


0t  A,  a,  h 


(4.36) 
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and 

-(Rii  +  h)  tin  (Ri^  +  h)  cos ^ cos  A  0  cos  0  sin  A 

{Ri.  +  h)  coi0  0  0  sintf.  (4.30) 

-(Aw  +  h)  sin  ^  cob  A  -{Ri:  +  h) cot <f> tin X  0  cos^cosA 

Equations  (4.31),  (4.34),  and  (4.36)  are  combined  to  obtain  [11]: 

6x i<:  (Ric  4-  h)66j:  1 

6yv  <  -  L  '  (R,.;  +  h)66v  >  (4.37) 

Si"  Sh  j 

where 


Lu 

= 

sin  a  cos  A  4-  cos  a  sin  A  Bin  <p 

Lu 

= 

cos  a  cos  A  •-  sin  a  sin  A  sin  <j> 

i|:i 

as 

sin  A  cos  <f> 

L-i  i 

=1! 

-  cos  a  cos  </> 

Ln 

ss 

sin  a  cos  <j> 

Ln 

sin  <j> 

Lm 

=s 

-  sin  a  sin  A  +  cob  a  cob  A  sin  <j> 

L'vi 

=2 

-  cos  a  sin  A  -  sin  a  cos  A  sin  <j> 

Lxi 

- 

cos  A  cos  $ 

Rearranging  the  the  right-hand  side  of  Equation  (4.37)  results  in: 


( R,;;  +  h)60u 
-(«/..  +  h)69r 
fih 


(4.38) 
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The  matrix  [  C",  ]  li  the  transformation  matrix  which  converts  a  vector  from  the  wander 
azimuth  frame  to  the  EC£F  frame,  and  Is  written  as; 


a  v  turn 

xv 

=  Cl  ^  (4.39) 

. 

The  individual  elements  of  Cl  are  given  as  [20] 


Cu 

— 

cos  q  cos  A  -  sin  a  sin  A  sin  <j> 

C*ia 

= 

-  sin  a  cos  X  -  cos  a  sin  X  sin  0 

0,3 

= 

sin  X  cos  0 

Oa, 

= 

sin  a  cos  <f> 

Cn 

ES 

cos  a  cos  0 

Crs 

= 

sln0 

Cm 

= 

-  sin  A  cos  a  -  cos  A  sin  0  sin  a 

0*, 

= 

sin  A  sin  a  -  cos  A  sin  <j>  cos  a 

cm 

= 

cos  A  cos  <j> 

Then  the  filter  estimate  of  the  range  between  the  user  and  the  ground  transponder  Is: 

®  -  (xr  +  Sx.,.  -  6xr )  -  (x,j  +  ( c;;,  '|ffxjf  -  [  c;„  ]dx„)  -  SRa„„  -  fkllr  (4.40) 
—  I  R:  +  EH., |  +  SR, ,im  +  6Rt,r  (4.41) 

where: 

{R,.;  +  h)66u 

6x„  =  -  ~(R,,:  +  h)$04,  •  ,  and  «x„ 

Sh 

In  the  truth  model,  the  measurement  Is  given  by: 


Z  =  I  X  , .  -  X,  I  -  6ll„lw  -  SRbr 

(4.43) 

=  ]  R|  |  +  6R„i,„  +  SRi,v 

(4.44) 

(4.45) 
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()?,/.;  +  h)50u 

-(£,.; +  S)M,  •  (4.42) 

Sh 


y  . .  . 

j; 

i 

1  Then  the  residual,  formed  u  the  difference  between  Equation  (4.44)  and  (4.41),  results  in 

the  true  whole  range  magnitude  being  cancelled: 

A«  =  z  -  8  (4.46) 

i 

where  z  and  8  are  speciAed  in  Equations  (4.44)  and  (4.41).  It  is  the  “new  information"  con- 
|  tained  in  the  scalar  residual  which  is  scaled  by  time- varying  Kalman  Alter  gains  (discussed 

i  in  the  next  section)  to  update  the  Alter  state  estimates. 

4  >7  Kalman  Fitter  Gain 

The  Kalmar  Alter  gain  matrix  K(t)  used  in  Equation  (4.18)  determines  which  of  the 
elements  in  the  Alter  state  vector  are  affected  by  measurements  that  are  processed  by  the 
Alter.  The  magnitude  of  the  “correction”  to  individual  states  at  update  time  is  determined 
by  the  magnitude  of  the  elements  in  the  K(t)  matrix.  The  K (t)  matrix  is  computed  prior 
to  measurements  from  the  following  equation: 

K(t.)  =  P(tr)Hr(t,)[H(t,)P(t(-)Hr(t,)  +  ».(*()]“*  (4-47) 

where  the  H  matrix  is  the  Jacobian  gradient  of  the  noiseless  measurement  vector  [3]  and 
the  P  matrix  is  the  Alter-computed  covariance  at  time  (t~).  In  the  case  of  the  RRS 
range  measurements,  six  transponders  are  used  to  form  moasurements.  For  the  purpose 
of  software  efficiency,  the  six  measurements  are  assumed  to  occur  cyclically  at  one  second 
intervals  (as  opposed  to  actual  hardware  measurements  which  may  occur  aiynchronously). 
Thus,  a  complete  RRS  measurement  cycle  takes  six  seconds,  with  one  transponder  mea¬ 
surement  ocruring  every  second. 

The  elements  of  the  H  matrix  have  previously  been  derived  in  terms  of  the  ECIiF 
frame  measurement  equation.  However,  in  order  to  update  the  INS  error  states,  the  coeffi¬ 
cients  which  precede  the  SX„  terms  must  be  transformed  back  into  the  error-angle  space 
in  which  INS  error  model  computations  take  place.  Note  that  the  difference  measurement 
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equation  may  be  written  In  it*  functional  form  which  it  [ll]t 


T*  -  h(*u,  yin  i 

’/.(»  >  kri  Jr >  ^nfnn  6R hr )  +  V 

(4.48) 

If  the  expression  above  is  Written  as  a  Taylor  series  In  which  the  terms  above  first  order 

are  neglected  (3,  22],  the  result  is: 

1  t  . 1 

, ,,  y  §5*  MO 

'  +  '  W 

,  Ste. 

(4.49) 

.  _  *Y 

where:..  •  .  .  .[ 

(‘Jnom 

(‘)tium 

and 

(4.60) 

(•)  ss  (a(/,  p„, 

*i»  ®yi  Vi’*  Mlbr) 

(4.61) 

Tl^e  noiseless  part  of  fz  ia  the  original  h  function.  The  second  term  In  Equation  (4.49) 

‘  '  '  t  »  \  t 

ia  evaluated  by  taking  eight  partial  derivative*  a*  indicated  by  Equation  (4.60).  Only  the 
non-zero  partial1  derivatives  Indicated  by  Equation  (4.60)  are  discuaaed  below;  most  of 
the  elements  of  the  corresponding  H  row  will  also  be  zero  [11]  due  to  the  fact  that  the 
measurement  depends  upon  only1  eight  states  [represented  by  the  arguments  contained  In 
(-)  above].  Tho  non-  zero  elements  of  H  are  shown  in  the  following  row  vector: 


Ml 

hrv 


l)h_  Oh 
0~>r  IJfllaim 


where  the  flrat  three  terms  {after  transformation  back  into  wander  azimuth  coordinates) 
will  affect  (INS)  state*  1,  2,  and  10.  The  6Ri,r  term  ia  the  EES  “common”  bias  state 
which  affects  all  transponder  measurements  and  -('presents  the  /f(5lvv+  l)  element.  The 
remaining  four  terms  are  specific  to  the  transponder  which  is  being  “interrogated"  and 
affect  the  respective  ®,y,  z,  and  atm  error  states  In  the  EES  model.  As  an  example, 
for  the  98-state  model  presented  in  Chapter  VI  and  Appendix  F,  these  four  terms  are 
Jf(75),  1/(76),  JZ(77),  and  H{ 78)  for  a  measurement  from  EES  transponder  number 
one. 


4-16 


Noting  that  the  h{>)  function  li  linear  hi  the  error  terms  due  tc  equipment  Induced  Idm 
and  the  atmospheric  error,  the  applicable  row  of  the  H  matrix  may  be  written) 


»h  Oh  Oh  M.  .ML  ML  i  i 
uvtr  TfSv  o»v  1  1 


Evaluating  the  partial  derivatives  of  h(-)  with  respect  to  the  transponder  position  error 
states  reveals  that  these  terms  are  simply  the  components  of  the  unit  linc*of-slght  from  the 
user  to  the  transponder,  expressed  in  the  ECEF  frame  [11].  The  H  row  may  be  written 
as: 

ft  ft  ft  VJOS,  yi05„  ULOS,  1  1 

The  following  section  is  adapted  from  a  development  in  [11].  Referring  to  Equa¬ 
tion  (4.40),  the  flrtt  order  term  (evaluated  with  respect  to  the  user  position  error  states 
only)  may  be  expressed  In  the  ECEF  frame  as: 


H(>)  St)  = 

(')ruMM 


dh  ,  « 

0*;,  Hutu  *u<' 


(4.B2) 


Again,  the  arguments  represented  by  (•)  in  the  equation  above  are  now  limited  to  the 
user  poaiiion  slates  [  y„,  s,,  ].  A  3-by-3  Identity  matrix  may  be  inserted  between  the 

“partial”  term  and  the  error  vector  without  changing  the  result: 


ff(.)  I.  =  gjr 


(4.53) 


Now  replace  the  identity  matrix  by  two  direction  cosine  matrices  [which  when  multiplied 
together  would  return  the  Identity  matrix]: 


H(.)  8.  =  c;„  cr  Sxv'r 

\ ' )  HfJMI 


(4.54) 


Evaluating  the  partial  derivative  and  multiplying  out  the  product  C )!'  Sx^r  leaves: 


H(‘)  St) 

( *) ....... 


-ULOS,-  -  VLOS„  -  ULOSs\  C'tl.  txv* 


(4.55) 


4-16 


If  the  row  vector  and  the  DCM  are  multiplied  together,  a  row  vector  i»  obtained.  Tor 
convenience,  the  remit  ii  defined  as  [  a  b  cj.  The  right  half  of  the  product  can  be  rewritten 
uiing  the  results  of  Equation  (4.42): 


*(•) 


ix  = 


(•)« 


a  6  c] 


(Re  +  h)S9v 
-(Re  +  h)60, 

Sh 


(4.66) 


Finally,  the  expression  above  may  be  rewritten  to  yield: 


*(•) 


6x 


(')nnm 


[  -b(Ri,:  +  h)  a(R,,:  +  h)  c] 


’  66,  ' 

1  seu 

6h 


(4.67) 


Then  the  first  three  H  elements  in  the  row  are  - b(R /.;  +  h),  a(R/;  +  /i),  and  c. 
Equations  for  all  other  non-sero  elements  in  H  are  developed  above.  These  equations  are 
programmed  into  MSOFE  for  the  RRS  measurement  model. 


4-8  Transponder  Physical  Locations 

The  six  transponders  which  are  modeled  in  the  NRS  composite  error  model  are  those 
which  are  on  and  near  Holloman  AFB,  NM.  The  locations  of  these  transponders  have  been 
surveyed  (typically  to  within  6  feet  in  each  of  the  3  axes),  and  the  coordinates  are  given 
in  Table  4.1. 


Table  4.1.  RRS  Transponder  Locations  [31] 


Transponder  ID 

Latitude 

Longitude 

Altitude 

005 

33°  0V  36.1472”' 

-  i06u08720J4047r 

4339  ft 

102 

32°  65'  68.5986" 

- 106'’  08'  50,3339" 

4074  ft 

181 

33°  44'  58.035" 

- 100°  22'  14.630" 

7932  ft 

211 

33°  17'  55.999" 

- 106°  3T  44.311" 

8842  ft 

212 

32°  47'  16.418" 

-105°  49' 16.474" 

9202  ft 

216 

32°  42'  12.235" 

-106°  07' 38.907" 

4481  ft 
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4 >9  Summary 

Tills  chapter  present*  the  basic  concept*  related  to  the  RES  transponder  syitam  and 
include*  the  ERS  transponder  system  error  model  equations  a*  well.  In  addition,  the  range 
measurement  model  equation  is  developed  in  detail.  The  conventional  difference  measure¬ 
ment  approach  it  presented  and  discussed,  chiefly  to  point  out  that  the  method  results  in 
the  whole-valued  quantity  being  canceled  in  the  "delta-measurement”.  Finally,  the  calcu¬ 
lations  required  to  generate  H  matrix  elements  for  a  typical  EES  range  measurement  sure 
demonstrated. 
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V.  Global  Positioning  System  Model  and  Space  Vehicle  Orbit  Calculations 
5.1  lntrodur.tion 

The  Global  Portioning  Syitem  (GPS)  li  designed  to  be  a  highly  accurate,  stand-alone 
navigation  system.  However,  for  this  research,  GPS  is  used  as  a  subsystem  to  improve  the 
navigation  solution  of  the  LN-93  based  NRS.  In  a  manner  somewhat  reminiscent  of  the 
HRS  transponder  system  discussed  in  Chapter  IV,  GPS  navigation  Information  1s  obtained 
from  electromagnetic  signal  propagation  through  the  media  (space  and  atmosphere)  be- 
tween  the  user  (NRS)  and  each  of  the  space  vehicles  (SVs)  which  the  user  "locks”  Into  a 
reception  channel  of  the  GPS  receiver.  Navigation  Information  is  obtained  by  receiving 
GPS  SV  ephemerls  data  which  are  broadcast  continuously  from  each  active  (“locked- on”) 
SV,  correlating  the  phase  of  the  signal  with  a  matching  signal  in  the  GPS  receiver,  and 
correcting  for  known  error  sources  to  produce  a  highly  accurate  range  estimate  between 
the  user  and  each  SV  which  is  monitored.  Although  not  used  in  this  thesis,  range-rate 
information  may  be  obtained  from  GPS  ephemerls  in  a  similar  manner.  As  in  RES,  GPS 
range  (and,  if  used,  range-rate)  measurements  make  refinements  to  the  NRS  navigation 
solution  possible. 

Figure  6.1  depicts  the  addition  of  the  GPS  subsystem  to  the  RRS  and  baro-aided 
INS  diagram  presented  In  Chapter  IV.  This  is  the  final  module  addition  and  completes  the 
full  NRS  model. 

A  dynamic  error  model  for  the  GPS  system  is  developed  in  the  Solomon  thesis  [33] 
and  revised  in  [11].  Portions  of  those  bodies  of  work  are  summarized  in  this  chapter. 
However,  substantial  changes  to  the  basic  GPS  model  are  made.  In  the  references  cited 
[11, 33]  a  simplified  GPS  model  was  assumed,  It  consisted  of  four  stationary  space  vehicles 
(SVs)  and  did  not  perform  geometry  optimization  calculations,  In  this  thesis,  a  24-SV 
“optimal”  constellation  based  on  a  paper  by  Green  j  1 2]  is  modeled.  The  new  model  includes 
orbital  calculations  for  all  SVs,  and  simulates  GPS  receiver  operation  as  well.  These 
enhancements  are  discussed  in  detail  in  subsequent  sections  of  this  chapter.  As  in  the 
case  of  RRS  measurements,  the  lever-arm  effect  is  important  to  consider  in  actual  GPS 
hardware  applications.  Once  again,  the  lever-arm  effect  is  avoided  by  assuming  that  the 


6-1 


Figure  5.1.  NRS:  LN-93  INS  Aided  by  Baro-  Altimeter,  RRS,  and  GPS 

INS  and  the  GPS  antenna  are  collocated.  The  assumption  of  collocation  is  made  In  the 
interest  of  generality  and  NRS  software  efficiency.  The  measurement  model  equations  for 
the  GPS  system  follow  a  parallel  development  to  that  of  the  RRS  measurement  model  in 
Chapter  IV. 

5.2  GPS  Range  Measurements 

Once  again,  in  ARS  or  the  proposed  NRS,  GPS  range  measurements  aid  in  estimating 
position  errors  of  the  reference  system.  The  GPS  range  measurement  is  derived  from 
decoding  ephemeris  data  which  are  broadcast  continuously  by  each  active  SV.  The  mer’n 
GPS  receiver  (now  considered  to  be  a  subsystem  in  NRS)  processes  signals  which  are 
received  from  the  GPS  SVs  to  determine  pseudo-range  between  the  user  and  the  SV.  The 
range  measurement  thus  obtained  is  corrupted  by  several  error  sources  which  must  be 
determined  and  compensated. 

In  its  simplest  form,  a  range  measurement  between  a  single  GPS  SV  and  the  user 
(in  this  case,  NRS)  may  be  determined  as  the  product  of  propagation  speed  of  the  electro¬ 
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tirngtirMr  (EM)  signal  and  elapsed  time  during  such  propagation.  Stated  mathematically, 
the  range  (typically  called  pseudo-range  due  to  inherent  inaccuracies  in  the  estimate)  is 
given  by: 

Rv  =  c  ■  tr  (6.1) 

where  R,,  Is  pseudo-range,  c  is  the  speed  of  light  in  vacuum,  and  t,.  represents  the  elapsed 
time  for  transit  of  the  EM  signal. 

However,  two  major  problems  exist  in  using  this  simple  definition  for  range.  First, 
the  EM  signal  Is  not  propagated  entirely  in  vacuum.  The  signal  originates  In  space  (where 
the  assumption  of  vacuum  Is  acceptable),  but  must  subsequently  propagate  some  distance 
through  the  earth’s  atmosphere  as  well.  Naturally,  the  signal  delay  introduced  by  atmo¬ 
spheric  propagtion  must  be  taken  Into  account.  Second,  in  order  to  preserve  any  hope  of 
accurately  determining  range,  it  is  critical  to  determine  the  EM  signal  transit  time  with 
an  extremely  high  degree  of  accuracy.  Recalling  that  light  (or  any  EM  signal)  propagates 
on  the  order  of  3  x  108m/s,  it  is  readily  apparent  that  even  a  very  small  error  In  deter¬ 
mining  the  EM  signal  propagation  time  can  wreak  havoc  on  attempts  to  use  pseudo-range 
Information  to  Improve  the  navigation  system  solution  In  NRS  (or  any  other  such  system). 

As  a  consequence  of  the  concerns  above,  it  is  Imperative  to  develop  a  much  higher 
fidelity  model  for  range  estimation.  A  typical  GPS  receiver  models  the  range  between  the 
user  and  space  vehicles  with  the  following  equation  (32|: 

®'f ~  d"  d R,  I  +  d  R; , -f'  6  R,(!;i  T  dH.yru  +  +  V  ( d .  2 ) 


where 


GPS  pseudo-range  measurement,  from  SV  to  user 

a, 

True  range,  from  SV  to  user 

8  Hr/  « 

Range  error  due  to  codeloop  error 

d®. irop  — 

Range  error  due  to  tropospheric  delay 

d  R|i>n  — 

Range  error  due  to  ionospheric  delay 

d  R-Srik  ~ 

Range  error  due  to  SV  dock  error 

dH. llrlk  “ 

Range  error  due  to  User  clock  error 

V  = 

zero-mean  white  Gaussian  measurement  noise 

As  In  the  comparable  equation  for  RRS,  the  GPS  pseudo-range  equation  above  In¬ 
cludes  the  true  range  [which  can  never  be  known  exactly]  along  with  terms  which  reflect 
sources  of  error  and  uncertainty  Inherent  to  GPS  range  measurements.  GPS  error  sources 
and  models  are  discussed  later. 


5.3  GPS  Pseudo- Range  Calculation  From  INS  Data 

As  in  the  RRS  subsystem,  It  is  desirable  to  formulate  a  difference  measurement  In  the 
GPS  model  as  well.  Once  again,  two  sources  of  range  information  must  be  obtained.  Like 
the  RRS  case  presented  in  Chapter  IV,  the  first  source  is  the  range  measurement  will  i 
comes  from  the  GPS  subsystem  and  which  is  modeled  by  Equation  (5.1).  The  second 
range  estimate  is  constructed  by  differencing  INS -indicated  position  and  SV  (broadcast) 
positions  to  calculate  the  range.  (Note  that  tile  indicated  INS  position  and  SV  position 
contain  uncertainties  which  must  be  considered.)  In  Chapter  IV,  the  user  (INS)  position 
is  represented  by  an  Jl1  vector  expressed  in  the  Litton  ECEF  which  is  repeated  below  for 
convenience.  SV  position,  Xv,  is  represented  in  like  fashion. 
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Then  user  to  SV  pseudo-range  may  be  calculated  as: 


► 

v. 

1  ' 

•v 

xs 

*,ss  = 

x„  -  X, 

< 

y. 

►  -  i 

Vs 

.  , 

.  **  . 

Equation  (6.4)  may  also  be  written  in  the  equivalent,  form: 


(5.3) 


(5.4) 


R,NS  -  \[Wc  -  +  (.'/,  -  y.  )-’  +  (*,;  -  Z.sY 


(5.5) 
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Invoking  perturbation  theory  [3,  22],  the  equation  above  is  written  as  a  first-order  Taylor 
series  to  approximate  the  (INS  derived)  calculation  of  user  to  SV  range.  The  truncated 
(to  first-order)  series  is  expressed  as: 


=  Rt  + 

+ 


0fl,w,(X„X„) 


0X, 


(X,,Xt;) 


«xs 


0X„ 


(Xs..x„)M 


sx, 


(6.6) 


When  equation  (5.5)  is  substituted  into  Equation  (5.6)  and  the  partial  derivatives  evalu¬ 
ated,  the  INS-derived  pseudo-range  approximation  becomes: 
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4*,  (5.7) 


Now  the  GPS  pseudo-range  difference  measurement  is  formed  as: 


6z  —  R,n,s  X 
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-  [ljfifl,.,  -  [l}6Rln>ll  -  +  #  (5-8) 


As  noted  in  the  RRS  case,  the  true  whole- valued  range  ( R ,)  formerly  present  in  the 
individual  pseudo-range  representations  (/t,iV..  and  R, is  cancelled  in  the  differencing 
operation.  The  bracketed  coefficients  in  the  equation  above  will  again  be  used  in  the  H 
matrix  development  discussed  in  a  later  section. 

In  order  to  form  the  GPS  difference  measurement,  it  is  assumed  above  that  the  INS 
and  SV  coordinates  are  expressed  in  the  same  frame.  However,  the  assumption  and  reality 
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arc  at  odds  once  again.  Tlie  SV  whole-valued  quantities  are  generally  expressed  In  either 
their  own  orbital  frame  ox  possibly  in  the  Inertial  frame  (as  defined  by  Brltting  [3]).  In 
Chapter  IV,  a  transformation  Is  presented  that  takes  the  whole- valued  INS  position  vector 
triplet  [latitude,  longitude,  altitude]  and  converts  it  to  the  (Litton)  ECEF  frame  as  shown 
below: 


(5.9) 


' 

*</ 

V 

> 

(RN  4-  h )  cos  i^sinA 

Vv 

>  -  < 

[(JUV)(1  -  ea)  +  *]  sM 

*" 

(RN  +  h)  cos  0  cos  A 

Equation  (6.9)  depends  on  knowing  the  user  position  in  terms  of  latitude,  longi¬ 
tude,  and  altitude.  The  discussion  in  Chapter  III  presents  the  LN-93  error-angle  vector 
[  6$,  ]  which  may  be  transformed  into  latitude,  longitude,  and  altitude  error 

space.  The  approach  for  obtaining  the  INS  position  error  in  ECEF  coordinates  (from  the 
f  SL ,  SX,  6h  ]  triplet)  is  presented  in  detail  In  Section  4.6. 


The  SV  positions  are  routinely  defined  In  terms  of  orbital  parameters  which  may 
then  be  defined  in  terms  of  inertial  space  coordinates.  Additionally,  the  “flight”  of  the 
SVs  in  their  orbital  planes  and  the  effect  of  earth  rotation  must  be  taken  Into  account.  SV 
positions  are  (carefully)  transformed  to  the  Litton  ECEF  as  discussed  below. 


Each  of  the  24  SVs  Is  assigned  an  initial  position  (refer  to  Sections  6.8  and  6.9)  based 
on  optimisation  of  global  GPS  coverage  using  the  24  SV  constellation  [12].  Four  SVs  are 
assigned  to  each  of  six  orbital  planes;  each  orbital  plane  is  inclinod  at  56°  with  respect  to 
the  inertial  space  p,  plane.  The  six  lines  of  nodes  [also  called  longitudes  of  ascending 
nodes  or  LANs]  are  equally  spaced  around  the  circumference  of  the  earth. 

A  single  (representative)  orbital  plane  is  shown  in  Figure  5,3.  The  orthogonal  axes  [/, 
J ,  R]  are  identical  to  the  [»,,  p,,  *,■]  axes  defined  In  the  Brltting  inertial  frame  discussed  in 
Chapter  II,  The  LAN  is  the  point  at  which  an  SV  in  this  orbit  crosses  the  equator  traveling 
from  south  to  north,  and  is  represented  in  the  figure  as  n.  Additionally,  the  angle  between 
the  LAN  and  the  inertial  frame  axis  x,  is  designated  11.  The  angular  momentum  vector 
( h  in  the  figure)  is  normal  to  the  orbital  plane,  ami  the  angle  between  h  and  K  is  called 
the  angle  of  inclination  and  is  given  the  symbol  (3  for  this  document.  Besides  the  orbital 


Figure  5.2.  Orbital  Path  for  a  Typical  SV 


plane  angle  of  inclination  (0)  and  the  LAN  angular  displacement  from  the  a,  axis  (ft),  it 
is  also  necessary  to  know  the  position  of  each  SV  in  its  orbital  path  at  some  point  in  time, 
say  to.  If  initial  conditions  are  known.,  it  Is  a  straightforward  calculation  to  determine  SV 
position  in  the  orbital  plane  at  any  time  subsequent  to  to. 

Although  the  figure  depicts  elliptical  orbits,  the  model  assumed  in  this  research  is 
that  of  circular  orbits.  This  assumption  is  based  on  the  optimal  SV  constellation  described 
in  [12].  In  the  optimal  constellation  (which  is  currently  being  Implemented  in  practice), 
SV  orbital  semi-major  axes  are  greater  (26,609  km)  than  previous  constellations,  which 
were  on  the  order  of  26,560  km.  The  net  effects  are  (a)  a  more  circular  orbit,  (b)  improved 
coverage,  and  (c)  reduced  elllpticity  in  the  orbital  paths.  While  the  ellipticity  1b  tabularized 
as  zero  in  the  Green  paper  [12],  the  actual  ellipticity  is  non-zero,  but  its  effect  is  sufficiently 
small  to  be  neglected  for  the  purpose  of  this  research.  If  later  investigation  proves  to  the 
contrary,  then  a  change  to  the  SV  orbit  calculations  may  be  made  to  account  for  ellipticity. 
At  this  point  the  SV  orbital  model  is  developed. 

Consider  for  a  moment  that  each  GPS  SV  travels  in  a  (nearly)  circular  path  in  a 
single  plane,  Irrespective  of  the  orientation  of  that  plane.  In  a  situation  such  as  this,  the  SV 
position  is  depicted  in  Figure  5.3  as  a  function  of  orbital  radius  and  an  angle  a.  The  angle 
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LAN 


Figure  5.3.  Planar  Representation  of  SV  Orbital  Path 

alpha  is  defined  (conveniently)  with  respect  to  the  LAN,  and  is  described  mathematically 
as: 

a(t)  =  Q„  4-  ui„  •  t  (5.10) 

Mow,  define  an  orthogonal  frame  in  R:1  which  has  its  origin  coincident  with  the  orbital 
plane  origin  shown  in  Figure  6.3.  The  new  frame  x„  axis  is  colinear  with  the  orbital  plane 
LAN;  the  new  yt,  axis  is  90  degrees  counterclockwise  from  and  coplanar  with  ®ut  and  z tJ  is 
pointing  out  of  the  page.  The  newly  defined  frame  is  shown  in  Figure  5.3.  Note  that  the 
SV  position  as  defined  in  this  frame  is  still  “two-dimensional,"  in  that  the  zu  component 
of  the  SV  position  vector  is  identically  zero. 

Keeping  in  mind  that  a  is  time- varying,  the  SV  position  in  the  orthogonal  orbital 
frame  depicted  in  Figure  5.3  is  given  as: 

R  CUH  (> 

R.  sin  a  '  (5.11) 

0 

Now  it  remains  to  rotate  the  SV  position  expressed  in  the  orbital  frame  above  into 
the  inertial  frame,  taking  care  to  account  for  the  orbital  plane  angle  of  inclination  and 


5-8 


xot  LAN 

Figure  5.4.  Orthogonal  Orbital  Frame  Definition 

the  angular  displacement  of  the  LAN  from  the  inertial  frame  a,-  axis.  Two  rotations  are 
needed  to  accomplish  the  transformation.  First,  the  orbital  frame  is  rotated  clockwise 
about  its  own  axis  through  an  angle  of  55°  (the  orbital  plane  angle  of  inclination,  0). 
Second,  the  orbital  frame  is  rotated  counterclockwise  about  the  z„  axis  through  an  angle 
equal  to  the  original  angle  from  the  a*  axis  to  the  line  of  ascending  nodes  (LAN),  fi.  After 
these  two  rotations,  the  orthogonal  orbital  frame  is  coincident  with  the  inertial  frame. 
Thus,  an  SV  position  vector  expressed  in  the  orthogonal  orbital  frame  is  transformed  to 
the  inertial  frame  via: 

' 

cosfl  -cos/3sinn  sin 0  cos  ft  xu 

sinH  sln/Jcosfl  -sin/3cosfi  >  <  yu  (6.12) 

0  sin/3  cos/3 

/V  / 

where  0  is  the  orbital  plan  angle  of  inclination  (55"),  and  n  is  previously  defined.  The 
result  thus  obtained  is  a  step  closer  to  getting  SV  position  vectors  in  terms  of  ECEF 
coordinates.  Two  more  steps  must  be  taken.  First  a  transform  from  inertial  space  to  the 
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(Brit l  ing)  ECEF  ii  given  by  [3]t 


I' 

cos(w,«t)  sin(u>,„t)  0 
-  sin(wlf;t)  eos(u>„.<)  0 
0  0  1 


(5.13) 


where  ui(fl  is  the  earth’s  rotation  rate  and  t  is  the  elapsed  time  since  fa.  The  same  value 
for  t  is  used  here  and  in  the  defining  equation  for  the  SV  orbital  angle  a.  Finally,  the 
difference  between  the  Britting  ECEF  and  Litton  ECEF  must  be  taken  into  account  to 
complete  the  transform.  Using  the  matrix  specified  in  Chapter  II,  SV  positions  in  Litton 
ECEF  terms  are  given  by: 


«r 

Vv 


C\  •  C't> 


t.illon 


•  Cl  .  |  Vo 


(5.14) 


At  last  the  S  V  position  is  expressed  in  the  desired  frame.  However,  recall  that  the  INS 
position  vector  above  is  expressed  in  terms  of  latitude,  longitude,  and  altitude.  Once  again, 
the  INS  position  must  be  expressed  in  terms  of  its  error-angle  states  but  coordinatized  in 
the  Litton  ECEF.  The  approach  is  identical  to  that  explained  in  Chapter  IV  and  is  omitted 
here. 


5.4  GPS  Error-State  Model  Equations 

The  GPS  error  state  vector  is  composed  of  30  elements  (shown  in  Table  A.6,  Ap¬ 
pendix  A).  The  GPS  states  occupy  the  thirty  “uppermost”  states  in  the  NRS  error  state 
model.  The  first  two  GPS  states  model  the  user  s«t  [receiver]  clock  bins  and  drift  errors, 
respectively.  The  error  state  model  equation  for  these  states  1b  [11,  33]: 

0  l  I  ZlrU;, 

0  0  |  X)  YU-,,,. 

where 
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viicU-t,  =  range  equivalent  of  user  set  clock  bias 

BUctkj,  =  velocity  equivalent  of  user  set  clock  drift 


(5.16) 


(5.17) 

Note  the  large  uncertainties  associated  with  the  user  clock  states.  Until  the  user  clock 
error  is  determined,  it  Is  the  single  largest  source  of  error  In  GPS  range  measurements. 
While  the  two  states  discussed  above  apply  to  all  GPS  measurements,  there  exist  five 
sources  of  errors  which  are  unique  to  each  Individual  SV. 

One  error  source  specific  to  each  SV  is  codeloop  range  quantization  error.  At  the 
heart  of  any  GPS  receiver  exists  a  pair  of  interacting  tracking  loops  [21],  One  of  these 
loops,  the  “code  tracking  loop”  is  the  source  of  pseudo-range  estimation  error  which  is 
modeled  as  a  first  order  Markov  process  [21]  with  an  exponential  autocorrelation  function. 
Other  significant  error  sources  Include  the  tropospheric  and  ionspheric  propagation  delays. 
Both  of  these  error  sources  are  identified  and  corrected  to  a  large  degree  by  the  GPS 
receiver.  However  the  uncompensated  error  contribution  of  these  error  sources  is  still 
significant.  Both  of  these  are  also  modeled  as  first  order  Markov  processes  (with  different 
time  constants).  Still  other  sources  of  error  which  must  be  included  in  the  GPS  model  are 
SV  dock  error  and  SV  f?:l  position  error,  each  of  which  is  treated  as  a  random  bias  state. 
The  reader  should  note  the  strong  possibility  for  observability  problems  in  a  model  such  as 
tills.  In  this  case,  three  position  error  states  are  used,  but  the  measurements  provide  new 
information  only  along  the  linc-oj-aight  vector  between  the  user  and  SV.  All  error  sources 
discussed  above  which  are  unique  to  each  SV  are  included  in  an  error  state  vector  shown 


The  Initial  state  estimates  and  covariances  for  these  states  are  [31]: 
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and 


(5.20) 


E{  w  «J\v(0WG’f\V  (*  +  r)} 


0.5  0  0  0  0  0  0 

0  0.004  0  0  0  0  0 

0  0  0.004  0  0  0  0 

0  0  0  0  0  0  0 

0  0  0  0  0  0  0 

0  0  0  0  0  0  0 

0  U  0  0  0  0  0 


6(r) 


(5.21) 


Once  again,  the  set  of  equations  above  apply  to  u  single  SV.  There  are  four  such  sets 
of  matrix  equations  for  GPS  SV  errors  modeled  in  this  thesis.  The  error-state  vector  is 
completely  specified  in  Appendix  A. 


6-12 


Figure  5.5.  Two-Dimensional  Representation  of  User  and  SV  Position  Errors. 

5.5  Predicted  GPS  Measurement  Equation 

The  Kalmau  filter  combines  range  measurements  [2]  that  are  generated  in  the  GPS 
truth  model  with  its  own  prediction  of  the  measurements  [8“]  to  calculate  an  optimal 
estimate  of  the  state  vector.  The  general  Kalman  filter  update  equation  [22]  defined  In 
Chapter  IV  is  repeated  belowt 


S+  *i“  +  K(i-8“)  (5.22) 

When  GPS  measurements  are  processed  sequentially,  the  measurement  vector  2  be¬ 
comes  a  series  of  scalar  range  measurements  z.  The  Kalman  filter  prediction  of  each 
incoming  measurement  is  also  in  the  form  of  sequential  scalar  ranges  measurements  z~ . 
The  indicated  positions  of  the  SV  and  user  are  depicted  in  two  dimensions  in  Figure  5.5. 
The  position  errors  depicted  in  Figure  5.5  are  extended  into  R3  space  and  are  represented 
mathematically  by: 

i,  =  x,  +  i'x,  (5.23) 

x(;  =  x,.  +  fxr  (5.24) 

where 
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za 

indicated  SV  position  vector 

= 

true  SV  position  vector 

6xs. 

= 

true  error  in  indicated  SV  poiltlon 

*u 

== 

indicated  user  position  vector  (LN-B3) 

- 

true  user  position  vector 

Sxu 

= 

true  error  in  indicated  user  position 

Subtracting  the  Kalman  Alter  poiltlon  error  eitlmates  Sxs  and  Sxu  from  the  applicable 
equations  yield*  the  Alter'*  belt  eitimate*  of  the  SV  and  u*er  poiltlon*: 


*V  -  6ks. 

(6.26) 

x,  +  Sxs  -  dxs. 

(5.26) 

x,.  -f  error x 

(6.27) 

x„  -  6x~ 

(6.28) 

x„  +  6x„  -  6x„ 

(5.29) 

x(J  -f  error,. 

(5.30) 

The  Kalman  Alter  estimate  of  the  upcoming  measurement  is  given  by: 

=  h;NS  -  R;  (6.31) 

or 

-  i.v)2  +  (Vo  -  y,  )2  +  (*U  -  zj2  -  SRrl  ~  6Rtro,>  -  b%on  “  6&SM  (5.32) 

where 


2„,  y(l,  ?„  --  filter  estimate  of  user  position 

2y,  y.¥,  *!,  ’■=  filter  estimate  of  SV  position 
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ff.l?  OPS  Measurement  Estimation 

In  Section  5.3  it  Is  assumed  that  uaer  position  and  position  errors  are  expressed 
in  the  ECEF  frame.  This  section  describes  the  transformations  needed  to  process  OPS 
measurements  in  the  ECEF  frame  and  update  the  INS  position  error  states  in  error-angle 
space  (the  LN-93  error-model  space).  Using  the  exact  same  approach  presented  in  Chapter 
IV,  the  INS  position  is  written  as: 


u  I  v 

Yv  |  =  OS,  Yv 


(5.33) 


Consequently,  the  filter  estimate  of  range  from  user  to  SV  may  be  expressed  as: 


z~  -  (xs  +  6xs  -  6xs~ )  -  (x„  +  [  Crw  ]dx't;  -  [  Crw  ]fx~v  ) 

—  —  —  ^Srlk 


(5.34) 

(5.35) 


where: 


(R/i  +  h)6fff/ 

tfx|,  =  <  -.(Jt,.;  +  h)60,r  ,  and  Sx~' 


(RE  +  h)n-u 
-(/?,.;  +  h)ST,r 


(5.30) 


In  the  truth  model,  the  measurement  is  given  by: 


z  —  |  x„  xs.|  —  6Rci  SR) f„f)  —  SRion  —  SRsdk  ~~  GRudis 


(5.37) 

(6.38) 


In  the  residual  below,  the  true  whole  range  magnitudes  are  again  cancelled: 


Az  =  z  --  7. 


(5.39) 
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It.  Is  the  “new  information”  contained  in  the  scalar  residual  which  is  scaled  by  time- varying 
Kalman  filter  gains  (discussed  in  the  next  section)  to  update  the  filter  state  estimates. 

5.7  GPS  Measurement  Matrix 

The  observation  matrix  H(t)  elements  associated  with  the  OPS  measurements  must 
be  obtained  in  order  to  complete  the  Kalman  filter  update  process.  Four  OPS  measure¬ 
ments  are  assumed  to  occur  simultaneously  at  0.25 Hz.  Thus,  a  complete  GPS  measure¬ 
ment  cycle  takes  four  seconds,  with  all  four  SV  measurements  being  processed  sequentially 
at  update  times. 

The  GPS  elements  of  the  H  matrix  have  previously  been  derived  in  terms  of  the 
ECEF  frame  measurement  equation.  As  demonstrated  in  Chapter  IV,  the  coefficients 
which  precede  the  terms  in  the  measurement  difference  equation  are  in  fact  the  elements 
of  the  H  row  for  the  measurement  being  processed.  Referring  to  Equation  (6.8),  the  first 
three  bracketed  coefficients  map  into  the  user  position  states,  the  next  three  coefficients 
map  into  the  SV  position  states,  and  the  last  four  (ones)  map  into  the  remaining  GPS 
error  states. 

For  example,  if  SV|  is  used  to  form  a  GPS  range  measurement,  coefficients  1,  2,  and 
3  in  Equation  (5.8)  become  elements  1,  2,  and  10  in  the  H  row.  These  clement  numbers 
correspond  to  the  Jt:i  position  errors  for  the  user  INS.  Similarly,  coefficients  4,  5,  and  6 
in  Equation  (5.8)  become  elements  126,  127,  and  128  (corresponding  to  the  position  error 
states  for  SV  number  four).  The  remaining  (unity)  coefficients  map  into  the  remaining 
GPS  error  states;  continuing  with  SV  4  as  the  example,  elements  122,  123,  124,  and  125 
In  the  H  row  are  set  to  unity.  Element  99  in  the  H  row  is  unity  for  oil  four  SVs  because 
the  user  clock  affects  all  GPS  measurements  by  the  same  magnitude. 

5.8  5V  Measurement  Set  Selection 

Although  only  four  SVs  we  used  for  measurements  at  any  given  time,  positions  for 
all  24  SVs  in  the  constellation  are  continuously  updated  in  the  truth  model  for  simula¬ 
tion  purposes.  This  is  necessary  in  order  to  emulate  the  function  of  an  operational  GPS 
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Figure  5.8.  Optimal  SV  Set  Geometry  [4] 


receiver,  Ai  previously  noted,  the  GPS  receiver  select#  the  set  of  four  SV#  from  all  possi¬ 
ble  combinations  of  SV#  in  view  to  provide  the  best  possible  geometry  for  enhancing  the 
navigation  nolution. 

When  SV  and  user  positions  are  expressed  in  a  common  reference  frurne,  the  “correct” 
group  (or  set)  of  four  SVs  must  be  chosen  from  all  possible  combinations  of  SVs  In  view. 
The  approach  is  to  solve  an  Iterative  algorithm  [27]  which  determines  the  optimal  SV  set 
based  on  geometric  consideration#.  Ideally,  the  GPS  receiver  selects  a  set  of  SVs  arranged 
such  that  one  SV  is  directly  overhead  with  respect  to  the  user;  the  remaining  three  SVs 
are  as  low  on  the  horizon  as  possible  (while  still  permitting  clear  reception  of  their  signals 
by  the  user),  and  they  are  spaced  120°  apart  in  azimuth  as  seen  by  the  user. 

However,  ideal  geometry  is  seldom  (If  ever)  achieved.  Therefore,  the  goal  is  to  deter¬ 
mine  which  set  of  four  SVs  offers  the  best  geometric  configuration  available  at  measurement 
times.  The  algorithm  used  solves  the  problem  by  maximizing  the  volume  of  a  tetrahedron 
defined  by  the  intersection  of  unit  vectors  (from  user  to  SV)  with  a  hemispherical  surface 
above  the  user.  The  concept  is  illustrated  in  Figure  5.6.  While  it  is  relatively  simple  for  a 
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GPS  receiver  to  determine  whether  or  not  a  particular  SV  is  in  view  (either  the  SV  signal 
is  present,  or  It  is  not],  life  Is  not  quite  so  simple  in  the  simulation  environment. 

Because  signals  from  SVs  must  be  simulated,  it  is  necessary  to  determine  which 
SVs  are  “in  view”  (meaning  that  their  signal  is  not  “masked”  from  the  user  by  physical 
obstructions  such  as  the  earth).  Therefore,  once  SV  and  user  positions  are  determined 
in  a  common  frame  at  any  given  time  of  Interest  (say  a  measurement  time  t, ),  then  the 
elevation  of  each  SV  (with  respect  to  the  user)  must  also  be  determined.  The  method  is  to 
compute  the  elevation  angle  based  on  the  unit  line  of  sight  vector  from  the  user  to  each  of 
the  SVs  in  the  optimal  constellation  at  each  measurement  time.  If  the  computed  elevation 
angle  is  greater  than  a  predetermined  mask  angie  (5°  for  this  research),  then  the  SV  is  “in 
view”  to  the  user. 

The  optimal  SV  constellation  Is  arranged  such  that,  depending  on  UBer  position  and 
time,  a  minimum  of  five  and  a  maximum  of  eleven  SVs  may  be  in  view  at  any  given  time. 
Further  consider  that  any  four  SVs  can  form  a  measurement  set.  Then  the  number  of 
combinations  from  which  the  best  set  is  selected  at  each  measurement  time  may  range  from 
5  to  as  many  as  330.  Obviously,  in  the  latter  case,  severe  computational  loading  can  result 
and  “intelligent"  algorithms  for  set  selection  are  well  worth  investigation  (particularly  for 
simulations  conducted  an  serial  processors). 

5.9  SV  Initial  Orbital  Parameters 

Initial  conditions  (orbital  parameters)  for  the  optimal  SV  constellation  are  extracted 
from  [12]  and  are  shown  in  Table  5.1.  There  are  six  orbital  planes  in  the  constellation; 
each  plane  contains  four  SVs.  As  previously  noted,  orbital  semi-major  axes  are  26,609 
kilometers  for  all  SVs  and  the  orbital  period  Is  1  I  hours,  69  minutes,  and  57  seconds.  An 
added  benefit  to  this  increased  orbital  radius  is  that-  SVs  are  expected  to  maintain  tholr 
nominal  trajectories  more  successfully,  resulting  in  a  lower  requirement  for  station-keeping 
manuvers  and  a  higher  percentage  of  “in  commission"  time. 
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Table  6.1.  GPS  Optimal  Constellation  Initial  Conditions  [12] 


Satellite  ID 

KTTWnTTfiil 

LAN 

Alpha 

01 

325.730284 

190.96 

13 

145.730284 

312.30 

02 

325.730284 

220.48 

14 

145.730284 

340.93 

03 

325.730284 

330.17 

15 

145.730284 

87.06 

04 

325.730284 

83.58 

16 

145.730284 

209.81 

05 

25.730284 

249.90 

17 

205.730284 

11.90 

06 

25.730284 

352.12 

18 

205.730284 

110.76 

07 

25.730284 

25.25 

19 

205.730284 

143.88 

08 

25.730284 

124.10 

. 

20 

205.730284 

246.11 

09 

85.730284 

286.20 

21 

265.730284 

52.42 

10 

85.730284 

48.94 

22 

265.730284 

165.83 

11 

85.730284 

155.08 

23 

265.730284 

275.62 

12 

85.730284 

-  . - . 

183.71 

24 

265.730284 

305.04 

5, 10  Summary 

This  chapter  presents  the  basic,  concepts  related  to  the  GPS  subsystem  and  introduces 
the  GPS  dynamics  error  model.  In  addition,  the  GPS  pseudo-range  measurement  model 
equations  are  developed  in  detail.  The  conventional  difference  measurement  approach  is 
used  once  again.  Finally,  S  V  orbital  calculations  are  presented  along  with  initial  conditions, 
and  the  criteria  for  set  selection.  The  approach  is  a  considerable  departure  from  previous 
research  [11,  19,  32,  33]  due  to  the  addition  of  a  full  GPS  constellation  in  which  orbits  are 
modeled  explicitly,  and  because  of  the  addition  of  crucial  functions  (PDOP  calculations, 
set  selection,  and  set  switching  algorithms)  which  emulate  the  operation  of  a  GPS  receiver. 
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VI.  Results 


Prior  to  a  moderately  detailed  analysis  of  the  results  obtained  for  the  considerable 
number  configurations  tested,  several  items  of  general  interest  which  affect  most  or  all  of 
the  simulations  are  discussed.  In  addition,  some  of  the  major  differences  between  this  and 
previous  research  are  highlighted. 

6.1  Monte  Carlo  Analysis  of  the  NRS  Error  Model 

The  research  conducted  in  this  thesis  hinges  on  the  development  and  employment 
of  a  system  ( truth  model)  which  is  mainly  constructed  from  models  for  the  LN-93  INS, 
the  RRS,  and  the  GPS  subsystems.  The  fact  that  this  truth  model  is  constructed  in 
software  (embedded  In  MSOFE  [5])  represents  a  significant  extension  to  previous  research 
[31,  33]  in  which  the  truth  models  were  not  installed  in  the  MSOFE  “system"  to  create 
an  environment  in  which  various  Kalman  filter  designs  could  be  tested.  (Rather,  previous 
truth  models  were  tested  in  the  MSOFE  covariance  mode.  Previous  researchers  used  actual 
data  obtained  from  empirical  sources  related  to  CIRIS  operation  to  drive  the  filter  designs 
contained  in  [31,  33].)  The  truth  model  developed  hi  this  research,  along  with  trajectory 
data  generated  by  PROFGEN  [1]  and  the  new  SV  orbit  calculation  software,  generates 
measurement  data  as  well  as  reference  variables  which  are  used  to  test  the  performance 
of  the  (full-order)  Kalman  filter,  and  constitutes  an  environment  in  which  a  variety  of 
full-order  and  reduced-order  Kalman  filters  may  be  tested  against  a  common,  high-fidelity 
standard. 

Several  system-level  configurations  (presented  below)  are  tested  and  smalyzed  in  this 
research.  Generally,  analysis  of  the  configuration  of  interest  consists  of  performing  a  series 
of  10  alignment  runs,  followed  by  a  series  of  10  Hight  simulation  runs,  and  observing 
the  stochastic  time  history  of  error-state  variables  of  interest.  (In  the  case  of  very  large 
dimension  models  with  many  measurement  updates,  fewer  runs  may  be  used,  resulting 
in  lower  confidence  in  the  data  sample  statistics  thus  derived  [22],  Specific  instances  are 
clearly  identified  where  they  exist.) 


8-1 


At  the  beginning  of  each  alignment  run,  the  truth  model  itate  vector  is  initialized  in 
Monte  Carlo  fashions  each  truth  state  Is  set  to  a  random  value  based  on  the  state’s  initial 
covariance  (and  a  random  number  from  a  pseudo-random  number  generator).  At  the  end 
of  the  alignment  runs,  terminal  conditions  for  the  truth  and  Alter  state  vectors  and  the 
Anal  covariances  are  written  to  a  data  Ale.  This  data  Ale  is  then  used  to  provide  initial 
conditions  for  flight  runs.  Thus,  a  true  Monte  Carlo  fashion  simulation  is  preserved  for  the 
entire  analysis  sequence. 

6.2  Feedback  and  Platform  Alignment 

As  noted  above,  the  truth  model  state  vector  elements  are  initially  randomized  in  a 
manner  intended  to  represent  actual  error  conditions  which  may  be  present  in  INS,  RRS, 
and  GPS  subsystems.  The  Intention  is  to  determine  the  effect  of  the  randomisation  of  truth 
model  states  on  the  Kalman  Alter’s  ability  to  perform  its  estimation  task.  This  (stochastic) 
stress  test  for  the  Kalman  Alter  is  another  major  difference  in  this  and  previous  research. 

As  a  result  of  this  stochastic  stress  test,  an  interesting  (but  crucial)  discovery  was 
made  concerning  the  truth  model.  Recall  that  the  INS  model  consists  of  error  states.  For 
instance,  state  number  four  is  the  error  in  north  platform  tilt.  Having  initialized  the  truth 
model  error  states  in  the  Monte  Carlo  fashion  described  above,  the  system  is  allowed  to 
propagate  its  states  for  the  period  of  the  eight-minute  alignment.  During  the  alignment,  the 
Kalman  Alter  is  provided  with  measurements  to  improve  the  estimation  process.  However, 
by  the  end  of  the  alignment  period,  the  truth  model  has  typically  developed  state  variables 
of  an  inappropriately  large  magnitude.  For  example,  the  latitude  and  longitude  error  states 
grow  to  the  order  of  20,000  feet  during  an  eight-minute  alignment. 

The  natural  question  to  raise  here  is  whether  or  not  such  behavior  should  be  expected. 
In  reply,  this  sort  of  behavior  might  easily  be  observed  in  a  physical  INS  that  is  turned  on 
but  NOT  torqued  to  local  level.  Typically,  errors  of  this  magnitude  are  undesirable  after 
having  completed  an  alignment.  Even  though  the  Kalman  Alter  can,  with  tome  forniB  of 
measurement  updating,  maintain  a  reasonably  good  estimate  of  the  misalignment  error 
states,  the  INS  would  subsequently  begin  a  Right  run  (navigation  mode)  with  large  errors 
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rather  than  with  minimal  errors  that  are  seen  in  physical  systems  that  are  “leveled”  during 
the  alignment  process. 

One  might  be  inclined  to  consider  that  the  error  magnitude  is  immaterial  because 
the  interest  is  simply  to  maintain  a  good  estimate  no  matter  what  the  error  magnitude 
happens  to  be.  Two  counterpoints  can  be  made  to  that  argument.  First,  serious  numerical 
complications  invariably  arise  when  attempting  to  propagate  tightly  coupled  error  models 
in  which  state  values  are  becoming  increasingly  large.  In  such  a  case,  numerical  precision  is 
crucial  but  often  impossible  because  computer  wordlength  becomes  a  serious  limiting  factor 
[22].  In  addition,  the  linear  perturbation  model  adequacy  is  placed  in  serious  jeopardy  as 
error  magnitudes  become  excessively  large  [26]. 

How  can  the  issue  be  resolved?  Two  approaches  are  used  successfully  in  this  research. 
First,  one  might  choose  to  use  feedback  during  the  alignment  process  [22].  In  this  approach, 
the  Kalman  filter  state  estimates  are  fed  back  to  the  truth  model.  This  technique  has  the 
effect  of  “leveling"  the  platform  by  reducing  the  magnitudes  of  the  truth  model  error  states. 
The  approach  does  not  appear  to  work  as  well  using  full  state  feedback  as  it  does  when 
using  only  partial  state  feedback  (l.e.,  only  position,  velocity,  and  tilt  states  are  fed  back). 
This  behavior  may  stem  from  the  fact  that  there  exists  a  problem  with  observability  in 
the  INS  error  model.  For  instance,  velocity  measurements  affect  only  the  “basic"  nine 
error  states  (position,  velocity,  and  tilts),  and  provide  no  new  information  to  improve 
the  Kalman  Alter  estimates  of  the  remaining  eighty-four  states.  Additionally,  one  must 
consider  that  feedback  of  some  states  is  impractical  due  to  the  inability  to  correct  errors 
which  may  be  estimated  correctly  but  which  are  physically  inaccessable.  As  an  example, 
the  atmospheric  error  associated  with  RRS  transponders  can  be  estimated  very  well,  but 
one  cannot  change  the  atmospheric  properties  which  create  the  error;  feedback  of  states  in 
this  category  is  not  only  impractical,  it  can  actually  exacerbate  estimation  errors. 

The  second  approach  is  to  use  an  “impulsive  reset”  [22 1  of  the  truth  model  error 
states.  The  first  requirement  is  for  the  Kalman  filter  to  acquire  high  quality  estimates  of 
error  variables.  Subsequently,  instantaneous  corrections  are  applied  to  the  truth  model 
states  (meaning  that  data  registers  are  reset,  based  on  the  available  feedback  states).  It 
is  similar  to  the  feedback  approach  discussed  above,  and  differs  primarily  in  the  frequency 
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at  which  feedback  ii  applied. 

The  “impulsive  reiet”  method  it  adopted  for  general  ute  in  thii  itudy.  Inveitlgatlom 
reveal  that  the  two  method!  provide  comparable  performance  in  reducing  the  truth  model 
error  magnitudes,  but  alio  (not  surprisingly)  revealed  a  noticeable  improvement  in  soft¬ 
ware  execution  time  when  the  second  approach  is  taken.  Consequently,  a  tingle  impulsive 
reset  is  applied  at  the  end  of  each  alignment  run  in  order  to  initialize  flight  runs  with  an 
appropriately  “level”  platform.  Once  again,  only  the  basic  nine  error  states  are  fed  back. 

6.3  Performance  Analyst*  Baseline 

A  complete  error-state  model  for  the  Navigation  Reference  System  (NRS)  Is  devel¬ 
oped  in  Chapters  III  through  V.  The  foundation  is  laid  with  the  93-state  LN-93  RLG-baied 
INS  model.  The  basic  INS  model  contains  a  single  baro-altlmeter  state  which  is  subse¬ 
quently  revised  to  improve  model  fidelity.  Next,  the  RRS  transponder  model  is  added  to 
the  INS  model.  The  NRS  model  is  completed  with  the  addition  of  a  OPS  model  that 
Incorporates  a  full  constellation  of  24  moving  space  vehicles. 

The  results  achieved  from  simulation  and  analyses  at  each  of  the  junctures  described 
above  are  now  presented.  Key  variables  (such  as  position,  velocity,  and  platform  tilt  errors) 
are  plotted  for  each  of  the  test  configurations.  The  plots  are  contained  in  Appendices  C 
through  I,  Salient  features  are  discussed  In  the  sections  which  follow. 

It  is  important  to  begin  research  of  this  type  with  a  “baseline"  to  which  subsequent 
performance  changes  (resulting  from  changes  to  the  model)  are  compared.  Consequently, 
the  first  configuration  which  it  tested  and  analyzed  is  that  of  the  basic  LN-93  INS.  The 
Litton  error  model  contains  93  states,  including  a  single  baro-altimeter  state. 

The  complete  93-stat.e  error  model  Is  programmed  Into  (lie  MSOFE  truth  model  using 
initial  conditions  suggested  by  Litton  [20j,  and  n  full-order  Kalman  filter  is  constructed 
based  oil  the  same  model.  The  model  is  then  tested  In  two  distinct  simulations.  First,  an 
“alignment"  series  (10  Monte  Carlo  runs)  is  performed.  Initial  INS  position  is  assumed  to 
be  45°  North  latitude  and  0°  longitude.  [These  conditions  are  chosen  for  consistency  with 
Litton  conditions.] 
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The  second  test  series  is  a  10- run  Monte  Carlo  simulation  of  the  expected  INS  per¬ 
formance  under  flight  conditions.  Variables  from  the  fighter  flight  profile  presented  in 
Chapter  III  are  used  as  the  nominal  conditions  about  which  the  truth  model  equations 
are  relinearised  during  the  propagation  cycle.  The  extended  Kalman  filter  is  relinearized 
about  its  best  estimates  of  the  trajectory  variables. 

This  phase  of  the  model  analysis  is  conducted  to  ensure  the  baseline  model 
performs  In  a  manner  which  Is  consistent  with  the  performance  stipulated  by  Litton  for 
the  LN-93.  During  the  alignment  series,  both  velocity  and  baro  altitude  measurements 
are  provided  to  the  Kalman  filter.  For  the  flight  series,  only  baro  altitude  measurements 
are  used. 

hi  referring  to  Appendix  C.l,  the  reader  will  find  several  plots  depicting  the  error 
behavior  of  the  93-state  INS  during  alignment.  As  noted  previously,  the  LN-93  error  states 
are  "randomised"  at  the  start  of  each  alignment  -  thus  10  different  random  initial  state 
vectors  are  used  to  seed  the  10  Monte  Carlo  alignment  runs.  The  center  plot  trace  (■■••) 
represents  the  mean  error  time  history  for  the  indicated  state.  Mathematical  descriptions 
for  data  time  histories  are  Included  at  the  beginning  of  Appendix  C. 

The  true  staudard  deviations  of  the  indicated  mean  error  variable  is  represented  by 
<Ti, The  tracei  (•  •  <)  which  bound  and  “track”  the  mean  error  time  history  represent 

the  mean  error  plus  and  minus  crfn„..  The  final  pah  of  traces  ( - )  represents  the  filter- 

computed  for  the  error  variables  Indicated.  They  are  symmetrical  about  zero 

because  the  Kalman  filter  “assumes"  that  Its  errors  are  zero-mean  [22],  The  filter-computed 
error  standard  deviation  magnitudes  (07,71,,,)  may  be  compared  to  similar  plots  in  the 
Litton  reference  [20].  In  the  case  of  the  position  states  (latitude,  longitude,  and  altitude) 
the  comparison  is  excellent!  the  magnitudes  of  the  07,71, plots  for  this  research  compare 
very  closely  to  the  Litton  results  [20]. 

The  only  significant  difference  is  In  the  vertical  velocity  error  state.  The  magnitude 
oflts  0  is  somewhat  lower  here  than  In  the  Litton  document  [20].  This  occurs  due  to  the 
fact  that  Litton  uses  only  horizontal  velocity  updates  during  alignment.  In  this  research, 
vertical  velocity  and  baro-altitude  measurements  are  used  in  conjunction  with  horizontal 
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velocity  measurements. 

Initially,  vertical  velocity  measurement*  are  adopted  because  of  very  large  magnitude 
errors  (discussed  previously)  that  developed  during  alignments.  It  is  postulated  that  the 
random  Initial  condition  effects  could  be  reduced  or  eliminated  by  providing  additional 
measurement  information.  Although  the  referenced  errors  are  not  adequately  quelled  by 
these  additional  measurements,  a  significant  Improvement  in  the  vertical  channel  states 
was  noted  (i.e.,  vertical  velocity  a  decreased  by  about  50  percent  compared  to  the  Lit¬ 
ton  results).  Consequently,  vertical  velocity  and  baro-altltude  measurements  are  used  in 
alignment  runs  throughout  this  research. 

Two  noteworthy  features  appear  In  the  93-state  alignment  plots  contained  in  Ap¬ 
pendix  C.l.  First,  the  latitude  and  longitude  error  time  histories  tend  to  drift  away  from 
the  expected  zero-mean.  This  result  occurs  because  of  the  lack  of  horizontal  position  In¬ 
formation  during  these  alignments  (notice  that  the  mean  altitude  error  looks  “reasonably” 
zero-mean  -  a  consequence  of  Incorporating  baro-altlmetcr  measurements  which  is  another 
departure  from  the  Litton  alignment  procedure). 

The  second  feature  which  is  likely  to  catch  the  reader’s  attention  is  the  flat-line  (zero) 
response  of  A5|.  This  is  the  characteristic  behavior  of  this  vertical  channel  aiding  state 
during  INS  alignment.  This  response  occurs  because  AS\  Is  a  function  of  altitude  rate 
(which  Is  nominally  zero  during  alignment).  (The  behavior  of  AS|  Is  far  more  Interesting 
In  the  flight  runs.) 

In  Appendix  C.2,  the  reader  will  And  the  mean  error  and  standard  deviation  time 
histories  for  several  LN-93  variables  which  result  during  a  lighter  flight  simulation.  Once 
again,  a  10-run  Monte  Carlo  series  is  run  for  the  flight  simulation.  Initial  conditions  for 
the  (light  run*  are  those  final  conditions  which  resulted  from  the  alignment  runs  described 
above. 

The  reader  may  note  that  the  horizontal  and  platform  tilt  error  states  (still)  do  not 
look  zero-mean.  This  Is  again  due  to  the  lack  of  horizontal  position  information.  [For  this 
series  of  runs,  baro-altltude  measurements  are  used  to  bound  vertical  channel  errors.  No 
other  measurements  are  available  to  Improve  the  Kalman  filter  state  estimates.] 
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Ill  Addition  to  the  Urge  magnitude  mean  errors,  the  filter-computed  error  standard 
deviation*  for  many  itate*  (vertical  channel  excepted)  are  roughly  twice  a*  large  ai  the  same 
quantltie*  In  the  Litton  document.  Note  that  this  statement  does  not  imply  that  the  truth 
model  behavior  in  aberrant;  it  means  simply  that  the  93-state  Kalman  Alter  baaed  on  the 
truth  model  is  not  performing  well  as  a  state  estimator.  This  is  not  unexpected;  without 
adequate  measurement  information,  and  considering  the  previously  noted  observability 
issues  Inherent  to  this  problem,  a  Kalman  Alter  cannot  be  expected  to  perform  well  under 
such  liighly  dynamic  conditions  as  those  simulated  in  the  fighter  flight  profile.  However, 
performance  can  be  Improved  as  Indicated  in  subsequent  simulations. 

Although  not  contained  in  the  appendix,  the  93-stat«  truth  model  covariances  were 
plotted  and  compared  to  Litton  results  [20].  The  excellent  agreement  between  those  results 
and  Litton  predictions  constitutes  additional  validation  that  the  truth  model  performs  in 
a  manner  which  is  consistent  with  expected  INS  performance.  The  reader  should  note 
that  this  performance  Is  only  achieved  by  using  feedback  during  the  alignment  -  results 
are  much  different,  without  feedback. 

As  promised,  the  behavior  of  state  13,  AS.{,  is  quite  interesting  in  the  flight  runs. 
The  mean  error,  true  standard  deviation,  and  filter-computed  standard  deviation  ull  have 
the  appearance  of  switching  on  and  off.  This  is  in  fact  the  case.  Recall  that.  A5 1  is  a  strong 
function  of  altitude  rate.  It  Is  “on”  during  intervals  of  altitude  change,  and  “off”  during 
flight  segments  at  constant  altitude,  The  state  values,  and  consequently  the  standard 
deviation,  become  zero, 

6,4  96-State  Error  Model  Performance 

After  adding  the  additional  baro-altitude  states  discussed  in  Chapter  III,  the  same 
sequence  of  testing  (alignments  followed  by  flight  runs)  is  applied.  The  goal  in  this  case 
is  merely  to  assure  that,  by  enhancing  fidelity  through  the  addition  of  states  in  the  truth 
model  and  Kalman  filter,  no  performance  degradation  is  induced  in  the  vertical  chan¬ 
nel  states.  This  step  is  viewed  as  a  “quality  control  check”  prior  to  making  subsequent 
additions  of  subsystem  module  (i.c.  RRS  and  (IPS).  Comparing  the  alignment  plots  ob¬ 
tained  here  to  those  in  Appendix  C.l  (93-state  alignment),  the  reader  is  apt  to  conclude 
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(Incorrectly)  that  the  plot*  are  Identical.  Closer  inspection  reveals  the  differences,  Note 
particularly  the  azimuth  error  state  (Figures  C-3  and  D.3)  and  the  scale  differences  on 
the  vertical  velocity  state  (Figures  C.4  and  D.4),  and  the  altitude  states  (Figures  C.5  and 
D.4). 

The  slight  differences  In  the  vertical  states  are  a  direct  reflection  of  the  change  in  the 
baro-model,  while  the  azimuth  error  anomaly  is  simply  due  to  a  difference  in  the  initial 
condition  for  the  state  which  had  no  effect  on  steady  state  operation  during  alignment 
runs.  The  additional  baro  states  have  no  apparent  effect  on  horizontal  position  (latitude, 
longitude)  errors  during  the  alignment  runs. 

However,  the  differences  are  slightly  more  obvious  in  plots  (tom  flight  runs  for  the  96- 
state  model.  In  this  case  the  filter- computed  standard  deviations  match  Litton  predictions 
quite  closely.  This  Is  primarily  due  to  using  slightly  increased  noise  magnitudes  for  the 
Kalman  filter  (as  opposed  to  the  noises  strengths  set  forth  in  the  Litton  document  which 
are  stiff  used  in  the  truth  model).  In  this  research,  it  1b  determined  that  an  experimentally 
tuved  Kalman  filter  appears  to  achieve  better  altitude  channel  performance  when  a  factor 
of  1.2  to  1.6  Is  applied  to  the  nominal  noise  strength  suggested  by  Litton. 

An  obvious  difference  is  apparent  in  the  vertical  channel  states  (compare  Figure  CM  1 
to  D.ll).  The  large  changes  In  the  true  baro-altimeter  orror,  depicted  in  Figure  D, 11(b), 
drive  similar  changes  in  the  INS  altitude  error,  shown  hi  Figure  D. 11(a). 

The  baro-altimeter  changes  stem  directly  from  the  addition  of  the  new  alimeter  error 
states.  The  most  significant  source  of  the  true  error  is  the  scale  factor  error  state  whicli  is 
included  in  the  truth  model  but  Is  not  included  in  the  filter.  Early  Investigations  revealed 
that  the  addition  of  this  state  to  the  filter  model  resulted  in  severe  numerical  difficulties 
in  the  filter  covariance  progatlon  carried  out  by  MSOFE. 

In  many  instances,  particularly  during  rapid  altitude  transitions,  either  integration 
errors  or  negative  variances  (albiet  extremely  small  magnitude)  resulted.  This  distressing 
result  (negative  variances  arc  mathematical  anomalies  [22])  is  directly  traceable  to  the 
vertical  channel  aiding  state,  AS|  which  varies  strongly  with  altitude  and  altitude  rate. 
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The  remit  is  that  A5|  behaves  quite  erratically  during  rapid  altitude  change  or  at  liigh 
altitude. 


Three  options  are  open  to  resolve  the  problem.  First,  greater  numerical  precision  [he. 
double  precision  variables]  would  likely  reduce  or  eliminate  the  tendency  to  obtain  negative 
variances  [25j.  Unfortunately,  increasing  numerical  precision  also  results  in  significantly 
increased  computational  burden  [22].  Second,  software  may  be  written  to  detect  negative 
covariances  and  (cautiously)  reset  them  to  zero.  This  approach  is  also  undesirable  because 
it  introduces  the  possibility  of  masking  other  problems  which  may  be  present  in  the  Kalman 
filter  [25].  Therefore,  the  third  option  is  taken.  The  problematic  scale  factor  state  is 
removed  from  the  Kalman  filter  while  leaving  the  truth  model  unchanged. 

The  result  is  plainly  evident.  Baro-altimeter  error  in  the  truth  model  depends  on 
scale  factor  error  while  the  Kalman  filter  omits  this  error  source.  Consequently,  during 
high  altitude  portions  of  the  flight  runs,  large  estimation  errors  occur  In  the  baro-altimeter 
channel,  with  attendant  estimation  error  in  the  INS  altitude  channel.  Fortunately,  the 
estimation  error  is  not  unacceptably  "bad"  and  does  not  seriously  degrade  the  overall 
performance  of  the  Kalman  filter,  as  evidenced  by  the  consistency  of  the  horizontal  position 
and  velocity  states. 

With  a  revised  baro-altimeter  model  now  integrated  Into  the  truth  model  (although 
not  in  the  Kalman  filter),  the  next  step  is  to  add  RRS  transponder  aiding.  In  preparation 
for  that  step,  the  96-state  model  is  re-run,  This  time,  however,  different  Initial  coordinates 
are  used  for  alignment  and  flight  run  initializations.  Previously,  Initial  latitude,  longitude, 
and  altitude  were  set  to  45°  north,  0°,  and  0  feet  respectively.  For  this  set  of  runs,  the  initial 
conditions  are  those  of  Holloman  AFB,  (Latitude  is  north  32.78°,  longitude  is  105.99°  West, 
and  altitude  is  4100  feet.)  This  change  is  necessary  due  to  the  fact  that  the  actual  locations 
of  the  RRS  transponders  used  In  subsequent  con  fi  mi  rut  ions  are  on  or  near  Holloman  APB. 

The  differences  between  mean  error  time  histories  from  this  (refer  to  Appendix  E) 
and  the  previous  case  (Appendix  D)  are  generally  unremarkable.  Performance  is  very  con¬ 
sistent  with  that  achieved  previously,  Once  again  the  filter-computed  mean  error  standard 
deviations  are  also  consistent  with  Litton  predictions  [20].  Just  as  in  the  previous  case, 
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Figure  0.1.  96-State  Model,  2-Hour  Flight:  Latitude  ( — ),  and  Longitude  (•••)  Filter- 
Computed  <r  Plots 

latitude  and  longitude  mean  errors  are  not  “white",  duo  to  lack  of  horizontal  position 
information. 

In  order  to  facilitate  a  direct  comparison  of  the  horizontal  positions  states  between 
this  and  subsequent  configurations,  Figure  6.1  is  included.  In  Figure  6.1,  latitude  error 
is  indicated  by  the  solid  trace  and  longitude  error  is  depicted  by  the  dashed  trace.  The 
behavior  exhibited  in  these  plots  is  quite  comparable  to  latitude  and  longitude  error  plots 
in  the  Litton  LN-93  documentation  [20]. 

6.5  72-State  Error- Model 

Occasionally,  a  step  forward  must  be  preceded  by  a  step  (or  two)  backward.  Initially, 
the  RRS  transponder  model  (see  Chapter  IV)  was  successfully  added  to  the  Ofi-state  INS 
model  to  form  a  122-state  truth  model.  Following  that  step,  the  GPS  model  (Chapter 
V)  was  added  to  the  INS/RRS  model  to  form  a  162  state  truth  model.  A  quick  dis¬ 
covery  was  made  at  that  point.  The  accelerator  board  (referenced  In  Chapter  I)  cannot 
handle  a  problem  of  this  size.  Recall  that  with  152  stales,  the  U  -  D  covariance  matrix 
has  152(153)/2  =  11,628  nomedundant  terms  which  must  be  continuously  calculated. 
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VVIk*I  her  because  of  data  storage  needs  or  executable  code  size  requirements  for  random 
access  memory  (RAM),  the  accelerator  board  1b  unable  to  compile  and  link  the  software 
to  create  executable  code, 

This  proble?  i  is  peculiar  to  the  board.  This  point  is  verified  by  successfully  compiling 
and  running  under  a  different  operating  system.  However,  in  the  interest  of  preserving 
the  speed  advantage  gained  by  using  the  accelerator  board,  a  decision  is  made  to  reduce 
the  truth  model  and  Kalman  filter  dimensions  sufficiently  to  circumvent  the  hardware 
limitations  which  arc  encountered. 

References  to  previous  work  [11,  19,  31,  33]  suggested  that  the  INS  states  (originally 
numbered  70  to  93  by  Litton)  are  “expendable.”  A  series  of  Monte  Carlo  alignment  runs 
are  made  to  verify  that  these  states  could  legitimately  be  eliminated  without  significantly 
affecting  the  error  behavior  of  the  truth  model. 

The  plots  contained  in  Appendix  F  directly  compare  the  performance  of  the  reduced 
truth  model  with  the  full  96-state  model.  In  all  cases  the  reduced-order  truth  model 
appears  to  have  “captured  the  essence”  of  the  full-order  model  behavior.  Minor  variations 
are  evident  in  the  altitude  and  vertical  velocity  states.  However,  the  difference  is  mainly 
an  artifact  of  the  scale.  Virtually  all  of  the  difference  may  be  attributed  to  the  simple  fact 
that  these  results  come  from  10-run  simulations.  If  more  runs  are  used  in  the  Monte  Carlo 
simulations  (possibly  26  to  50  [6,  22]  runs),  the  differences  may  diminish  significantly. 
With  just  10  runs,  the  agreement  between  the  two  candidate  truth  models  is  generally 
excellent. 

As  an  aside,  the  96-state  truth  model  covariances  which  are  plotted  in  Appendix  F 
may  be  compared  to  Litton  results  [20],  The  excellent  agreement  between  these  results 
and  Litton  predictions  constitutes  additional  validation  that  the  truth  model  is  performing 
in  a  manner  which  is  consistent  with  expected  INS  performance.  The  reader  should  note 
that  this  performance  is  only  achieved  by  using  feedback  during  the  alignment  -  results 
are  much  different  without  feedback. 

The  obvious  extrapolation  (given  the  excellent  agreement  between  the  96-state  and 
72-state  models)  is  that  the  72-state  model  is  performing  its  INS  emulation  function  with 
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adequate'  fidelity.  Therefore,  the  72-state  model  iu  adopted  to  represent  the  INS  error  ulnte 
subsystem  which  is  integrated  into  the  overall  NR.S  truth  model.  The  additional  states 
(above  the  69  Litton  states  which  are  retained)  are  the  new  baro-altimeter  states. 

6. 6  98-State  Error-Model 

Combining  the  26-state  RRS  model  with  the  72-state  INS  model  produces  the  98- 
state  model  which  is  tested  to  determine  the  baseline  INS/RRS  performance.  Appendix 
G  contains  the  mean,  error,  mean  error  ±<rir,„.t  and  0  ±  cr jm<r  time  history  plots  similar 
to  those  shown  in  previous  configurations.  In  addition  to  the  variables  that  are  plotted 
In  previous  configurations,  results  that  are  representative  of  RRS  transponder  error  Btates 
are  also  plotted  in  Appendix  G. 

The  impact  of  providing  RRS  measurements  to  the  Kalman  filter  is  dramatic  indeed, 
An  obvious  improvement  for  the  alignment  simulations  is  evident  in  the  horizontal  position 
channel  estimates  (refer  to  Figure  G.9).  Latitude  and  longitude  mean  errors  are  now  much 
closer  to  the  zero-mean  predicted  by  Kalman  filter  theory.  In  other  cases,  where  the  mean 
error  is  not  zero-mean  (l.e.,  the  azimuth  error  state,  Figure  G.3),  the  result  is  very  likely  to 
be  due  to  the  fact  that  these  data  are  obtained  from  10-run  Monte  Carlo  simulations.  For 
Improved  statistical  validity,  25  to  50  or  more  Monte  Carlo  runs  are  generally  considered 
necessary  [5,  22].  Careful  comparison  of  the  vertical  scales  for  horizontal  channel  states 
shows  that  the  addition  of  RRS  transponder  measurents  during  alignment  reduces  the 
true-  and  filter-computed  one  sigma  values  by  more  than  50  percent.  Additionally,  there 
are  no  significant  biases  or  ramps  evident  in  the  plots. 

A  far  more  substantial  effect  is  evident  in  the  flight  runs.  With  RRS  transponder 
measurements,  horizontal  position  errors  are  reduced  from  several  thousand  feet  (in  pre¬ 
vious  configurations)  to  less  than  40  feet  during  all  phases  of  the  flight.  The  significant 
variations  in  the  true  and  filter  one  sigma  values  are  due  to  two  sources.  First,  significant 
transitions  occur  in  the  altitude,  velocity,  and  till  states  during  periods  of  high  dynamic 
manueverlng.  Second,  the  aircraft  range  from  the  transponders  has  an  obvious  effect  in 
horizontal  position  errors.  Flight  regimes  in  which  the  aircraft  is  at  low  altitude  or  is  a  long 
distance  from  the  transponders  result  in  increasing  positional  uncertainty.  When  the  air- 
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craft  in  high  overhead  with  respect  to  the  transponder  positions,  much  better  estimations 
are  possible. 

Because  problems  persist  in  the  vertical  channel,  and  with  the  publication  of  a  paper 
by  Litton  engineers  on  the  topic  [2],  the  scale  factor  state  is  (in  this  and  subsequent 
configurations)  essentially  eliminated  from  both  the  truth  and  filter  models.  Although 
obviously  an  undesirable  approach,  the  scale  factor  slate  was  simply  initialized  to  zero  for 
all  runs.  This  change  accounts  for  the  improved  appearance  of  the  altitude  error  estimation 
in  the  later  configurations  (note  the  improvement  in  the  altitude  channels  in  Appendices 
G.2,  H.2,  and  1.2  as  compared  to  Appendices  D.2  and  E.2).  This  step  la  taken  in  the 
interest  of  continuing  research  with  minimal  Interference  from  the  newly  installed  and 
possibly  invalid  [2]  baro  model, 

In  all  cases,  the  filter  appears  to  be  reasonably  well  “tuned”  in  that  the  filter  one 
sigma  values  bound  the  mean  error  traces  the  majority  of  the  time,  without  being 

excessively  conservative,  [Conservative,  as  used  here,  denotes  the  condition  in  which  the 
magnitude  of  the  filter  computed  mean  error  <r  is  larger  than  the  true  <r  magnitude.] 

Ill  order  to  facilitate  a  direct  comparison  of  the  horizontal  positions  states  between 
this,  previous,  and  subsequent  configurations,  Figure  6.2  is  included.  In  Figure  6.2,  latitude 
error  is  indicated  by  the  solid  trace  and  longitude  error  is  depicted  by  the  dashed  trace, 
Note  the  significantly  reduced  <r  magnitudes  in  this  configuration  compared  to  Figure  6.1. 
Performance  is  Improved  by  roughly  three  orders  of  magnitude  when  compared  to  the  INS 
performance  when  aided  by  baro-altitude  only. 

6.7  126-State  INS/RRS/GPS  (NRS)  Model 

In  the  final  phase  of  integration ,  the  CPS  model  is  added  to  the  98-state  TNS/RRS 
mode],  bringing  the  total  number  of  NRS  stales  in  128.  [This  number  was  previously 
determined  to  be  below  the  threshold  at  which  I  lie  accelerator  board  stops  cooperating.] 

Two  modes  of  testing  arc  performed  for  the  CPS  addition.  First  the  alignment  and 
flight  runs  are  performed  with  GPS  measurements;  iu>  RRS  measurements  are  included. 
Second,  the  alignment  and  flight  series  are  re-run  using  both  GPS  and  RRS  measurements. 
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Figure  6.2.  98-State  Model,  2-Hour  Flight!  Latitude  ( — ),  and  Longitude  (••■)  Filter- 
Computed  <t  Plot*  (Note:  RRS  Measurements  Only) 

6.7.1  123-State  NRS  Model:  GPS  Measurements  Only.  Result*  from  the  10-run 
alignment  simulation  of  NRS  aided  by  GPS  (and  the  usual  velocity  and  baro)  measurements 
are  shown  in  Appendix  H.l.  Performance  is  slightly  degraded  compared  to  the  DH-state 
RRS-alded  alignment  (Appendix  G.l).  Even  though  the  values  are  not  substantially 
different  from  the  previous  case,  the  filter  is  estimating  its  own  errors  less  effectively.  The 
difficulty  seems  to  be  limited  to  the  horizontal  position  states.  Other  states  appear  to 
be  equally  well-tuned  as  in  the  previous  case.  The  first  reaction  may  be  to  adjust  the 
driving  noises  associated  with  the  latitude  and  longitude  states.  This  approach  may  work, 
but  extreme  caution  is  warranted.  Another  postulation  is  that  single  precision  numerical 
operations  are  inadequate  for  a  problem  of  this  type.  Such  may  well  be  the  case.  However, 
double  precision  calculations  require  a  significant  increase  in  software  execution  time.  This 
point  is  addressed  in  more  detail  in  Chapter  VIT. 

Two  factor*  must  be  considered.  First,  the  nine  basic  error  states  include  error  angle 
states  in  the  LN-93  model.  The  correlation  between  error  angle  states  and  navigation 
frame  position  errors  does  not  facilitate  the  tuning  process.  (Recall  the  position  error 
transformation  presented  In  Chapter  IV  to  obtain  j  5\,6L,  Sh  ]  from  the  error  angle  vector 
(  S9.n  S$u,  50t  ].)  The  hazard  is  that  "tuning"  one  of  the  error  angle  states  generally  affects 
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nil  of  the  position  stntes  *  sometimes  In  an  undesirable  fashion.  The  extent  of  she  hae-anl 
goes  beyond  that,  however.  Because  of  the  highly  coupled  nature  of  the  basic  nine  error 
states,  adjusting  the  noise  level  on  one  state  may  cause  adverse  effects  on  several  others. 
Caution  and  patience  are  the  only  solutions  to  the  tuning  problem. 

Three  one-hour  flight  runs  are  performed  for  the  128-state  NRS  model  (compared  to 
10  two-hour  runs  in  previous  cases).  The  execution  time  for  a  problem  of  this  slsc  is  a  major 
concern.  As  an  example,  the  10-run  Monte  Carlo  simulation  for  the  alignment  runs  takes 
approximately  15  to  18  hours.  Consider  that  the  alignments  (480  seconds  each)  are  about 
one-fifteenth  as  long  as  the  flight  runs.  Additionally,  flight  trajectory  data  interpolation 
increases  flight  simulation  times  significantly.  As  a  result,  a  10-run  flight  simulation  series 
is  estimated  to  require  a  minimum  of  225  hours  (more  than  nine  days)  to  complete  on  the 
fastest  processor  generally  available  at  AFIT,  [There  is  hope  for  the  future.  This  issue  is 
addressed  again  in  Chapter  VII. j 

The  reduction  In  flight  time  is  really  not  e  serious  drawback.  In  earlier  cases,  the 
two-hour  flight  profile  is  used  simply  to  facilitate  direct  comparisons  of  the  configuration 
under  test  and  the  Litton  reference  [20].  Such  comparisons  may  just  as  easily  be  made  for 
one-hour  flight  times.  The  validity  of  this  point,  is  evident  if  the  reader  observes  the  error 
behavior  exhibited  in  the  RRS  aided  INS  case  (Figure  G.9).  The  maximum  estimation 
ertors  occur  at  approximately  2000  seconds  (during  rapid,  high-g  manueverlng  far  from 
the  transponder  sites),  and  have  returned  to  lower  magnitudes  by  the  end  of  the  first  hour 
of  flight. 

However,  the  reduction  to  three  runs  in  the  Monte  Carlo  series  Is  cause  for  concern 
[22].  In  cases  such  as  this,  confidence  in  the  validity  of  the  statistical  data  is  reduced.  (In 
fact,  10-run  aeries  for  problems  of  this  type  are  considered  to  be  marginally  adequate  in  a 
statistical  sense  [22],)  Therefore,  quantitative  comparisons  of  performance  based  on  three 
runs  should  probably  be  avoided.  On  the  other  hand,  cautious  qualitative  comparisons  arc 
probably  still  reasonable. 

With  that  qualification,  a  comparison  between  tills  and  the  previous  configuration 
(pre-GPS)  is  presented.  The  3-run  flight  results  In  which  GPS  (but  no  RRS)  measurements 
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are  used  to  old  the  INS  ore  shown  In  Appendix  II. 2. 

The  general  characteristic  which  seems  to  be  evident  In  virtually  all  error  states  (at 
least  those  whose  time  histories  are  plotted)  Is  that  the  filter-computed  a  magnitudes  ap¬ 
pear  to  he  consistently  lower  for  the  GPS-aided  INS,  compared  to  the  RRS-aided  INS. 
However,  the  true  er  values  are  somewhat  larger  than  previously  and  the  mean  error  time 
histories  are  not  satisfactorily  bounded  by  the  filter- computed  <r  magnitudes,  This  result 
is  attributed  to  two  factors.  First,  data  plots  are  based  on  3-run  Monte  Carlo  flight  sim¬ 
ulations.  Second,  numerical  problems  which  sometimes  result  from  using  single  precision 
versus  double  precision  are  quite  possibly  at  fault.  (This  issue  is  alio  addressed  further  in 
Chapter  VII.) 

In  order  to  facilitate  a  direct  comparison  of  the  horisontal  positions  states  between 
this,  previous,  and  subsequent  configurations,  Figure  6.3  is  included.  In  Figure  6.3,  latitude 
error  Is  indicated  by  the  solid  trace  and  longitude  error  is  depicted  by  the  dashed  trace, 
Note  the  significantly  reduced  <r  magnitudes  in  this  configuration  compared  to  Figure  6.2. 
In  some  flight  regimes  (notably  those  far  from  the  tranponder  sites,  or  those  with  depressed 
elevation  angles  between  the  user  and  transponders)  performance  appears  to  have  improved 
by  about  one  order  of  magnitude  when  compared  to  the  INS  performance  achieved  with 
llRS  aiding  alone. 

6.8  128-Slate  INS/RRS/GPS  (NRS)  Model 

6.8.1  128-State  NRS  Model:  GPS  and  RRS  Measurement s  For  this  configuration, 
both  RRS  and  GPS  measurements  are  used.  For  the  alignment  simulations,  10  Monte  Carlo 
runs  are  performed.  However,  only  a  single  Monte  Carlo  run  is  performed  for  the  flight 
simulation.  (Therefore,  extreme  caution  must  be  exercised  in  any  analytical  conclusions 
drawn  from  the  flight  results.) 

Alignment  results  look  quite  similar  to  results  obtained  with  either  GPS  or  RRS 
alone.  The  general  trend  is  a  very  slight  reduction  in  the  filter  estimates  of  the  overall 
error  magnitudes  during  the  alignment  phase;  no  remarkable  behavior  is  observed. 

On  the  other  hand,  flight  results  are  not  quite  as  optimistic.  Although  it  is  clear  that 
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Figure  0,3.  128-State  Model,  2-Hour  Flight:  Latitude  ( — ),  and  Longitude  (••  •)  Filter- 
Computed  <r  Plot*  (Note:  GPS  Measurements  Only) 


the  filter  “tlilnks”  that  It  li  doing  a  better  job  of  estimating  the  variables  of  interest,  it  is 
in  fact  not  doing  as  well  as  it  thinks.  The  filter,  aided  by  both  RRS  and  GPS,  “thinks” 
that  the  estimation  errors  are  smaller  than  they  are.  Its  estimates  of  error  variances  art- 
consistent  with  expectations ,  but  are  inappropriate  for  the  situation,  For  all  error  states, 
the  mean  error  time  history  trace  should  be  bounded  (08  percent  of  the  time)  by  the  filter- 
computed  c r,  This  is  not  the  case  for  some  states,  notably  those  of  horizontal  position 
errors. 

The  apparent  difficulty  in  this  case  is  attributed  to  two  factors.  First,  one  must 
return  to  the  issue  of  using  single-precision  versus  double- precision  for  a  problem  of  this 
nature.  Such  a  choice  may  easily  account  for  some  of  the  unexpected  filter  behavior. 
Additionally,  the  fact  that  a  single  flight  run  is  presented  is  cause  for  carefully  considering 
the  data  which  are  thus  obtained.  Many  of  the  "temporary  divergences”  which  occur  at 
various  high  dynamic  points  of  the  flight  regime  are  statistically  “smoothed”  when  a  larger 
number  of  runs  are  used. 

In  order  to  facilitate  a  direct  comparison  of  thy  Kalman  Alter  horizontal  positions 
states  between  tills,  previous,  and  subsequent  configurations,  Figure  6.4  is  included.  In 
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Figure  0.4.  128-State  Model,  2‘Hour  Flight!  Latitude  ( — ),  and  Longitude  (•  •  •)  Filter- 
Computed  9  Plote  (Note:  RRS  and  GFS  Measurement*) 


Figure  0.4,  latitude  error  ie  Indicated  by  the  aolld  trace  and  longitude  error  la  depicted  by 
the  daahed  trace.  The  9  magnitude!  In  thi*  configuration  are  (lightly  reduced  compared  to 
Figure*  6.2  and  6.2.  The  alight  ramping  effect  ia  exacerbated  by  the  plot  scale.  Although 
such  behavior  ia  not  consistent  with  theory  [IT],  It  may  be  explained  by  the  small  magnitude 
pseudo-noises  which  are  used  in  the  Alter  model  for  the  GPS  position  states.  The  noise 
strength  may  be  somewhat  higher  than  necessary  to  preclude  the  filter  covariances  from 
decreasing  to  aero.  Experimentation  with  these  "tuning’'  parameters  is  warranted. 

After  reviewing  the  configuration*  Involving  GPS  measurements,  two  conclusions  are 
drawn.  First,  GPS  aiding  alone  appears  to  offer  a  noticeable  improvement  over  simple 
tranponder-aided  INS.  In  several  flight  regimes,  the  GPS-aided  solution  is  significantly 
better  than  that  of  the  RRS-alded  INS.  This  fact  is  a  direct  result  of  two  features  which 
are  inherent  to  this  simulation  (and  which  may  imt  necessarily  be  true  In  general).  (The 
flight  path  extendi  far  beyond  the  optimal  coverage  areas  for  the  fixed  transponders,  and 
only  six  transponders  are  used;  many  actual  flights  use  twenty  or  more  transponders.) 
Second,  the  true  Kalman  filter  estimation  errors  are  not  adequately  bounded  by  the  filter- 
computed  9  (as  they  are  in  previous  cases)  for  the  simulations  in  which  GPS  measurements 
are  used. 
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The  Apparent  problem  with  the  true  estimation  error*  may  bo  attributed  to  nu¬ 
merical  preclilon  problem*  reaultlng  from  single-precision  computatlona,  The  dealre  to 
avoid  doublc-predalon  simulations  la  prevloualy  dlacuaaed.  However,  It  la  poaalble  that 
thia  reticence  reaulted  In  leaa  than  apectacular  filter  ;,erformance,  A*  In  moat  engineering 
•ltuatlona,  the  “trade-offi"  muat  be  carefully  coniidered  and  the  potential  for  problema 
underatood. 

Another  lmportan*  conclualon  which  is  tentatively  drawn  la  that  the  configuration  In 
which  OPS  and  RRS  meaaurementa  arc  uaed  together  offers  the  beat  overall  performance. 
Although  the  trua  error  behavior  doe*  not  appear  to  be  conalatent  with  filter  predlctlona, 
rlmple  Alter  tuning  01  ai'Wcrease  In  numerical  preclilon  may  be  rufflclert  to  reduce  or 
eliminate  the  dlaparity. 

6,9  Summary 

Thl*  chapter  present!  the  reiult*  from  a  large  number  of  configuration- dependent 
Monte  Carlo  ilmulationi,  Each  configuration  la  teated  in  both  an  alignment  and  flight 
simulation.  General  results  from  each  case  are  diliu.sed,  with  reference  to  particular  plots 
(contained  In  appendices  and  distributed  in  the  text)  for  emphaala  of  important  points. 
The  reader  la  cautioned  of  the  potential  for  observability  and  numerical  problema  which 
may  stem  from  simulations  of  the  type  conducted  here. 
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VII.  Conclusions  and  Recommendations 


Thii  effort  li  concluded  with  aome  thoughti  on  Ita  objective*,  the  degree  of  lucceaa 
achieved,  end  the  potential  to  epawn  future  reaearch.  The  major  goal*  of  thii  research 
are  closely  tied  to  the  ultimate  objective  of  providing  the  Central  Inertial  Guidance  Test 
Facility,  Holloman  AFB,  NM,  with  a  superior  navigation  reference  system  that  can  be  used 
as  the  standard  to  which  all  other  navigation  systems  are  compared.  This  is  a  non-trivlal 
objective  which  demands  the  “optimal"  Integration  of  several  highly  complex  systems  and 
which  encompasses  a  large  body  of  knowledge  from  several  engineering  disciplines.  The 
variety  of  disciplines  and  technologies  needed  to  construct  such  a  reference  system  success¬ 
fully  is  at  once  exhilarating  and  staggering, 

7,1  Bare- Altimeter  Model 

In  the  attempt  to  further  the  broader  goal  stated  above,  the  goal  of  improving  the 
baro-altimeter  portion  of  the  LN-93  error  model  is  undertaken.  Though  Innocuous  at  a 
glance,  tills  portion  of  the  project  turns  out  to  be  far  more  troublesome  than  anticipated, 
confounds  a  great  deal  of  effort  and  consumes  an  inordinate  amount  of  time  which  would 
have  been  better  devoted  to  other  pursuits. 

It  appears  that  this  research  validates  a  result  that  Litton  has  recently  published 
in  n  paper  presented  at  the  National  Technical  Meeting  of  the  Institute  of  Navigation, 
January  1991  [2].  In  essence,  the  paper  states  that  the  vertical  channel  model  Included  in 
the  original  LN-93  document  was  not  fully  evolved,  and  a  revised  model  is  presented. 

In  addition  to  pointing  out  how  the  new  model  differs  from  the  previous  Litton 
model,  the  paper  also  describes  a  scale  factor  model  that  Is  totally  Incompatible  with  the 
scale  factor  state  developed  In  the  revised  bant  model  presented  in  Chapter  III.  Because 
the  vertical  channel  issue  had  already  been  “pul  to  rent"  In  this  research  when  this  fact 
came  to  light,  the  model  was  not  redesigned.  Rather,  the  scale  factor  state  was  simply 
Initialised  to  sero  for  all  runs.  This  change  accounts  for  the  improved  appearance  of  the 
altitude  error  estimation  in  the  later  configurations  (note  the  Improvement  in  the  altitude 
channels  In  Appendices  G.2,  H.2,  and  1,2  as  compared  to  Appendices  D.2  and  E.2).  After 
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“wrolng"  the  scale  factor  state,  the  truth  and  filter  models  are  ngnlit  [effectively]  defined 
identically;  neither  includes  the  scale  factor  error  model.  The  choice  to  set  the  scale  factor 
contribution  to  sero  is  not  correct;  it  is  merely  expedient  to  preclude  interference  with  the 
other  major  objectives  of  this  work.  It  is  fully  recognised  that  a  more  appropriate  model 
is  required. 

Consequently,  revision  of  the  baro  model  is  both  a  success  and  a  failure  in  terms  of 
thesis  objectives.  The  baro  model  was  indeed  revised  and  was  implemented  in  the  truth 
model  and,  to  a  lesser  extent,  In  the  filter  model.  Its  operational  characteristics  Including 
the  tendency  toward  numerical  problems  in  the  Kalman  filter  were  identified.  However, 
new  facts  have  now  surfaced  which  indicate  that  the  proposed  model  was  not  entirely 
appropriate  [2];  the  result  Is  that  the  new  model  is  revised  to  contain  three  states  rather 
than  four. 

The  solution  to  the  difficulties  with  the  vertical  channel  states  Is  (hopefully)  quite 
simple.  Future  researchers  in  this  area  must  be  provided  access  to  a  current  error  model 
for  the  LN-93,  Including  the  newly  revised  vertical  channel  mechanisation. 

7J  GPS  Model 

Another  goal  established  at  the  outset  was  to  devise  software  to  calculate  SV  orbital 
positions  for  the  interval  of  the  simulated  aircraft  flight  run.  This  goal  is  completed.  In 
fact,  software  of  a  general  nature  is  fully  developed  and  tested  to  solve  the  SV  orbital 
mechanics  problem  for  arbitrary  intervals  of  time,  anu  to  find  and  select  the  best  set  of 
SVs  for  performing  GPS  measurements  In  a  moving  vehicle  (user).  These  achievements 
represent  significant  extensions  of  previous  research. 

7,3  Truth  Model 

The  major  goal  of  assembling  and  testing  a  high-fidelity  NRS  truth  model  which  is 
embedded  in  MSOFE  is  also  completed.  The  complete  error  model  is  thoroughly  docu¬ 
mented  and  is  fully  operational  in  software.  This  provides  an  environment  within  which 
future  Kalman  filter  designs  may  be  evaluated.  This  achievement  is  another  major  step 
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beyond  prevloui  research  and  ii  particularly  important  because  it  provides  an  environ¬ 
ment  in  which  reduced-order  Kalman  Alter  designs  may  be  tested  prior  to  operational  use 
at  CIGTF. 

As  a  first  step  in  the  future,  AFIT  research  should  focus  on  achieving  a  signifi¬ 
cant  reduction  in  the  number  of  states  in  the  Kalman  filter  used  to  estimate  NRS  errors. 
This  step  is  not  only  important  for  potential  operational  considerations,  but  will  also  aid 
tremendously  by  “creating  space”  for  other  enhancements  to  the  NRS  model.  (One  such 
enhancement  is  differential  GPS.)  Additionally,  state  reduction  in  the  Kalman  Alter  1b 
expected  to  have  the  advantage  of  reducing  the  observability  problems  which  have  been 
noted  to  exist  in  the  full-order  Kalman  filter.  Observability  analyses  (see  Section  7.5)  will 
assist  in  state  reduction  efforts. 

Inherent  to  the  state  reduction  process  is  the  need  to  establish  the  adequacy  of  the 
reduced  INS  truth  model.  (Recall  that  the  INS  truth  model  is  reduced  from  96  to  72  states 
In  this  research.)  Although  a  10-run  Monte  Carlo  series  is  performed  to  show  that  the  72 
state  model  fidelity  is  adequate,  further  analysis  may  be  warranted.  It  is  recommended  that 
future  researchers  run  the  72  state  Kalman  filter  against  the  96  state  INS  truth  model  and 
examine  the  residuals  for  appropriate  statistical  properties  [26]  to  establish  conclusively  if 
the  chosen  state  reduction  Is  adequate. 

7.4  Monte  Carlo  Analyses 

To  the  maximum  extent  possible,  the  requirement  to  determine  performance  char¬ 
acteristics  (via  Monte  Carlo  analyses  )  for  each  of  the  previously  defined  subsystem  and 
system  level  combinations  has  been  satisfied.  A  series  of  10-run  Monte  Carlo  series  is  per¬ 
formed  for  fen  of  twelve  configurations.  In  the  other  two  cases  (the  largest  dimensions  in 
the  truth  model  and  Kalman  filter),  a  three  run  series  is  performed  for  qualitative  analysis. 
It  seems  fair  to  say  that  the  goal  Is  largely  sutisllvd,  although  it  cannot  be  claimed  that 
the  two  series  which  contained  only  three  runs  are  conclusive  proof  that  NRS  functions  as 
well  as  desired. 


7-3 


7.5  Opportunities  for  Future  Research 

Some  topics  which  are,  necessarily,  given  cursory  treatment  in  this  research  might 
easily  form  the  basis  for  extensive  research  in  the  future.  One  such  area  is  that  of  error 
budgeting  [22]  and  optimal  tuning  of  the  NR£>  Kalman  filter.  As  noted  previously,  tuning 
a  Alter  of  the  sise  represented  here  is  a  formidable  task.  However,  it  is  one  which  should  be 
undertaken.  The  major  requirement  is  access  to  significantly  more  computational  power 
than  that  which  is  generally  available  to  AFIT  students.  Options  include  the  new  Sun 
Sparc- Stations  (operating  at  28  MIPS)  which  are  expected  to  be  “on-line”  in  the  near 
future.  In  addition,  every  effort  should  be  made  to  use  double-precision  variables  and  to 
perform  all  calculations  in  double-precision  as  well.  Many  of  the  numerical  problems  which 
are  evident  in  this  research  may  be  eliminated  by  using  double-precision  variables. 

Another  suggestion  that  goes  hand-in-hand  with  error-budgeting  is  that  of  Kalman 
filter  state  reduction.  As  noted  in  Section  6.5,  a  highly  desirable  performance  characteri¬ 
sation  test  is  to  run  any  proposed  reduced-order  Kalman  filter  against  the  full  order  truth 
model  and  observe  the  statistical  properties  of  the  residuals  [26].  Investigations  should 
focus  on  determining  the  degree  to  which  the  residuals  satisfy  the  properties  of  zero-mean, 
white,  Guasslan  statistics  In  order  to  determine  the  adequacy  of  the  reduced-order  model 
[26], 

After  simulation  efforts  complete  the  initial  design  and  tuning  of  a  Kalman  filter, 
another  form  of  testing  is  typically  applied  to  ensure  that  the  filter  will  perform  well  in 
the  job  for  which  it  was  designed.  Consequently,  the  Kalman  filter  should  be  supplied 
actual  data  which  has  been  collected  from  operational  ARS  or  NRS  hardware.  The  only 
data  of  this  type  which  has  been  collected  to  date  was  obtained  from  CIRIS  using  the 
LN-15  or  LN-39  INS.  Although  the  system  errors  in  the  two  INS  models  are  believed  to 
be  roughly  equivalent,  the  effect  of  "driving"  the  l.,N  -IK)  INS  with  data  which  is  collected 
from  the  LN-39  INS  has  not  been  verified  prior  to  (his  thesis.  It  is  assumed  that  the  effect 
is  negligible  [17].  This  assumption  can  be  verified  or  refuted  in  future  research. 

Another  area  which  seems  intriguing,  highly  valuable,  and  under  developed  is  that 
of  Kalman  filter  state  reduction  techniques.  In  many  cases,  a  filter  construction  based 


7-4 


on  a  truth  model  hat  potential  observability  problems.  Generally)  such  problems  can 
be  identified  via  observability  analyses  [22].  In  fact,  such  analysis  is  attempted  in  this 
research  but  the  only  tool  suitable  for  the  task  (MATRIX, \)  is  not  capable  of  calculating 
the  observability  Gramian  matrix  of  the  dimensions  necessary  for  this  thesis.  It  can  be 
an  extremely  valuable  project  to  develop  robust  software  which  is  specifically  designed  to 
produce  a  high-accuracy  (i.e.  double-precision)  observability  Gramian  for  high  dimensional 
systems  such  as  the  NRS.  This  would  provide  an  extremely  valuable,  time-saving  and  cost¬ 
saving  tool  to  future  Kalman  filter  designers  who  are  faced  with  design  considerations 
which  necessitate  system  dimension  reduction. 
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Appendix  A.  NRS  Error  Model  State  Definitions 

The  complete  NRS  error  model  is  composed  of  96  LN-93  INS  error  states  (including 
the  baro-altimeter  states),  26  RRS  transponder  states,  and  30  GPS  error  states,  for  a  total 
of  152  states  in  the  truth  model.  Naturally,  there  are  152  states  in  the  full-order  filter 
model  as  well.  The  non-zero  elements  of  the  F  and  Q  matrices  of  the  LN  -93  error  model 
are  presented  In  Appendix  B. 

The  reader  should  note  that  several  errors  exist  in  the  Litton  reference  [20].  Many 
such  errors  are  identified  in  [13].  Others  are  noted  in  [19]. 

A.J  INS  and  Baro-altimeter  Error  States 

Tables  A.l  through  A.4  describe  the  LN-93  error  model  (93  states)  as  defined  in 
the  Litton  LN-93  CDRL  [20].  Note  that  three  additional  barometric  altitude  states  which 
were  not  explicitly  included  in  the  Litton  model  are  added  for  a  total  of  96  states.  A 
detailed  discussion  on  the  need  for  and  development  of  the  additional  baro-altimeter  states 
is  included  in  Chapter  3.  These  additional  baro-altimeter  states  are  Inserted  in  the  error 
state  vector  sequence  at  numbers  24,  25,  and  26.  Additionally,  the  original  state  23  lh,-  is 
now  labeled  Shy.  Consequently,  INS  error  states  which  were  originally  numbered  24  -+  93 
are  now  renumbered  as  27  -*  96. 

A. 2  RRS  Transponder  Error  States 

Table  A. 5  defines  the  RRS  transponder  error  states  as  they  are  modeled  in  the  NRS. 
These  states  are  defined  in  and  extracted  from  [31].  A  total  of  26  states  are  included  to 
model  the  error  characteristics  of  six  ground  transponders  plus  interrogator  error  sources, 

A. 3  GPS  Error  States 

Table  A.6  defines  the  GPS  error  states  as  Ihe.v  are  modeled  In  the  NRS.  These 
states  are  defined  in  and  extracted  from  [11].  The  definitions  are  believed  to  be  ultimately 
traceable  to  the  paper  by  D.B.  Cox  [6].  A  total  of  30  states  are  included  to  model  the  error 
characteristics  of  4  space  vehicles  plus  user  equipment  error  sources. 
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Table  A.l.  NES  System  Model:  INS  States  1  -*  22 


State 

Number 


State 

Symbol 


Definition 


1 


69* 


X-component  of  vector  angle  from  true  to  computer  frame 

i  i*  i  4*  i  i  5* 


Y» component  of  vector  angle  from  true  to  computer  frame 
Z-component  of  vector  angle  from  true  to  computer  frame 


ee , 


X-component  of  vector  angle  from  true  to  platform  frame 


Jk, 


Y-component  of  vector  angle  from  true  to  platform  frame 


<t>t 


Z-component  of  vector  angle  from  true  to  platform  frame 


6Vv 


X-component  of  error  In  computed  velocity 


SVu 


Y-component  of  error  in  computed  velocity 


SV, 


Z-component  of  error  In  computed  velocity 


10 


Sh 


Error  in  vehicle  altitude  above  reference  ellipsoid 


11 


ShL 


Error  In  lagged  Inertial  altitude 


12 


6S:i 


Error  in  vertical  channel  aiding  state 


13 


iS.t 


Error  in  vertical  channel  aiding  state 


14 


X-component  of  gyro  correlated  drift  rate 


15 


Y-component  of  gyro  correlated  drift  rate 


16 


Z-component  of  gyro  correlated  drift  rat e 


17 


X-component  of  accelerometer  and 
velocity  quantizer  correlated  noise 


18 


Y-component  of  accelerometer  and 
velocity  quantizer  correlated  noise 


19 


Z-component  of  accelerometer  and 
velocity  quantizer  correlated  noise 


20 


S9j- 


X-component  of  gravity  vector  errors 


21 


22 


Y-component  of  gravity  vector  errors 


6g, 


Z-component  of  gravity  vector  errors 


r 

it 


Table  A.2.  NRS  System  Model;  INS  States  23 


Table  A.3.  NRS  System  Model;  INS  States  51  72 


State 

Number 

Definition 

51 

X-component  of  accelerometer  bias  repeatability 

52 

mjm 

Y-component  of  accelerometer  Irian  repeatability 

63 

Vf„ 

Z-component  of  accelerometer  bias  repeatability 

54 

Sa, 

X* component  of  accelerometer  and  velocity 
quantlier  icale  factor  error 

65 

Sav 

• 

Y-component,  of  accelerometer  and  velocity 
quantiser  icale  factor  error 

66 

s , 

Z-component  of  accelerometer  and  velocity 
quantiser  icale  factor  error 

57 

Sqa, 

X-component  of  accelerometer  and  velocity 
quantiser  icale  factor  asymmetry 

58 

SQAy 

Y-component  of  accelerometer  and  velocity 
quantiser  Scale  factor  asymmetry 

59 

Sqa, 

Z-component  of  accelerometer  and  velocity 
quantiser  icale  factor  aiymmetry 

60 

lll^g 

Coefficient  of  error  proportional  to  iquare 
of  meaiured  acceleration 

61 

fvu 

Coefficient  of  error  proportional  to  iquare 
of  meaiured  acceleration 

62 

fit 

Coefficient  of  error  proportional  to  square 
of  measured  acceleration 

63 

: 

Coefficient  of  error  proportional  to  products  of  acceleration 
along  and  orthogonal  to  accelerometer  sensitive  axis 

64 

Coefficient  of  error  proportional  to  products  of  acceleration 
along  and  orthogonal  to  accelerometer  sensitive  axis 

65 

m i 

Coefficient  of  error  proportional  to  products  of  acceleration 
along  and  orthogonal  to  accelerometer  sensitive  axis 

mi 

Coefficient  of  error  proportional  to  products  of  acceleration 
along  and  orthogonal  to  accelerometer  sensitive  axis 

67 

fm 

Coefficient  of  error  proportional  to  products  of  acceleration 
along  and  orthogonal  to  accelerometer  sensitive  axis 

mu i 

u 

Coefficient  of  error  proportional  to  products  of  acceleration 
along  and  orthogonal  to  accelerometer  sensitive  axis 

69 

_ _ 

X  accelerometer  misalignment  about  Z-oxis 

Y  accelerometer  misalignment  about  Z-axls 

71 

Z  accelerometer  misalignment  about  Y-axis 

72 

1 

Z-accelerometer  misalignment  about  X-axis 

Table  A.4.  NRS  Syetem  Model:  INS  State*  73  -»  90 


State 

Number 

State 

Symbol 

Definition 

m 

m 

X-  component  of  accelerometer  biaa 

thermal  tranalent 

74 

B 

Y-component  of  accelerometer  bias 
thermal  transient 

75 

Z-eomponent  of  accelerometer  bias 
thermal  transient 

70 

X-component  of  initial  gyro  drift  rate 
bias  thermal  transient 

Y-component  of  initial  gyro  drift  rate 
bias  thermal  transient 

78 

Z-component  of  initial  gyro  drift  rate 
bias  thermal  transient 

79 

mmm 

X  gyro  compliance  term 

80 

X  gyro  compliance  term 

81 

F (W 

X  gyro  compliance  term 

82 

X  gyro  compliance  term 

83 

F„, 

X  gyro  compliance  term 

84 

Fr« 

X  gyro  compliance  term 

85 

Y  gyro  compliance  term 

86 

F„„ 

Y  gyro  compliance  term 

87 

Fv>v 

Y  gyro  compliance  term 

IBi 

mam. 

Y  gyro  compliance  term 

89 

B'JIPB 

Y  gyro  compliance  term 

90 

Y  gyro  compliance  term 

91 

kAh 

'■m 

Z  gyro  compliance  term 

92 

Z  gyro  compliance  term 

93 

Ft.,., 

Z  gyro  compliance  term 

94 

HUfflHI 

Z  gyro  compliance  term 

95 

Fiyy 

Z  gyro  compliance  term 

96 

Z  gyro  compliance  term 

l‘5 


Table  A. 5.  NRS  System  Model:  RRS  Error  State* 


State 

Number 

■rmji 

Definition 

[NOTE:  Sins  -  Total  INS  State*]  . 

Sins  +  1 

mt£m 

Range  error  due  to  equipment  bia* 

Sins  +  2 

Svb 

Sins  +  3 

mum 

Traniponder  1  x-component  of  potition  error 

Sins  +  4 

Miasm 

Tramponder  1  y-component  of  poiitlon  error 

Sins  +  5 

miasm 

Transponder  1  ■•component  of  poiitlon  error  \ 

Sins  +  8 

mmm 

TVansponder  1  range  error  due  to  atm  propagation 

Sins  +  7 

SPn. 

Traniponder  2  x-component  of  position  error 

Sins  +  8 

musm 

Transponder  2  y-component  of  position  error 

Sins  +  9 

miasm 

Transponder  2  ■•  component  of  position  error 

Sins  + 10 

mmm 

Transponder  2  range  error  due  to  atm  propagation 

Sins  +  11 

mum 

Sins  + 12 

SPn  •„ 

Sins  + 13 

JPn, 

Transponder  3  s-component  of  position  error 

Sins  +  14 

Mmm 

Transponder  3  range  error  due  to  atm  propagation 

Sins  +  18 

mam 

Transponder  4  x-component  of  position  error 

Sins  + 18 

SPn* 

Transponder  4  y-component  of  position  error 

Sins  +  17 

SPn. 

Transponder  4  s-component  of  position  error 

Sins  + 18 

SP'I'i, 

Traniponder  4  range  error  due  to  atm  propagation 

Sins  +  19 

mam 

Traniponder  8  x-component  of  position  error 

Sins  +  20 

■Z2E3V 

Traniponder  6  y-component  of  position  error 

Sins  +  21 

SPn. 

Traniponder  6  s-component  of  position  error 

Sins  +  22 

SPTin 

Traniponder  6  range  error  due  to  atm  propagation 

Sins  +  23 

SPn. 

Transponder  6  x-component  of  position  error 

Sins  +  24 

SPn* 

Traniponder  8  y-component  of  position  error 

Sins  +  26 

mjmm 

Traniponder  6  s-romponent  of  position  error 

Sins  +  28 

Traniponder  6  range  error  due  to  atm  propagation 

Table  A.6.  NRS  Syitem  Model:  QPS  Error  State* 


State 

Number 


ms  +  S/ihs  +  1 


Sms  +  Srhs  +  2 


ins  +  Suns  +  3 


ms  +  Suns  +  4 


ins  +  Smts  +  5 


ins  +  S  hhs  +  ® 


Sms  +  Suns  +  7 


Sms  +  Smts  +  8 


Sms  +  Smts  +  9 


Sms  +  Sims  +  10 


Sms  +  Smts  +  11 


Sms  +  Suits  +  12 


ms  +  Situs  +  13 


Sms  +  Suns  +  14 


Sms  +  Suits  +  18 


com 


6RultM) 


BRt 


usm 


Sms  +  Suits  +  18 


SR,iw; 


Sms 

Q 

D 

20 

Sms 

+ 

Suns 

+ 

21 

Sms 

+ 

Situs 

+ 

22 

Sins 

+ 

Situs 

+ 

23 

Sins 

+ 

Suits 

+ 

24 

Sms 

+ 

Situs 

+ 

25 

Sins 

+ 

Suits 

+ 

26 

Definition 

NOTE:  Situs  =  Total  RRS  State*] 


User  clock  bla* 


U*er  clock  drift 


SV  1  code  loop  error 


1  tropotpherlc  error 


'  1  lonoipherlc  error 


SV  1  clock  error 


SV  1  x*  component  of  poiitlon  error 


SV  1  y-component  of  poiitlon  error 


SV  1  lycomponent  of  poiitlon  error 


SV  2  code  loop  error 


V  2  tropotpherlc  error 


V  2  lonoipherlc  error 


SV  2  clock  error _ 

SV  2  x-component  of  poiitlon  error 


SV  2  y-component  of  poiitlon  error 
SV  2* i* component  of  poiitlon  error 


SV  3  code  loop  error 


SV  3  tropoipherlc  error 


SV  3  lonoipherlc  error 


8R(7rx) 


hraiMi 


Sms  +  Smts  +  29 


Sms  +  Suits  +  30 


SV  3  y-component  of  poiitlon  error 


SV  3  *- component  of  poiitlon  error 


SV  4  code  loop  error 


SV  4  tropoipherlc  error 


SV  4  lonoipherlc  error 


SV  4  clock  error 


SV  4  x-compunent  of  poiitlon  error 


SV  4  y-component  of  poiitlon  error 


SV  4  z-component  of  poiitlon  error 


Appendix  B.  Litton  LN-93  Error-State  Model  Dynamics  Matrix 


The  LN-93  error-itate  dynamlci  matrix  (F)  ai  provided  by  Litton  li  a  93-by-93 
array  that  containi  a  large  number  of  dementi  that  are  Identically  lero.  Litton  partition! 
the  F  matrix  into  thirty-iix  lubarrayi  [20]  reflecting  the  logical  diviiiom  of  error  murcei 
dlicuued  in  Chapter  III. 

The  reader  ihould  note  that  only  the  NON-ZERO  element!  are  included  In  the  table! 
which  follow,  and  ihould  further  note  that  the  reviled  baro-altlmetcr  model  itatei  are  NOT 
Included  in  thli  let  of  ORIGINAL  F  matrix  elementi  extracted  from  the  Litton  document 
[20]. 

A  notationai  convention  [13]  it  to  label  elementi  of  the  C|  ,  lenior-to-true,  matrix 
ai  Cij  where  i  li  the  row  and  ii  the  column  in  the  traniformation  matrix. 
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Table  B.2.  Element!  of  the  Dynamic*  Submatrix  Pu  [20, 13] 


Element  Term 


Smmm 

■ded 


■OSDBEMqtDMKagi 


MlKfflillMESlIlMBM 

— WSmwmMMmMmMmm 

i«  1 1 1 


Table  B.3.  Elements  of  the  Dynamics  Submatrlx  Fki  [20,  13] 


B™ 


ement 


m  hiMB 

mmmwBmmmi 


Element 


wwmmi 


imimm 


c 


CviVfi, 


O.oCiv 


Cn 


nilHCKSl 

^aam 


~C713u^ 


C  i  ,iu»; 


0.5(7, 


•Cnw,b, 


lm 

mmm 

in 

S 

'■CEDI 

■GSQll 


-C'jnuijfc 


C*a:*u»j 


0.5Caii 


C;,m 


Cfoiuij 


.5  C:|I| 


Table  B,4.  Elements  of  the  Dynamic*  Submatrix  Fu  [20,  13] 


m 


Element  Term 


7,48)  |  Ci, 


KEMCnOMM 

■ULW^HIK 


ement  I  Term  ||  Element  Term 


Cu 


CyiA 


CialA 


cviAr 


III 


uA?A 


jjl 

)W0£Hi! 


8,83)  |  CnA 


8,56 


CvAjtAV 


lESEfftMWI 

■ 

n 

mm: 

in 


Cn\A 


8,88 


8,61)  CjiAtlAi 


B,64)  |  Ca;i j4 


mmi 


Cm  A?  A 


jEaiMQCB»Emoeiei(3nMB3BMi 

■■!  II I  !■■  lilMiM 


CXi  A 


Cm  A" A" 


Cm  A"  A 


Cm  A 


Table  B.5.  Elements  of  the  Dynamics  Submatrix  Fi#  [20,  13] 


II 


Element  Term 


Element  Term 


■ms 


Cvi 


8,71)  |  Cn 


Table  B.6.  Element!  of  the  Dynamics  Submatrix  Fi0  [20, 13] 


Table  B.7.  Element*  of  the  Dynamics  Submatrix  Fan  [20,  13] 


Element 

Term 

■ftTfcl 

eebi 

mmm 

warn 

KGXDB 

mm 

E02ED1 

"A. 

Ran 

mam 

mmm 

ggi 

BM 

mm 

m 

wm 

9ESa22KS 

iM 

■Hi 

■BBS1 

Table  B.8.  Elements  of  the  Dynamics  Submatrix  Fns  [20,  13] 


i  mm 

Eiszai 

mm 

(70,70 T 

Bn 

1020)1 

mm 

Eggp 

gag 

mmm 

BHEfflB 

BOH 

mmm 

ES 
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B.1  Elements  of  the  Process  Noise  Motrin 


The  Litton  document  [20]  includei  *  93-by-93  process  noiee  matrix  (Q)  for  the  LN-93 
error  model,  Like  the  7  matrix,  the  Q  matrix  i*  partitioned  into  eubarrayi  which  corre¬ 
spond  to  the  error-state  aubvecton  discussed  in  Chapter  III.  The  vast  majority  of  the 
elements  in  Q  are  identically  eero.  Only  the  non-zero  elements  of  Q  are  shown  below. 

Table  B,9.  Non-seero  Elements  of  Process  Noise  Submatrix  Qn  [20, 13] 


Element 

Term 

Element 

Term 

I  mssm 

b-m 

■QM 

EH! 

IKEH 

mmm 

EM! 

m 

EH 

IK2U 

EHI 

Table  B.10.  Non-sero  Elements  of  Process  Noise  Submatrix  [20,  13] 


Element 

Term 

Element 

Term 

(14,14) 

(15,15) 

2K*L 

(18,18) 

20b.„<. 

(17,17) 

(18,18) 

Ti»,»r 

2/^,X.„ 

(20,20) 

2/3«K*flr*Kr 

(21,21) 

2^*V<r« Ku 

""(22,22)’ 

"723-^T 

2  faliMh, 
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Appendix  C.  98-State  INS  Alignment  and  Flight  Simulations 


All  plot i  contained  in  this  and  subsequent  appendices  are  discussed  in  Chapter  VI 
(Results).  With  the  exception  of  Appendix  F,  all  plots  contained  in  this  and  subsequent 
appendices  contain  five  traces.  The  innermost  trace  (----)  on  each  data  plot  is  the  mean 
error  time  history  for  the  applicable  state  and  is  defined  by  [22]; 


&(<■•)  =  iE  «>(*,)  *  i  E  {aj(ti)  -  W<,)}  (c.i) 

Jy  iy 

where  «/(<,)  is  the  filter-computed  estimate  of  variable  i  and  aWijfai)  1»  the  truth  model 
value  of  the  same  variable,  at  time  for  sample  j,  and  N  is  the  number  of  time  histories 
in  the  simulation  (10  in  this  thesis). 

In  addition  to  the  center  trace,  two  more  pairs  of  traces  are  plotted.  The  first  pair 
(represented  by  •  •  •)  is  symmetrically  displaced  about  the  mean  and  as  a  result  follows 
the  “undulations"  of  the  The  locus  of  these  traces  is  calculated  from  [23]  A/„(t,)  ± 

\/FF(tTj ,  where  PE(ti)  is  the  true  error  covariance  at  time  t;.  The  true  standard  deviation 
is  calculated  from  [22]: 


<r„-M  =-  JPE{t7)  =  (C.2) 


where  N  is  the  number  of  runs  in  the  Monte  Carlo  simulation  (10  in  this  thesis),  and  Af,?(t,) 
is  the  mean-squared  value  of  the  variable  at  each  time  of  interest  (such  as  measurement 
times). 

The  last  pair  of  traces  ( — )  represents  the  filter  computed  ±  ,  values  for  the 

same  variables  of  interest  and  a/e  symmetrically  displaced  about  zero  because  the  filter 
“believes"  that  it  is  producing  zero  mean  errors  [24].  These  quantities  are  propagated  and 
updated  in  the  MSOFE  [8]  software  using  the  covariance  propagation  equation  shown  in 
Chapter  II.  These  traces  represent  the  filter’s  estimate  of  its  own  error. 
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C,  I  !)H-  State  Model;  Alignment  Using  Litton  Initial  Conditions 

The  plots  In  this  section  represent  results  of  a  10-run  Monte  Carlo  alignment  simula¬ 
tion,  In  these  alignments,  INS  aiding  consists  of  velocity  measurements  plus  baro-altimeter 
measurements. 

Note  that  the  plot  for  state  13,  A5.|  is  zero  for  all  time  during  all  alignment  simula¬ 
tions.  This  is  a  normal  condition  due  to  the  variable’s  dependence  on  altitude  rate.  T.he 
state  becomes  non-zero  during  flight  runs.  All  other  plots  are  discussed  in  Chapter  VI 
(Results). 

The  filter  computed  error  estimates  [07, are  compared  to  similar  plots  -ontained 
in  the  Litton  reference  documentation  [20].  Comparisons  are  qualitative  only  and  are 
intended  to  demonstrate  trends.  Such  comparisons  are  contained  in  Chapter  VI  (Results), 
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Figure  C.l.  Alignment;  93-State  Model  (a)  Latitude  and  (b)  Longitude  Errors. 
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Figure  C.4,  Alignment:  93-State  Model  (a)  North  Velocity  and  (b)  Vertical  Velocity 
Error  Statei. 
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Figure  C.5.  Alignment:  93-State  Model  (a)  INS  Altitude  Error  State  and  (H)  Baro- 
Altimeter  Total  Error. 
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Figure  C.8.  Alignment:  83-State  Model  (a)  i  and  (b)  A5|  Vertical  Channel  Aiding 
Error  State*. 
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C.S  93  State  Model)  Fighter  Flight  Uting  Litton  Initial  Condition* 

The  plot i  in  this  section  represent  reiulti  of  a  10-run  flight  simulation  in  which  a 
flight  proftle  (m  described  in  Chapter  III)  li  uied  to  character!**  the  LN-93  performance 
for  a  typical  fighter  mieiion.  (A  diacusiion  of  the  Litton  flight  profile  ii  provided  in  Chapter 
III.)  For  this  group  of  rune,  baro-altlmeter  aiding  it  the  only  meaeurement  need. 

The  purpose  of  this  eet  of  rune  li  to  establiah  that  the  software  implementation 
is  consistent  with  the  Litton  error  model  and  that  performance  result*  are  similar  to 
those  obtained  by  Litton.  Consistency  in  these  respects  does  not  constitute  error  model 
validation,  but  does  lend  confidence  that  the  software  functions  a*  intended  by  Litton. 

The  Alter  computed  error  estimates  [o/mc(.]  sure  compared  to  similar  plots  contained 
in  the  Litton  reference  documentation  [20].  Comparisons  are  qualitative  only  and  are 
intended  to  demonstrate  trends.  Such  comparisons  sue  contained  in  Chapter  VI  (Results). 
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Figure  C.7.  Flight:  93-State  Model  (a)  Latitude  and  (b)  Longitude  Error  States. 
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Figure  C.8.  Flight:  93-State  Model  (a)  East  Tilt  and  (b)  North  Tilt  Error  States. 
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Figure  C.tO.  Flight:  93-State  Model  (a)  North  Velocity  and  (b)  Vertical  Velocity  Error 
State*. 
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Figure  C.ll.  Flight!  93-State  Model  (a)  INS  Altitude  Error  State  and  (b)  Baro- Altimeter 
Total  Error. 
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Figure  C.12.  Flight:  93-State  Model  (a)  A$:j  and  (b)  AS 1  Vertical  Channel  Aiding  iSrror 
States. 
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Appendix  D.  96-State  INS  Alignment  and  Flight  Simulations 


ThJi  appendix  contains  the  plot*  'rotn  two  gets  of  10-run  Monte  Carlo  simulations 
of  the  96-state  error  model  which  incorporates  the  revised  baro-altimeter  model.  Plots 
contained  in  this  appendix  are  defined  in  exactly  the  same  manner  as  discu«sed  in  Ap¬ 
pendix  C. 

DJ  9 6 -State  Model:  Alignment  Using  Litton  Initial  Conditions 

The  plots  in  this  section  represent  results  of  a  10-run  Monte  Carlo  alignment  simula¬ 
tion.  In  these  alignments,  INS  aiding  consists  of  velocity  measurements  plus  baro-altimeter 
measurements.  Initial  conditions  are  again  those  chosen  by  Litton. 

Note  that  the  plot  for  state  13,  A5.|  is  zero  for  all  time  during  all  alignment  simula¬ 
tions.  This  is  a  normal  condition  due  to  the  variable’s  dependence  on  altitude  rate.  The 
state  becomes  non-zero  during  flight  rims.  All  other  plots  are  discussed  in  Chapter  VI 
(Results). 

The  filter  computed  error  estimates  [07 ,•/(„,.]  are  compared  to  similar  plots  contained 
in  the  Litton  reference  documentation  [20j-  Comparisons  are  qualitative  only  and  are 
intended  to  demonstrate  trends.  Such  comparisons  are  contained  in  Chapter  VI  (Results). 
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Figure  D.l.  Alignment:  96-State  Model  (a)  Latitude  and  (b)  Longitude  Errors.  (Note: 
Initial  conditions  are  those  chosen  by  Litton.) 
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Figure  D.2.  Alignments  96-State  Model  (a)  East  Tilt  and  (b)  North  Tilt  Error  States. 
(Notes  Initial  conditions  are  thou  chosen  by  Litton.) 


Figure  D.5.  Alignment!  96-State  Model  (a)  INS  Altitude  Error  State  and  (b)  Baro- 
Altimeter  Total  Error,  (Note:  Initial  condition*  are  thoie  choien  by  Litton.) 


D,t  96  Staff  hfodali  Fightar  Flight  Using  Litton  Initial  Conditions 

The  ploti  iu  thii  Mellon  represent  results  of  a  10-run  flight  simulation  In  which  a 
flight  profile  (u  described  in  Chapter  IU)  is  used  to  characterise  the  LN-83  performance 
for  a  typical  fighter  mission.  For  this  group  of  runs,  only  baro-altlmeter  aiding  is  used, 
Note,  however,  that  a  revised  baro-altimeter  model  is  in  place  in  the  truth  model  (see 
discussion  in  Chapter  III),  which  accounts  for  the  increase  of  3  states  in  the  overall  error 
model. 

The  purpose  of  this  set  of  runs  is.  to  establish  that  the  software  function  is  not 
adversely  affected  by  the  addition  of  the  new  baro-altimeter  states. 

The  filter  computed  error  estimates  [o/uinr I  are  compared  to  similar  plots  contained 
in  the  Litton  reference  documentation  [20].  Comparisons  sue  qualitative  only  and  are 
intended  to  demonstrate  trends.  Such  comparisons  are  contained  In  Chapter  VI  (Results). 
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Figure  D.9.  FUght:  98-St*te  Model  (a)  Aaimuth  and  (b)  East  Velocity  Error  State*. 
(Noter  Initial  condition*  are  tho*<*  chosen  by  Litton.) 
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Figure  D.10.  Flight:  06-State  Model  (a)  North  Velocity  and  (b)  Vertical  Velocity  Error 
State*.  (Note:  Initial  condition*  are  those  chosen  by  Litton.) 
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Figure  Dll,  FUghti  96-State  Model  (u)  INS  Altitude  Error  State  and  (b)  Bnro- Altimeter 
Total  Error,  (Note:  Initial  conditions  are  those  chosen  by  Litton.) 


“  -  *  • 

Mean  Error  =  A/.r  -  (M.,.)(n,„ 

Mean  Error  ±  c r/,.,,. 

— 

0  ±  ff/iltrr 

D-13 


Figure  D.12.  Flight:  98-Sta,te  Model  (a)  AS: i  and  (h)  AS i  Vertical  Channel  Aiding  Error 
States.  (Note:  Initial  conditions  are  those  chosen  by  Litton.) 
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Appendix  E.  96~State  INS  Alignment  and  Flight  Simulations:  Part  II 


Thii  appendix,  contains  the  plots  from  two  sets  of  10-run  Monte  Carlo  simulations  of 
the  96-state  error  model  which  incorporates  the  revised  baro- altimeter  model.  It  differs 
from  Appendix  D  in  that  initial  conditions  are  chosen  to  reflect  an  alignment  at  Holloman 
AFB,  NM.  This  step  is  taken  in  order  to  establish  a  performance  baseline  using  the  INS 
with  baro-aiding  only.  Subsequent  tests  with  ground  transponders  (see  Appendix  F)  and 
CPS  (see  Appendix  G)  are  compared  to  results  from  simulations  in  this  set. 

Variables  plotted  in  this  appendix  are  defined  in  exactly  the  same  manner  as  those  in 
Appendix  C. 

E,1  96-State  Model:  Alignment  Using  Holloman  Initial  Conditions 

The  plots  In  this  section  represent  results  of  a  10-run  Monte  Carlo  alignment  simula¬ 
tion.  In  these  alignments,  INS  aiding  consists  of  velocity  measurements  plus  baro-altlmcter 
measurements. 

Note  that  the  plot  for  state  13,  AA'.|  is  zero  for  all  timo  during  all  alignment  simula¬ 
tions.  This  is  a  normal  condition  due  to  the  variable's  dependence  on  altitude  rate.  The 
state  becomes  non-zero  during  flight  runs.  All  other  plots  are  discussed  in  Chapter  VI 
(Results). 

The  Alter  computed  error  estimates  [cry, are  also  compared  to  similar  plots  con¬ 
tained  in  the  Litton  reference  documentation  [20].  Comparisons  are  qualitative  only  and 
are  intended  to  demonstrate  trends.  Such  comparisons  are  contained  in  Chapter  VI  (Re¬ 
sults). 
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Figure  E.l.  Alignment:  96-State  Model  (a)  Latitude  and  (b)  Longitude  Errors.  (Note: 
Initial  conditions  are  those  for  Holloman  AF3.) 
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Figure  E.2.  Alignment:  96-State  Model  (n)  East  Till  and  (1>)  North  Tilt,  Error  States. 
(Note:  Initial  conditions  are  those  for  Holloman  AFB.) 
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Figure  E.3.  Alignment!  96-State  Model  (a)  Azimuth  and  (b)  East  Velocity  Error  States. 
(Note:  Initial  onditions  are  those  for  Holloman  AFB.) 
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Figure  E.4.  Alignment:  96-State  Model  (a)  North  Velocity  and  (l>)  Vertical  Velocity 
Error  State*.  (Note:  Initial  conditions  are  those  for  Holloman  AFB.) 
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Figure  E.6.  Alignment:  98-State  Model  (a)  A5- 1  and  (b)  AS 1  Vertical  Channel  Aiding 
Error  State*.  (Note:  Initial  conditions  are  those  for  Holloman  AFB.) 


E.S  96-State  Model!  Fighter  Flight  Using  Holloman  Initial  Conditions 

The  plot*  in  this  lection  represent  result*  of  a  10-run  flight  simulation  In  which  a 
flight  profile  (as  described  in  Chapter  III)  is  used  to  characterise  the  LN-93  performance 
for  a  fighter  mission  originating  and  terminating  at  Holloman  AFB,  NM.  For  this  group  of 
runs,  only  baro-aitimeter  aiding  is  used. 

The  purpose  of  this  set  of  runs  is  to  establish  that  the  software  function  for  a  flight 
using  Holloman  coordinates  results  in  performance  which  is  comparable  to  that  achieved 
in  the  Litton  flight  runs  reported  in  Appendix  C. 

The  filter  computed  error  estimates  [tr/,'f(,.r]  are  also  compared  to  similar  plots  con¬ 
tained  in  the  Litton  reference  documentation  [20].  Comparisons  are  qualitative  only  and 
are  intended  to  demonstrate  trends.  Such  comparisons  are  contained  in  Chapter  VI  (Re¬ 
sults). 
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Figure  E.8.  Flight:  96-State  Model  (a)  East  Tilt  and  (l»)  North  Tilt  Error  States.  (Note: 
Initial  conditions  are  those  for  Holloman  AFB.) 
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Figure  E.9.  Flight:  96-State  Model  (a)  Azimuth  nttd  (b)  East-  Velocity  Error  Stales. 
(Note:  Initial  conditions  are  those  for  Holloman  AFB.) 
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Figure  E.10.  Flight!  98-State  Model  (a)  North  Velocity  and  (b)  Vertical  Velocity  Error 
States,  (Note!  Initial  conditions  are  those  for  Holloman  AFB.) 
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Figure  E.ll.  Flight:  98-State  Model  (a)  INS  Alt  itude  Error  State  and  (b)  Bnro- Altimeter 
Total  Error.  (Note:  Initial  conditions  are  those  for  Holloman  AFB.) 
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Figure  E.12.  Flight:  96-State  Model  (a)  tlS:\  and  (b)  AS 1  Vertical  Channel  Aiding  Error 
States.  (Note:  Initial  conditions  are  those  for  Holloman  AFB.) 
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Appendix  F.  72-State  Reduced  INS  Model  Validation 


Thi*  appendix  contains  the  plots  from  a  10-run  Monte  Carlo  alignment  simulation  of  a 
72-state  error  model  which  incorporates  a  reduced  INS  model  and  a  revised  baro-altimeter 
model.  The  72-state  truth  model  one  sigma  values  obtained  during  thiB  simulation  are 
compared  directly  to  the  06-state  truth  model  one  sigma  values  obtained  during  the  align¬ 
ment  runs  presented  in  Appendix  E.  The  purpose  of  this  comparison  is  to  ensure  that  the 
72  state  model  adequately  represents  the  96  state  model, 

Prior  to  the  error  model  reduction  discussed  above,  truth  models  consisting  of  first 
122  and  finally  162  states  had  been  programmed  into  MSOFE.  The  122-state  ver  ion 
(which  integrated  the  96-state  INS  and  the  26-state  RRS  models)  was  successfully  com¬ 
piled  and  run.  Subsequently,  the  GPS  error  model  was  integrated  Into  the  overall  error 
model,  bringing  the  truth  and  Kalman  filter  error  model  dimensions  to  162  states  each. 
At  that  point,  hardware/ compiler  limitations  prevented  compiling  MSOFE  to  run  on  the 
accelerator  board.  Inadequate  random  access  memory  (RAM)  is  believed  to  be  the  cause 
of  difficulty.  In  order  to  continue  to  use  the  accelerator  board,  the  size  of  the  problem  had 
to  be  reduced. 

As  a  result,  several  of  the  less  important  states  contained  in  the  LN-93  error  model  are 
eliminated.  Certain  states  are  chosen  for  elimination  based  on  previous  research  [11,  19], 
and  validation  of  the  selection  is  accomplished  in  a  series  of  Monte  Carlo  runs  which 
culminate  in  the  plots  presented  in  this  appendix. 

Although  24  states  have  been  eliminated  from  the  INS  error  model,  the  plots  con¬ 
tained  in  this  appendix  demonstrate  that  no  significant  loss  of  fidelity  results.  In  all 
important  states  (i.e. ,  position,  platform  tilts,  and  velocities),  the  reduced  system  model 
accurately  portrays  the  error  behavior  of  the  full  S)6-st.ate  truth  model. 

Variables  plotted  in  this  appendix  are  simply  the  standard  deviation  time  histories  for 
several  truth  model  states  of  interest. 
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F  l  7S~State  MOM:  Alignment  Using  Holloman  Initial  Conditions 

The  plot . in  thl,  section  represent  result,  of  a  10-run  Monte  Carlo  alignment  simu¬ 
lation  of  the  reduced  INS  error  model.  In  these  alignments,  INS  aiding  consists  of  velocity 
measurement#  plus  baro-altimeter  measurements  only. 

The  hlter  computed  error  estimate,  [<ry,7(nr]  a«  «1,0  compared  to  similar  plot,  con- 
tained  in  the  Litton  reference  documentation  [20].  Comparison,  are  qualitative  only  and 
are  Intended  to  demonstrate  trends,  Such  comparisons  are  contained  In  Chapter  VI  (Re- 
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Figure  F.2.  Alignment:  98-State  vs.  72-State  Model  (a)  East  Tilt  and  (b)  North  Tilt 
Error  States.  (Note:  Initial  conditions  are  those  for  Holloman  AFB.) 
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Figure  F.6.  Alignment:  96-State  v»,  72-State  Model  (a)  INS  Altitude  Error  State  and 
(b)  Baro- Altimeter  Total  Error.  (Note:  Initial  condition*  are  those  for  Hol¬ 
loman  AFB.) 
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Appendix  G,  98-State  INS-RRS  Alignment  and  Flight  Simulations 


Thii  appendix  contain*  the  plot*  one  10- run  Monte  Carlo  alignment  and  one  5-run 
Monte  Carlo  flight  simulation  of  the  98-state  error  model  which  Incorporates  the  truncated 
LN-93  INS,  the  revised  baro-altimeter  model,  and  ERS.  Initial  conditions  are  again  chosen 
to  reflect  an  alignment  at  Holloman  AFB,  NM.  This  step  Is  important  because  actual  sur¬ 
veyed  positions  for  the  RRS  transponders  are  used  in  this  study  (the  RRS  transponders  are 
physically  located  on  or  near  Hollomsm  AFB).  These  results  (using  ground  transponders) 
are  compared  to  results  from  simulation*  In  Appendix  F  and  Appendix  H. 

Variables  plotted  In  this  appendix  are  defined  exactly  in  the  same  manner  as  those  in  Ap¬ 
pendix  C;  all  statistical  quantities  are  calculated  in  the  manner  described  at  the  beginning 
of  Appendix  C. 

G.l  98-State  Model:  Alignment  Using  Holloman  Initial  Conditions 

The  plots  In  this  section  represent  results  of  a  10-run  Monte  Carlo  alignment  simu¬ 
lation.  In  these  alignments,  INS  aiding  consists  of  velocity  measurements,  baro-altimeter 
measurements,  and  RRS  transponder  range  measurements. 

The  purpose  of  this  set  of  runs  is  to  establish  that  the  software  function  for  the 
INS-RRS  combination  is  comparable  to  (or  better  than)  that  achieved  in  the  Litton  flight 
runs  reported  in  Appendices  C  through  E. 

Note  that  the  plot  for  state  13,  AS. i  is  zero  for  all  time  during  all  alignment  simula¬ 
tions.  This  is  a  normal  condition  due  to  the  variable’s  dependence  on  altitude  rate.  The 
state  becomes  non-zero  during  flight  runs.  All  other  plots  are  discussed  in  Chapter  VI 
(Results). 

The  filter  computed  error  estimates  )  arc*  also  compared  to  similar  plots  con¬ 
tained  in  the  Litton  reference  documentation  [20].  Comparisons  are  qualitative  only  and 
are  intended  to  demonstrate  trends.  Such  comparisons  are  contained  in  Chapter  VI  (Re¬ 
sults). 
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Figure  G,3.  Alignment:  98-State  Model  (a)  Azimuth  and  (b)  Ea»t  Velocity  Error  States. 
(Note:  Initial  conditions  are  those  for  Holloman  AFB.) 
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Figure  G.4.  Alignment:  98-State  Model  (a)  North  Velocity  and  (b)  Vertical  Velocity 
Error  Statei.  (Note:  Initial  conditions  are  those  for  Holloman  AFB.) 
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Figure  CJ. 5,  Alignment!  98-State  Model  (a)  INS  Altitude  Error  State  and  (b)  Baro- 
Altimeter  Total  Error.  (Note:  lirltial  condition*  are  those  for  Holloman 
AFB.) 
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Figure  G.8.  Alignment:  98-State  Model  (a)  Transponder  1,  Z  Axis  Position  (b) 
Transponder  1,  Atmospheric  Propagation  Error  States 
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G.g  98-State  Model:  Fighter  Flight  Uev  g  Holloman  Initial  Condition! 

The  plots  In  this  section  represent  results  of  a  10-run  flight  simulation  in  which  a 
flight  profile  (as  described  in  Chapter  III)  is  used  to  characterise  the  performance  of  the 
LN-03with  RRS  aiding  for  a  fighter  mission  originating  and  terminating  at  Holloman  AFB, 
NM. 

The  purpose  of  this  set  of  runs  is  to  establish  that  the  software  function  for  the 
INS-RRS  combination  is  comparable  to  (or  better  than)  that  achieved  in  the  Litton  flight 
runs  reported  in  Appendices  C  through  E.  For  this  group  of  runs,  baro* altimeter  aiding 
'.a  used  in  conjunction  with  RRS  aiding. 

The  Alter  computed  error  estimates  [ffjiUtr]  are  also  compared  to  similar  plots  con¬ 
tained  in  tiie  Litton  reference  documentation  [20].  Comparisons  are  qualitative  only  and 
are  Intended  to  demonstrate  trends.  Such  comparisons  are  contained  in  Chapter  Vl  (Re¬ 
sults). 
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Figure  G.9.  Flight:  68-State  Model  (a)  Latitude  unci  (b)  Longitude  Error  State*.  (Note: 
Initial  condition*  are  thoie  for  Holloman  AFB.) 
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Figure  G.10.  Flight:  98-State  Model  (a)  East  Tilt  and  (b)  North  Tilt  Error  States,  (Note: 
Initial  conditions  are  those  for  Holloman  AFB.) 
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FlgureG.il.  Flight:  98-State  Model  (a)  Azimuth  and  (b)  Ea«t  Velocity  Error  Staten. 
(Note:  Initial  conditions  are  those  for  Holloman  AFB.) 
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Figure  Q.12.  Flight:  98-State  Model  (a)  North  Velocity  and  (b)  Vertical  Velocity  Error 
States.  (Note:  Initial  conditions  are  those  for  Holloman  AFB.) 
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Figure  G.13.  Flight:  98-State  Model  (a)  INS  Alt  itude  Error  State  and  (b)  Baro- Altimeter 
Total  Error.  (Note:  Initial  conditions  are  those  for  Holloman  AFB.) 
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Figure  G.14.  Flight:  98-State  Model  (a)  A5.i  and  (b)  A5|  Vertical  Channel  Aiding  Error 
States.  (Note:  Initial  conditions  are  those  for  Holloman  AFB.) 
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Figure  Q. 15.  Alignment:  98-State  Model  (a)  Transponder  1,  X  Axis  Position  (b) 
Transponder  1,  Y  Axis  Position  Error  States 


Figure  G.lfl.  Alignment:  08-State  Model  (a)  Transponder  1,  Z  Axis  Position  (b) 
Transponder  1,  Atmospheric  Propagation  Error  States 
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Appendix  H.  128-State  INS-GPS  Alignment  and  Flight  Simulation/) 


This  appendix  contains  the  plots  from  a  10-run  Monte  Carlo  alignment  and  a  3- 
run  Monte  Carlo  flight  simulation  of  the  128-state  error  model  which  incorporates  the 
truncated  LN-93  INS,  the  revised  baro-altimeter  model,  RRS,  and  GPS.  Initial  conditions 
are  once  again  chosen  for  an  alignment  at  Holloman  AFB,  NM.  These  results  using  GPS 
measurements  only  are  compared  to  results  from  simulations  In  Appendices  F,  G,  and  I. 
RRS  (transponder)  measurements  are  NOT  used  in  this  configuration. 

Variables  plotted  in  this  appendix  are  defined  exacciy  in  the  same  manner  as  those 
in  Appendix  C.  All  statistical  quantities  are  calculated  in  the  manner  described  at  the 
beginning  of  Appendix  C. 

H,1  128-State  Model:  Alignment  Using  Holloman  Initial  Conditions 

The  plots  in  this  section  represent  results  of  a  10-run  Monte  Carlo  alignment  simula¬ 
tion.  In  these  alignments,  INS  aiding  consists  of  velocity  measurements,  baro-altimeter 
measurements,  and  GPS  pseudo-range  measurements.  (RRS  measurements  are  NOT 
used.) 

The  purpose  of  this  set  of  runs  is  to  establish  that  the  software  function  for  the 
INS-GPS  combination  is  comparable  to  (or  better  than)  that  achieved  in  the  Litton  flight 
runs  reported  in  Appendices  C  through  F. 

Note  that  the  plot  for  state  13,  AS{  is  zero  for  all  time  during  all  alignment  simula¬ 
tions.  This  is  a  normal  condition  due  to  the  variable’s  dependence  on  altitude  rate.  The 
state  becomes  non-zero  during  flight  runs.  The  GPS  User  clock  states  exhibit  large  mag¬ 
nitude  transients  during  the  initial  phase  of  the  alignment  runs.  As  a  result,  an  additional 
“window"  is  plotted  on  the  page  which  follows  the  first  user  clock  plots.  The  purpose  is 
to  demonstrate  the  steady  state  behavior  of  these  critical  error  states.  All  other  plots  are 
discussed  in  Chapter  VI  (Results). 

The  filter  computed  error  estimates  [<r/, are  also  compared  to  similar  plots  con¬ 
tained  in  the  Litton  reference  documentation  [20].  Comparisons  are  qualitative  only  and 
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arc  Intended  to  demonstrate  trends.  Such  comparisons  are  contained  in  Chapter  VI  (Re¬ 
sults). 
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Figure  II.l.  Alignment:  128-State  Model  (a)  Latitude  and  (b)  Longitude  Errors.  (Note: 
GPS  Measurements  Only) 
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Figure  H.2.  Alignment:  128-State  Model  (a)  East  Tilt  and  (b)  North  Tilt  Error  States. 
(Note:  GPS  Measurements  Only) 
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Figure  H.3.  Alignment:  128-State  Model  (a)  Azimuth  and  (b)  East  Velocity  Error  States. 
(Note:  GPS  Measurements  Only) 
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Figure  H.4.  Alignment;  128-State  Model  (a)  North  Velocity  and  (b)  Vertical  Velocity 
Error  States.  (Note:  GPS  Measurements  Only) 
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Figure  II. 5.  Alignment:  128-State  Mode!  (a)  INS  Altitude  Error  State  and  (It)  Bnro- 
Altlmeter  Total  Error.  (Note:  GPS  Measurement*  Only) 


.... 

—  .  —  ■-  - - - JW*; - —j-  —  - —  ■■ 

Mean  Error  -  M.r  -  (M.,)ir,„- 

Mean  Error  ±-  erlnil 

— 

0  ±  ft J liter 

H-7 


(12)  DEL  S3  Error 
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Figure  H.6.  Alignment!  128-State  Model  (a)  d.S’i  and  (b)  A5|  Vertical  Channel  Aiding 
Error  Stater.  (Note!  GPS  Measurement h  Only) 
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Figure  H. 7.  Alignment!  128-State  Model  (a)  TrRtmponder  1,  X  Axil  Poiition  (b) 
Trawponder  1,  Y  Axil  Poiition  Frror  State*  (Note:  GPS  MeasurcmoiitH 
Only) 
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Figure  H. 8.  Alignment:  128-State  Model  (a)  Transponder  1,  Z  Axis  Position  (b) 
Transponder  1,  Atmospheric  Propagation  Error  States  (Note:  GPS  Mea¬ 
surements  Only) 
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FigureH.il.  Alignment:  128-State  Model  (a)  SV  1,  Codeloop  and  (b)  SV  1,  Atmospheric 
Propagation  Error  Statei  (Note:  (IPS  Measurement*  Only) 
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Figure  11.12.  Alignment!  128-State  Model  (a)  SV  I,  Clock  and  (b)  SV  1,  X  Position 
Error  States  (Note:  GPS  Measurements  Only) 
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11,2  128-State  Model:  Fighter  Flight  Using  Holloman  Initial  Condition s 

The  plots  in  this  lection  repreient  remits  of  a  3-run  flight  simulation  in  which  a 
flight  profile  (at  described  in  Chapter  III)  Is  used  to  characterise  the  performance  of  the 
LN-93with  RRS  aiding  for  a  fighter  mission  originating  and  terminating  at  Holloman  AFB, 
NM. 

It  is  well  understood  that  a  3-run  simulation  has  reduced  statistical  validity  [22,  23, 
24],  However,  the  eweeution  time  for  a  10-run  simulation  of  this  site  is  prohibitive.  There¬ 
fore,  the  3-run  results  are  included  for  a  CAUTIOUS  comparison. 

The  purpose  of  this  set  of  runs  is  to  establish  that  the  software  function  for  the 
INS-RRS-GPS  combination  is  comparable  to  (or  better  than)  that  achieved  in  the  Litton 
flight  runs  reported  in  Appendices  C  through  F. 

The  filter  computed  error  estimates  [cr/m, , .]  are  also  compared  to  similar  plots  con¬ 
tained  in  the  Litton  reference  documentation  (20).  Comparisons  are  qualitative  only  ami 
are  Intended  to  demonstrate  trends.  Such  comparisons  are  contained  in  Chapter  VI  (Re¬ 
sults). 
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Figure  H.14.  Flight:  128-State  Model  (a)  East  Tilt,  and  (b)  North  Tilt  Error  States. 
(Note:  GPS  Measurements  Only) 
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Figure  11,17.  Flight:  128-State  Model  (a)  INS  Altitude  Error  State  and  (b)  Daro 
Altimeter  Total  Error.  (Note:  GPS  Measurements  Only) 
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Figure  H.20.  Alignment:  128-State  Model  (a)  Transponder  1,  Z  Axl»  Position  (b) 
Transponder  1,  Atmospheric  Propagation  Error  States  (Note:  GPS  Mea¬ 
surements  Only) 


IMBBa 

Mean  Error  ~  A/,  (A/, 

Mean  Error  ±  trh 

IB 

0  ±  ff  filhr 

H-23 


0  1000  2000  3000  4000 


Tlm«  (uc) 

(b) 

Figure  H.21.  Flight:  128-State  Model  (a)  User  Clk  Bla*  find  (b)  Uner  Clork  Drift  Staton 
(Note:  GPS  Meaeuremcnts  Only) 
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Figure  H.22.  Flight!  128-State  Model  (a)  SV  I,  Cndeloop  and  (b)  SV  1,  Atmospheric 
Propagation  Error  States  (Note:  CiPS  Measurements  Only) 
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Figure  H.23.  Flight:  128-State  Model  (a)  SV  I,  Clock  and  (h)  SV  t,  X  Position  Error 
State*  (Note:  GPS  Muaiurement*  Only) 
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Appendix  I.  1B8-State  INS-RRS-GPS  Alignment  and  Flight  Simulations 

This  Appendix  contain*  the  plot*  from  a  10-run  Monte  Carlo  alignment  and  a  10- 
run  Monte  Carlo  flight  limulation  of  the  128-state  error  model  which  incorporate*  the 
truncated  LN-03  INS,  the  reviled  baro-altimeter  model,  RRS,  and  GPS.  Initial  condition! 
are  once  again  choien  for  an  alignment  at  Holloman  AFB,  NM.  These  results  (using  GPS 
measurements)  are  compared  to  results  from  simulations  in  Appendix  P,  Appendix  G,  and 
Appendix  H. 

Variables  plotted  in  this  appendix  are  defined  exactly  in  the  same  manner  as  those  in  Ap¬ 
pendix  C;  all  statistical  quantities  are  calculated  in  the  manner  described  at  the  beginning 
of  Appendix  C. 

LI  lS8-State  Model t  Alignment  Using  Holloman  Initial  Conditions 

The  plots  In  this  section  represent  results  of  a  10-run  Monte  Carlo  alignment  simu¬ 
lation.  In  these  alignments,  INS  aiding  consists  of  velocity  measurements,  baro-altimeter 
measurements,  RRS  transponder  range  measurements,  and  GPS  pseudo-range  measure¬ 
ments. 

The  purpose  of  this  set  of  runs  is  to  establish  that  the  software  function  for  the 
INS-RRS-GPS  combination  is  comparable  to  (or  better  than)  that  achieved  in  the  Litton 
flight  runs  reported  in  Appendices  C  through  F.  In  addition,  results  are  expected  to  be 
somewhat  better  than  those  reported  in  Appendix  G  and  Appendix  H. 

Note  that  the  plot  for  state  13,  A5.|  is  sero  for  all  time  during  all  alignment  simula¬ 
tions.  This  is  a  normal  condition  due  to  the  variable’s  dependence  on  altitude  rate,  The 
state  becomes  non-sero  during  flight  runs.  The  GPS  User  clock  states  exhibit  large  mag¬ 
nitude  transients  during  the  Initial  phase  of  the  alignment  runs.  As  a  result,  an  additional 
“window”  is  plotted  on  the  page  which  follows  the  firsl  user  clock  plots.  The  purpose  is 
to  demonstrate  the  steady  state  behavior  of  these  critical  error  states.  All  other  plots  are 
discussed  in  Chapter  VI  (Results), 

The  filter  computed  error  estimate*  [oj,un)  are  also  compared  to  similar  plots  con¬ 
tained  in  the  Litton  reference  documentation  [20].  Comparisons  are  qualitative  only  and 
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are  Intended  to  demomtrate  trend*.  Such  comparisons  are  contained  in  Chaptor  VI  (Re¬ 
sults). 
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Figure  1.1.  Alignment:  128 -State  Model  (a)  Latitude  and  (b)  Longitude  Errors.  (Note: 
Both  GPS  and  RRS  Measurements) 
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Figure  1.2.  Alignment;  128-State  Model  (a)  East  Tilt  and  (b)  North  Tilt  Error  States. 
(Note:  Both  GPS  and  RES  Measurements) 
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Figure  1.3.  Alignment:  128-State  Model  (a)  Azimuth  and  (b)  East  Velocity  Error  States. 
(Note:  Both  GPS  and  RRS  Measurements) 
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Figure  1.4.  Alignment:  128-State  Model  (a)  North  Velocity  anti  (b)  Vertical  Velocity 
Error  State*.  (Note:  Both  GPS  and  RRS  measurement*  used.) 
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Figure  1.6.  Alignment:  128-State  Model  (a)  INS  Altitude  Error  State  and  (b)  Baro- 
Altimeter  Total  Error,  (Note:  Both  GPS  and  RRS  Mearurements) 
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Figure  1.6.  Alignment:  128-State  Model  (a)  and  (b)  A£|  Vertical  Channel  Aiding 
Error  Statea.  (Note:  Both  GPS  and  RR.S  Measurements) 
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Figure  1.7.  Alignment:  128-State  Model  (a)  Transponder  1,  X  Axle  Position  (b) 
Transponder  1,  Y  Axle  Position  Error  Spates  (Note:  Both  GPS  and  HHS 
Measurement*) 
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Figure  1.8.  Alignment;  128-State  Model  (a)  Traniponder  1,  Z  Axi*  Portion  (b) 
Transponder  1,  Atmospheric  Propagation  Error  States  (Note:  Both  GPS  and 
RES  Measurements) 
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Figure  1.9.  Alignment:  128-State  Model  («)  User  C'lk  Bias  and  (b)  User  Clock  Drift 
States  (Note:  Both  GPS  Measurements  Only) 
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Figure  1. 10.  Alignment:  128-St«te  Model:  Steady  State  (a)  Uier  Clk  Dias  and  (b)  User 
Clock  Drift  States  (Note:  GPS  ami  HRS  Measurements) 
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Figure  1.11.  Alignment:  128— State  Model  (a)  S  V  1,  Codeloop  and  (b)  SV  1,  Atmoiphcrlc 
Propagation  Error  State*  (Note:  Both  GPS  Meaeurementi  Only) 
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Figure  1.12.  Alignment:  128-State  Model  (a)  SV  l,  Clock  and  (b)  SV  l,  X  Position  Error 
State*  (Note:  Both  GPS  Measurement)!  Only) 
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I.S  ISA -Slate  Model:  Fighter  Flight  Using  Holloman  Initial  Conditions 

The  plot i  In  this  lection  represent  results  of  &  10-run  flight  simulation  in  which  a 
flight  profile  (ai  deicribed  in  Chapter  UI)  i«  used  to  characterise  the  performance  of  the 
LN-93wlth  RRS  aiding  for  a  lighter  mission  originating  and  terminating  at  Holloman  AFB, 
NM. 

The  purpose  of  this  set  of  runs  is  to  establish  that  the  software  function  for  the 
INS-RRS-GPS  combination  le  comparable  to  (or  better  than)  that  achieved  in  the  Litton 
flight  runs  reported  in  Appendices  C  through  F.  In  addition,  results  are  expected  to  be 
somewhat  better  than  thole  reported  in  Appendix  Q. 

The  Alter  computed  error  estimates  [(T/,/*,,,.]  are  also  compared  to  similar  plots  con¬ 
tained  in  the  Litton  reference  documentation  [20].  Comparisons  are  qualitative  only  and 
ore  Intended  to  demonstrate  trends.  Such  comparisons  are  contained  in  Chapter  VI  (Re¬ 
sults). 
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Figure  1.15.  Flight:  128-State  Model  (a)  Azimuth  and  (1>)  East  Velocity  Error  States. 
(Note:  Both  GPS  and  RRS  Measurements) 
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(9)  U  Velocity  Error  (fps)  (8)  N  Velocity  Error  ftps) 
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Figure  1.19.  Alignment;  128-State  Model  (a)  Transponder  1,  X  Axis  Position  (b) 
Transponder  1,  Y  Axis  Position  Error  States  (Note:  Both  GPS  and  RRS 
Measurements) 
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Figure  1.20.  Alignment:  128-State  Model  (a)  Transponder  1,  Z  Axil  Position  (b) 
Transponder  1,  Atmospheric  Propagation  Error  States  (Note:  Both  GPS 
and  RRS  Measurements) 
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Figure  1.23.  Alignment!  128-State  Model  (a)  SV  1,  Clock  and  (b)  SV  1,  X  Position  Error 
Statei  (Notes  RRS  and  Both  GPS  Measurements) 
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